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ABSTRACT 


A  computerized  forest  inventory  system  has  been  developed  to  meet  the 
needs  of  New  Hampshire  landovmers  and  foresters.   The  system  will  analyse 
all  sampling  designs  currently  in  wide  use  throughout  the  State  producing 
tables  easily  understood  by  landowners.   The  system  includes  a  rigorous 
statistical  analysis  of  volume  standard  errors  as  well  as  a  newly  developed 
height  double  sampling  technique.   The  system  is  available  through  the 
Institute  of  Natural  and  Environmental  Resources. 


KEY  WORDS:   Forest  inventory.   Computer  inventory.   Sampling.   Double  sampling. 
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INVENT: 
FOREST  INVENTORY  SYSTEM  FOR  NEW  HAMPSHIRE  LANDOWNERS 

by 
D.  S.  Linden  and  J.  P.  Barrett* 

INTRODUCTION 
INVENT  Is  a  computerized  forest  inventory  processing  system  designed 
for  New  Hampshire  landowners.   INVENT  was  designed  with  the  following  five 
goals  in  mind: 

1.  Produce  a  system  that  will  analyse  all  sampling  designs  and  methods 
presently  in  wide  use  throughout  the  State. 

2.  Produce  output  tables  easily  understood  by  foresters  and  landowners — 
output  tables  that  look  as  if  they  were  designed  by  a  forester,  not 

a  computer. 

3.  Find  and  implement  volume  equations  that  are  well  adapted  to  the 
New  Hampshire  species — calculate  both  board-foot  and  cubic-foot 
volumes . 

A.   Indicate  the  reliability  of  volume  estimates. 

5.   Produce  a  system  that  any  forester  can  use  and  understand — a  system 
available  to  all  foresters  at  low  cost. 

The  realization  of  these  goals  is  the  result  of  two  years  work  in  close  coopera- 
tion with  consulting  foresters,  county  foresters,  and  forest  industry.  As  a 
result.  New  Hampshire  foresters  have  "been  emancipated  from  the  drudgery  and 
limitations  of  computation  and  can  concentrate  on  forest  data  needs."   (Davis 
1966,  p.  282). 

The  use  of  INVENT  is  made  available  by  the  Institute  of  Natural  and 
Environmental  Resources  at  the  University  of  New  Hampshire.   Institute  personnel 
will  process  data  from  tally  sheets  submitted  by  subscribing  foresters.  Persons 
interested  in  using  the  system  should  contact  Professor  James  Barrett  at  I.N.E.R., 
University  of  New  Hampshire. 

This  publication  is  written  in  four  separate  sections.   Section  1  contains 
all  the  information  a  field  forester  needs  to  know  to  utilize  the  service. 
Section  2  describes  how  the  data  is  computer  coded  and  how  the  program  is  run. 
Section  3  documents  the  computer  program  itself  with  subroutine  descriptions 


*Graduate  Research  Assistant  and  Professor  of  Forest  Resources,  Institute  of 
Natural  and  Environmental  Resources,  respectively. 
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and  variable  dictionaries.   Section  4  discusses  the  statistical  techniques 
used  within  the  program.   The  appendices  include  tally  forms,  species  code 
list,  volume  equations,  a  guide  to  sample  size  determination,  a  program 
listing  and  other  supportive  material. 
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SECTION  1.   FIELD  FORESTER  INFORMATION 

The  first  step  in  a  forest  inventory  is  to  systematically  arrange  the 
forest  into  groups  or  categories  based  upon  some  definite  scheme.  We  shall 
refer  to  these  groups  as  either  compartments  or  strata,  using  these  terms 
interchangeably.   The  boundaries  of  these  compartments  may  be  determined  by 
forest  type,  stocking  level,  merchantability  class  or  administrative  conveni- 
ence.  Whatever  the  criteria  actually  used,  INVENT  will  assume  that  each 
compartment  should  be  processed  individually  and  will  produce  a  complete  set 
of  summary  tables  for  each  compartment.   Furthermore,  INVENT  will  produce  a 
set  of  summary  tables  for  all  compartments  combined  assuming  stratified 
random  sampling  (Cf.  Freese  1962).   INVENT  will  handle  up  to  999  compartments 
although  in  practice  one  would  seldom  have  more  than  20  compartments  within 
any  forest. 

For  each  forest  to  be  processed,  the  forester  must  supply  the  informa- 
tion requested  on  the  INVENT  Processing  Request  Form  (Appendix  A,  page  48), 
This  information  includes : 

1.  Name  and  address  of  forester  making  request. 

2.  Number  of  compartments  in  the  forest. 

3.  Alternate  form  class  specification. 

4.  Diameter  class  size. 

5.  Confidence  level. 

6.  Tract  description- 
Alternate^  Form  Class  Specification.   If  heights  are  measured  to  the 

merchantable  top,  form  class  equations  are  used  to  determine  board  foot 
volume.   INVENT  uses  the  following  form  classes  as  a  default: 

Hemlock  70         All  other  softwoods  76         All  hardwoods  74 
Alternate  form  classes  may  be  specified  using  the  space  provided  on  the  INVENT 
Processing  Request  Form.   These  form  classes  will  be  used  for  all  compartments 
in  the  forest. 
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Dlameter  Class  Size.   INVENT  assumes  that  trees  will  be  grouped  Into 
diameter  classes.  The  size  of  the  class  Is  up  to  the  user  within  the  follow- 
ing limitations.   The  class  size  must  be  a  whole  number.  The  smallest  diameter 
class  allowable  is  the  class  size;  i.e.,  if  the  class  size  is  2  in.  then  the 
first  class  must  be  2  inches.  All  other  diameter  classes  must  be  integer 
DHjltiples  of  the  class  size.  The  maximum  diameter  class  must  be  less  than  or 
equal  to  40  inches. 

Confidence  Level.   INVENT  uses  the  95  percent  confidence  level  as  a 
default.   Using  this  level,  there  is  only  a  one  in  twenty  chance  that  the 
true  volume  will  be  outside  the  confidence  interval  computed  by  INVENT  (Cf. 
Freese  1967,  p.  11).  This  is  the  most  common  confidence  level  used  in 
forestry  in  the  Northeast.   If  the  user  wishes  to  use  a  different  level,  he 
may  specify  that  level  as  long  as  it  is  in  the  range  from  65  percent  to 
99  percent. 

Tract  Description.  A  tract  description  of  up  to  70  characters  should  be 
specified  for  the  forest.  This  description  will  be  printed  at  the  top  of 
each  page  of  the  combined  compartment  summaries. 

For  each  compartment  within  the  forest,  the  forester  must  supply  the 
information  requested  on  the  INVENT  Compartment  Information  Form  (Appendix  A, 
page  49).  This  information  includes: 

1.  Compartment  name  and/or  description. 

2.  Sampling  method  and  sample  unit  size. 

3.  Tree  height  measurement  units. 

4.  Tree  top  specification. 

5.  Whether  or  not  the  heights  of  all  volume  trees  were 
measured. 

6.  Size  of  compartment  in  acres. 

Compartment  Name.   A  compartment  name  and/or  description  should  be 
specified  for  each  compartment.  The  description  may  be  up  to  70  characters 
in  length  and  is  printed  at  the  top  of  each  page  of  the  compartment  output 
tables. 
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Sampling  Method.   INVENT  will  process  point  cruises  (variable  plot) , 
plot  cruises  (fixed  radius  plot),  strip  cruises,  and  100  percent  tallies. 
The  sampling  method  need  not  be  the  same  for  each  compartment.   When  using 
point  or  plot  sampling,  a  complete  statistical  analysis  of  the  reliability 
of  the  volume  estimates  will  be  calculated  only  if  the  data  for  each  point 
or  plot  is  tallied  separately.   If  only  summary  data  is  available,  no  such 
statistical  analysis  can  be  made.  When  point  sampling,  the  BAF  of  the  angle 
used  must  be  specified.  When  plot  sampling,  the  plot  size  in  acres  must  be 
specified.   Strip  cruises  are  analysed  as  a  plot  cruise  where  only  one  large 
plot  is  sampled.   In  this  case,  the  plot  size  is  the  total  acreage  of  all 
strips.  There  is  no  statistical  analysis  of  reliability  on  strip  cruises. 
One  hundred  percent  tallies  require  no  statistical  analysis  of  reliability 
since  all  trees  in  the  population  have  been  measured.  One  hundred  percent 
tallies  are  also  treated  as  a  plot  cruise  where  only  one  large  plot  is  sampled. 
Here,  however,  the  plot  size  is  equal  to  the  total  acreage  of  the  compartment. 

Tree  Height  Measurement.   Tree  height  may  be  measured  in  feet  or  sixteen 
foot  logs  to  the  nearest  half -log. 

Tree  Top  Specification.  Trees  may  be  measured  to  a  merchantable  top  or 
measured  by  total  height.  There  is  also  a  mixed  option  where  softwoods  are 
measured  by  total  height  and  hardwoods  by  merchantable  height.  As  mentioned 
earlier,  form  class  volume  equations  for  board  foot  volume  are  used  when 
merchantable  height  is  specified. 

Were  Heights  of  All  Volume  Trees  Measured?   INVENT  allows  double  sampling 
for  height  (Cf.  Freese  1962,  p.  43).  Using  this  technique,  the  forester  only 
measures  height  on  a  subsample  of  the  trees.   INVENT  develops  a  height-dbh 
relationship  from  the  height  trees  sampled  and  uses  this  relationship  to 
estimate  height  for  the  remaining  trees.   This  method  must  be  used  with  some 
caution.   On  at  least  one-third  of  the  trees  in  every  species,  the  height 
should  be  measured.   At  least  three  height  trees  must  be  measured  for  every 
species  within  each  compartment.   An  analysis  table  of  the  height-dbh  relation- 
ship is  printed  for  each  compartment  using  the  option.  A  warning  message  is 
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printed  if  an  insufficient  niimber  of  trees  within  any  species  was  sampled. 
(See  table  17). 

Size  of  Compartment .   The  size  of  each  compartment  in  acres  must  be 
specified.   If  not  known  exactly  it  must  at  least  be  estimated. 

Tally  Forms  and  Product  Specification 

Tally  forms  are  available  for  use  with  INVENT  (see  Appendix  A,  page  46) . 
Trees  are  tallied  using  a  two  letter  species  code  (Appendix  A,  page  45), 
diameter  class,  height  and  product.   The  product  column  is  used  as  follows: 
INVENT  allows  trees  to  be  graded  as  sawlog  or  pulpwood.   If  you  do  not  wish 
to  grade  trees  simply  ignore  the  product  column  and  all  trees  will  be  assumed 
sawlog.   If  you  do  wish  to  grade  trees  there  are  two  methods  available.   You 
can  grade  the  entire  tree  as  sawlog  or  pulpwood  by  entering  an  S  in  the  product 
column  for  sawlog,  or  entering  a  P  for  pulpwood.  Alternatively,  you  may  spec- 
ify the  percent  sawlog  in  the  tree.   If  you  put  60  in  the  product  column,  60 
percent  of  the  volume  will  go  into  sawlog  and  40  percent  will  go  into  pulpwood. 

The  column  marked  #  is  used  if  two  or  more  trees  of  the  same  species,  size 
and  product  are  tallied  on  a  given  plot. 

Output  Tables 
Seventeen  tables  are  presented  as  an  example  of  INVENT 's  output.  Tables 
1-16  are  part  of  a  59  page  inventory  report  based  on  a  3  compartment  cruise 
of  Whitaker  Woods  in  North  Conway,  N.H.  Tables  1-9  are  for  compartment  1 
while  tables  10-16  are  for  all  3  compartments  combined.   Table  17  is  taken 
from  another  cruise  where  the  height  double  sampling  option  was  used.   This 
table  is  a  summary  of  that  relationship.   The  output  tables  are  further  dis- 
cussed in  section  3  on  page  34.  However,  the  tables  should  be  self-explanatory 
to  most  foresters. 
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Table  1.  Job  Control  Information  (Individual  Stratum) 


INVENT    VER. 


DiTE    OF    EON: 


11/1/78  I.N.E.R. 

7-Dec-78 


O.N.H. 


B.A.F.  =    10.00 

DIAHETER    CLASS    SIZE    =       2 

HEIGBTS    OF    ALL      VOLOHE    TREES    HEASOBED 

INPUT    FILE    NAME:     CONilT.CI 

CONFIDENCE    LEVEL    =    .95 


NOBBEE    OF    POINTS  SAHPLED    =      20 

HEIGHTS    HEASORED  BY    HERCH    HEIGHT    IN    LOGS 

STSATOB    I       1 

ACREAGE  OF  TRACT  =    a7.0 


Table  2.  Species  Stand  and  Stock  Table 


WHITE    PINE  FC    =    76    * 

* 
************************************************************************************************** 


BASAL    AREA 

TREES 

SAHLOG 

SAHLOG 

POLPHOOD 

DIAHETER 

PER 

PER 

COBIC    FOOT    VOLDME 

BOARD    FOOT    VOLOHE 

CUBIC    FOOT    VOLUHE 

CLASS 

ACRE 

ACRE 

PER    ACRE 

PER     ACRE 

PER    ACRE 

8 

1.5 

4.3 

0.0 

0.0 

30.7 

10 

1.0 

1.8 

16.8 

82.8 

0.0 

12 

«.5 

5.7 

91.1 

535.4 

0.0 

14 

5.0 

4.7 

109.9 

705.8 

0.0 

16 

7.5 

5.4 

168.7 

:           1115.0 

16.4 

18 

9.0 

5.1 

221.1 

1545.0 

9.6 

20 

7.0 

3.2 

196.5 

1405.3 

0.0 

22 

3.5 

1.3 

101.6 

747.0 

0.0 

2H 

2.5 

0.8 

70.3 

:              538.0 

0.0 

26 

a. 5 

1.2 

134.1 

1035.6 

0.0 

28 

1.0 

0.2 

33.7 

258.3 

0.0 

30 

1.0 

0.2 

28.4 

233.  1 

0.0 

32 

2.0 

0.4 

58.4 

483.9 

0.0 

TOTAL 

50.0 

34.3 

1230.6    ♦ 

35)1 

8685.2  i    35X 

56.7   1  83X 

HEAN  STAND  DIAHETER  =  16.3 


BOARD  FOOT  COEFFICIENT  OF  VARIATION  =   78X 


MERCHANTABLE  H.S.D.  =  16.3 


POLP-HOOD   COEFFICIENT  OF  VARIATION  =  183% 


PERCENT  CRUISE  =   6.2H 
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Table  3,  Species  Stand  and  Stock  Table 


*****•••*•**•***•••••••*•*•••***••****••**••***•••••••••••**•«•••*••***•••*•••**••*********••*•**** 


• 

BED 

OAK 

FC   =    74    » 

• 

• 

*•*********•*•**•****•*****************•«••*•••*•**•**••••*«•*•••*••**•**»********•**•*•******•**•* 

:       EASiL   AfiEA 

TBEES 

SAWLOG 

SAHLOG 

PULPBOOD 

DliHETEB 

PES 

PEB 

CUBIC    FOOT    yOLDHE 

BOABD    FOOT    VOLOHE 

CUBIC    FOOT    VOLOHE 

CL&SS 

ACRE 

ACSE 

PEB    ACBE 

PEB    ACBE 

PEB    ACBE 

6 

1.0 

5.1 

0-0 

0.0 

19.4 

8 

1.5 

14.3 

6.1 

26.3 

16.7 

10 

1.5 

8.3 

5U.  1 

:               273.4 

13.2 

12 

5.5 

7.0 

72.11 

417.1 

0.0 

14 

2.0 

1.9 

30.9 

183.9 

0.0 

16 

2.0 

l.t 

28.2 

179.2 

0.0 

18 

3.0 

1.7 

U1.8 

:              277.1 

0.0 

20 

1.5 

0.7 

25.1 

165.4 

0.0 

22 

0.5 

0.2 

U.6 

:                35.6 

0.0 

TOTAL 

21.5 

30.5 

263.3    * 

58X 

1557.9  1   56JI         :                49.3    ♦151X 

HEAN    STAND    DIADETEB    =    11.4 


BCAED    FOOT    COEFPICIEMT    OF    VABIATIOH    =    122* 


HEBCHAHTABLE    H.S.D.    =    12.2 


PDLP-BOOD      COEFFICIENT    OF    TABIATION    =    327% 


PEBCENT    CEOISE    =       3. OX 


Table  ^.     Softwood  Summary  Stand  and  Stock  Table 


«  AIL  SCFTBOODS 

* 


BASAL    ABEA 

TBEES 

SA8L0G 

SAMLOG 

POLPHOOD 

DIABETEB 

PEE 

PEB 

COBIC    FOOT    VOLUHE 

BOABD    FOOT    VOLOHE 

COBIC    FOOT    VOLOBE 

CLASS 

ACBE 

ACBE 

PEB    ACBE 

PEB    ACRE 

PER    ACRE 

6 

0.5 

2.5 

0.0 

0.0 

8.1 

8 

3.5 

10.0 

0.0 

0.0 

57.3 

10 

7.0 

12.8 

86.1 

392.6 

4.0 

12 

10.5 

13.4 

169.7 

9  56.7 

0.0 

14 

9.0 

8.4 

183.7 

1147.4 

0.0 

16 

11.5 

8.2 

238.4 

1574.8 

16.4 

18 

13.5 

7.6 

307.3 

2079.7 

9.6 

20 

8.5 

3.9 

226.6 

1593.2 

0.0 

22 

4.5 

:                 1.7 

125.5 

915.2 

0.0 

24 

2.5 

0.8 

70.3 

538.0 

0.0 

26 

:                   4.5 

1.2 

134.1 

1035.6 

0.0 

28 

:                   1.0 

0.2 

33.7 

258.3 

0.0 

30 

1.0 

:                  0.2 

28.4 

233.1 

0.0 

32 

2.0 

0.4 

58.4 

483.9 

0.0 

TOTAL 

:               79.5 

:                71.5 

1662.2  1 

29X 

11208. 3   ♦    30X 

95.4   ♦    63% 

SEAN  STAND  DIAHETEB  =    14.3 


BOARD  FOOT  COEFFICIENT  OF  VARIATION  =   65X 


BEBCBANTABLE  B.S.D.  =  14.5 


PDLP-BOOD      COEFFICIENT    OP    VARIATION    =    138X 


PEBCENT    CROISE    =       4.7X 
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Table  5,  Hardwood  Summary  Stand  and  Stock  Table 


•  AIL    BIBDROODS  • 


:       BASAL    AREA 

TSEES 

SAflLOG 

SAHLOG 

POLPiOOD 

OI&HETEB 

PER 

PES 

COBIC    FOOT    VOLOHE 

BOARD    FOOT    VOLOHE 

COBIC    FOOT    TOIDBE 

CLASS 

ACRE 

ACRE 

PER    ACRE 

PEB    ACRE 

PER    ACBE 

6 

3.0 

15.3 

0.0 

0.0 

48.6 

8 

3.5 

10.0 

6.1 

26.3 

47.8 

10 

6.0 

11.0 

60.0 

302.8 

25.2 

12 

7.5 

9.5 

77.3 

448.  5 

14.4 

^n 

3.0 

2.8 

1*3.8 

260.4 

0.0 

16 

2.5 

1.6 

35.8 

224.4 

0.0 

18 

3.0 

1.7 

41.8 

277.  1 

0.0 

20 

1.5 

0.7 

25.1 

165.4 

0.0 

22 

0.5 

0.2 

4.6 

35.6 

0.0 

TOTAL 

30.5 

53.0 

29U.7  1    52% 

1740.5   +    5U 

136.0   1   72* 

HEAH    STAHD   DTAHETEB    =    10.3 


BOARD    FOOT   COEFFICIEBT   OF   TABIATIOB    =    109X 


HEBCBAHTABLE    H.S.D.    =    11.6 


PDLP-BOOD      COEFFICIEHT   OF   VABIATIOM    =    157X 


PEBCEMT    CBOISE   =      2.  4 J 


Table  6.  All  Species  Summary  Stand  and  Stock  Table 


ALL  SPECIES 


BASAL    AREA 

DIAHETEB 

PEE 

CLASS 

ACRE 

6 

3.5 

8 

7.0 

10 

13.0 

12 

18.0 

14 

12.0 

16 

14.0 

18 

16.5 

20 

10.0 

22 

5.0 

24 

:                 2.5 

26 

:                 4.5 

28 

1.0 

30 

:                 1.0 

32 

:                 2.0 

TBEES 

FEB 
ACRE 


SARLOG 

COBIC    FOOT    VOLOHE 

PER    ACRE 


SANLOG 

BOARD    FOOT    VOLOHE 

PEB    ACRE 


POLPROOD 

COBIC    FOOT    VOLOHE 

PEB    ACBE 


17.8 

0.0 

20.1 

6.1 

23.8 

146.1 

22.9 

246.9 

11.2 

227.6 

10.0 

274.3 

9.3 

349.1 

4.6 

251.7 

1.9 

130.1 

0.8 

70.3 

1.2 

134.1 

0.2 

33.7 

0.2 

28.4 

0.4 

58.4 

0.0 

26.3 

695.4 

1405.2 

1407.8 

1799.2 

2356.7 

1758.6 

950.8 

538.0 

1035.6 

258.3 

233.1 

483.9 


56.7 
105.1 

29.2 

14.4 
0.0 

16.4 
9.6 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


TOTAL 


110.0 


124.5 


1956.8    +    21)1 


12948.9   *    22X 


231.4   *    50% 


HEAN    STAHD   DIAHETER    =    12.7 


BOABD    FOOT   COEFFICIEHT   OF    VABIATIOH    =      48X 


HEBCBAHTABLE    H.S.D.    =    13.5 


POLP-ROOD      COEFFICIEHT  OF    VASIATIOH    =    109X 


PERCEHT  CRDISE  =   3. 8X 
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Table  7,  Species  Composition  Table 


•  * 

*  SPECIES    COHPOSITIOH    BT    PEBCEHT  * 


S&WLOG 

SAWLOG 

PULPWOOD 

COBIC    FOOT 

BOARD    FOOT 

:          COBIC    FOOT 

SPECIES 

BASAL    ABEA 

as. 45 

TREES 

VOL ORE 

VOLOHE 
67.07 

VOLOHE 

UHITE    PINE 

27.59 

62.89 

24.50 

RED    PINE 

5.91 

5.15 

6.90 

6.76 

0.00 

HEMLOCK 

20.91 

211.67 

15.15 

12.72 

16.75 

SDGAB    HiPLE 

0.91 

3.20 

0.00 

0.00 

6.82 

RED    BAPLE 

0.91 

1.53 

0.26 

0.25 

3.68 

WHITE    BIBCH 

0.91 

2,56 

0.25 

:                  0.24 

3.92 

BEECH 

5.U5 

10.80 

1.09 

0.91 

23.03 

RED   O&K 

19.55 

2«.51 

13.46 

12.03 

21.31 

ilL    SOFTBOODS 

72.27 

57. U1 

61.94 

86.56 

41.25 

ALL    HARDWOODS 

27.73 

U2.59 

:                15.06 

13.44 

58.75 

_„_—_« 

_________«___«. 







Table  8.  Volume  Summary  for  All  Species 


*************************************************************************************************** 

*  * 

*  VOLOHE  TOTALS  FOR  ALL  SPECIES  ♦ 

*  * 


SAWLOG 

COBIC    FOOT    VOLOHE 

PER    ACRE 


SAWLOG                                 : 

BOARD    FOOT    VOLOHE                : 

PEE    ACRE                              : 

8685. 

2                            : 

876 

0                            : 

1647. 

2                          : 

0 

0                          : 

32. 

8                            : 

31. 

4                         : 

118. 

4                          : 

1557. 

9                          ; 

PDLPWOOD 

COBIC    FOOT    VOLOHE 

PEE    ACRE 


WHITE    PINE 
RED    FINE 
H EH LOCK 
SOGAE    HAPLE 
RED    HAPLE 
WHITE    BIRCH 
BEECH 
RED    OAK 


1230.6 

135.0 

296.5 

0.0 

5.1 

4.8 

21.4 

263.3 


56.7 

0.0 

38.8 

15.8 

8.5 

9.1 

53.3 

49.3 


ALL    SOFTWOODS 
ALL   EABDWOODS 


1662.2 
294.7 


11208.3 
1740.5 


95.4 
136.0 


ALL    SPECIES 


1956.8 


12948.9 
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Table  9.  Volume  Summary  for  All  Species  Expanded  by  Acreage 


t^t ********************** ************************************************************************* 

* 

VOLUHE  TOTALS  EXPSNDED  BY  iCBEiGE  * 

* 

ttttt ********************************  ************************************************************* 


S&RLOG 
CUBIC  FOOT  VOLDBE 


SKHLOG 
BOARD  FOOT  VOLOBE 


POLPHOOD 
CUBIC  FOOT  TOLOSE 


WHITE  PIHE 
BED  PINE 
HEMLOCK 
SDGAB  SAPLE 
BED  HAPLE 
WHITE  BIBCH 
BEECH 
BED  OAK 


57838 

63a7 

13936 

0 

2«0 

228 

looa 

12377 


408202 

U1170 

77al9 

0 

1542 

1478 

5563 

73222 


2664 
0 

1822 
742 
400 
426 

2505 

2317 


ALL  SOFTWOODS 
ALL  HARDWOODS 


78122 
13849 


526792 
81805 


4486 
6390 


ALL  SPECIES 


91971 


608597 


10876 


Table  10.  Job  Control  Information  (Stratified  Total) 

IHVENT    VER.     3  11/1/78  I.H.E.E.  O.N.H. 

DATE    OF    BON:        7-Dec-78 


DIAHETEE    CLASS    SIZE    =       2 

TOTAL    ACEEAGE    SAMPLED    =  105.0 

C08FIDEBCE    LEVEL    =    .95 


NOHBEE    OF    STRATA    SAMPLED    =       3 
TOTAL    POINTS    SAMPLED    =  56 
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Table  11,  Species  Volume  Table  (Stratified  Total) 

*  » 

*  WHITE    PINE  PC    =    76    » 

*  * 


DIAHETES 
CL&SS 


SiHLOG 
CDBIC  FOOT  VOLOBE 


SAHLOG 
BOARD  FOOT  VOLOHE 


POLPWOOD 
COBIC  FOOT  VOLBHE 


6 

8 
10 
12 

ia 

16 
18 
20 
22 
2"t 
26 
28 
30 
32 
3U 
36 


0 

823 

1«99 

5111 

5550 

9600 

1«5ai 

141)68 

10177 

7902 

92117 

2506 

2074 

4228 

568 

455 


0 

2718 

7465 

30058 

35604 

63484 

100826 

103635 

75100 

60263 

71309 

19544 

16148 

34874 

4555 

3920 


914 
1716 
0 
0 
0 
769 
453 
0 
0 
0 
0 
0 
0 
0 
0 
0 


TOTAL 


88750  +  26X 


629503  ♦  26X 


3853  ♦  73* 


Table  12.  Species  Volume  Table  (Stratified  Total) 


BED  OAK  FC  =  74 


DIAHETER 
CLASS 


SAHLOG 
CDBIC  FOOT  VOLOHE 


SAHLOG 
BOARD  FOOT  VOLOHE 


POLPHOOD 
COBIC  FOOT  VOLOHE 


6 
8 
10 
12 
14 
16 
18 
20 
22 
24 
28 


0 

1234 

23270 

42632 

16995 

29345 

19598 

9308 

1674 

0 

2035 


1635 

2981 

622 

310 

0 

0 
664 
294 

0 
291 

0 


24968  ♦  37X 


146091  ♦  36X 


6799  ♦  61X 
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Table  13.  Softwood  Summary  Volume  Table  (Stratified  Total) 


•  * 

*  AIL    SCFTHOODS  • 


DI4METEE 
CLASS 


SAHLOG 
CDBIC  FOOT  VOLDIJE 


SAILOG 
BOAED  FOOT  VOLOHE 


PULPHOOD 
CDBIC  FOOT  VOLOBE 


6 
8 
10 
12 
14 
16 
18 
20 
22 
24 
26 
28 
30 
32 
34 
36 
40 


0 

1053 

6105 

10944 

12957 

22644 

26139 

18292 

13103 

9151 

10723 

2799 

3291 

4407 

915 

45S 

459 


0 

3451 

28274 

61345 

79177 

141792 

172099 

127282 

95876 

68560 

81395 

21620 

24339 

36297 

7039 

39  20 

3204 


1874 

4161 

320 

107 

0 

769 

591 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


TOTAL 


143437  +  19X 


955670  ♦  19X 


7822  +  50X 


Table  14.  Hardwood  Summary  Volume  Table  (Stratified  Total) 


AIL  EAROHOODS 
************************************************************************************************* 


DIASETEE 
CLASS 


SAHLOG 
COBIC  FOOT  VOLOHE 


SAWLOG 
BOARD  FOOT  VOLOHE 


PDLPWOOD 
COBIC  FOOT  VOLOHE 


6 
8 
10 
12 
14 
16 
18 
20 
22 
24 
26 
28 


607 

1234 

48605 

56781 

24463 

34023 

19598 

9308 

2825 

0 

2009 

2035 


7887 
11231 

2785 

2805 
826 
488 

1212 
«58 
610 
607 
237 
0 


35115  +  26X 


201490  +  26% 


29146  +  27X 
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Table  15,  All  Species  Summary  Volume  Table  (Stratified  Total) 


*  ilL  SPECIES  » 


DIAHETEB 
CL&SS 


SAiLOG 
COBIC  FOOT  VOLDHE 


SAWLOS 
BOAED  FOOT  VOLOHE 


POLPUOOD 
COBIC  FOOT  VOLOHE 


6 
8 
10 
12 
14 
16 
18 
20 
22 
24 
26 
28 
30 
32 
34 
36 
40 


366 

1341 

15522 

20955 

17124 

28116 

29144 

19695 

13484 

9151 

11038 

3089 

3291 

II407 

915 

455 

459 


607 

4685 

76879 

118127 

103641 

175815 

191697 

136589 

98701 

68560 

83404 

23655 

24339 

36297 

7039 

3920 

3204 


9760 

15393 

3104 

2912 

826 

1257 

1803 

458 

610 

607 

237 

0 

0 

0 

0 

0 

0 


178552  +  14)1 


1157160  ♦  14J 


36968  ♦  22X 


Table  16.  Volume  Summary  for  All  Species  (Stratified  Total) 


»  • 

»  VOLOBE  TOTALS  FOE  ALL  SPECIES  * 

«  * 


SPECIES 


SAHLOG 
COBIC  FOOT  VOLOME 


SARLOG 
BOARD  FOOT  VOLOHE 


PDLPHOOD 
COBIC  FOOT  VOLOHE 


OBITE  PINE 
BED  PINE 
HEMLOCK 

ED-BL-WH  SPEOCE 
SOGAB  MAPLE 
BED  HAPLE 
OBITE  BIBCH 
VELLOS  BIRCH 
BEECH 
BED  OAK 


88750 

10094 

44593 

0 

0 

1072 

4474 

89 

4513 

24968 


629503 

65651 

260516 

0 

0 

6271 

25213 

649 

23266 

146091 


3853 

0 

3746 

224 

742 

6150 

3255 

0 

12200 

6799 


AIL  SOFTWOODS 
ALL  BABDHOODS 


143437 
35115 


955670 
201490 


7822 
29146 


ALL  SPECIES 


178552 


1157160 


36968 
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Table  17.  Regression  Summary  (If  Option  Used) 


JIHILISIS  OF  THE  HEIGHT  -  DBH  RELiTIOHSHIP 
aODEL:   HEIGHT  =  BO  +  B1/0BH 


_ 

BO 

. 

B1 

. 

B 

STAHOABD 

ATEEAGE 

VOLOBE 

HEIGHT 

SPECIES 

IHTEECEPT    : 

SLOPE 

: 

SQOABED 

EBBOB    % 

HEIGHT 

TBEES 

TBEES 

8HITE   PIHE 

3.296 

-21.321 

':" 

54.0 

5.1 

1.5 

92 

40 

HEBLOCK 

1.799 

-8.528 

; 

16.2 

:             9.6 

1.1 

48 

23 

BED    BIPLE 

3-106 

-29.436 

: 

62.1 

13.5 

0.9 

39 

11 

BLICK    CBERBI 

•  «« 

HABHIliG 

*** 

0.0 

0.0 

3 

0 

SHEET    BIBCB 

•  *• 

RABHIBG 

•  *• 

:             0.0 

0.0 

1 

0 

B»SS»OOD 

0.0 

2.0 

2 

2 

BED   OilK 

1.1«0 

: 

-5.U71 

; 

30.1 

:           19.7 

0.8 

:                 5 

4 

WHITE   OAK 

*•* 

KABNIBG 

••• 

0.0 

0.0 

3 

0 

Bl&CK   OXK 

••* 

RABHIBG 

•  •• 

0.0 

0.0 

2 

0 

ELH 

1.157 

z 

-«.992 

; 

18.7 

:           15.0 

0.8 

9 

6 

BICKOBT 

2.865 

• 

-23.076 

• 

73.0 

9,2 

1.3 

6 

3 

ilL    SOPTBOODS 

:             a.6 

1.4 

140 

63 

ALL    HABOROODS 

s             8.1 

1.0 

70 

26 

ALL  SPECIES 


4.1 


1.3 


210 


89 


••*  RABHIBG  ••• 


IHSDFFICIEHT  BDHBEB  OF  HEIGHT  TBEES  BEASOBED  FOB  THIS  SPECIES! 
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SECTION  2,   DATA  CODING  AND  PROGRAM  EXECUTION 

Program  Operating  Environment 

INVENT  was  developed  on  the  DEC-10  system  at  the  University  of  New  Hamp- 
shire.  INVENT  is  an  interactive  system  acquiring  all  of  its  job  control 
related  information  from  the  user  through  a  dialogue.   Unit  5  (TTY  on  DEC-10) 
is  used  as  the  input  and  output  device  for  all  of  INVENT 's  interactive 
dialogue. 

For  each  stratum  the  plot  and  tree  data  must  reside  on  a  disk  file  with 
a  unique  name.  The  name  of  the  file  is  specified  by  the  user  at  execution 
time  as  INVENT  processes  each  stratum.   Unit  21  is  used  to  input  each  disk 
file.   The  file  uses  a  fixed  format  as  described  below. 

If  the  user  wishes  to  use  alternate  form  classes,  these  must  reside  on 
a  disk  file  also.   The  user  specifies  the  name  of  this  file  at  execution 
time  and  unit  1  is  used  to  input  this  file.   The  file  uses  a  fixed  format  as 
described  on  page  18. 

All  output  is  written  to  unit  3  (LPT  on  DEC-10)  to  a  file  named  INVENT. LPT. 

Data  File  Format  and  Coding 
An  individual  data  file  must  exist  for  each  stratum.  The  file  consists 
of  card  images  coded  as  follows: 

Card  Type    Columns Format Variable Description 

Title       1-70      12A5       TITLE       Stratum  title 

Plot         1-2         12        POINT       Plot  number  (1) 

Number  of  trees  on  plot  (1) 

Tree        1-2        A2       SPEC       Species  code 

Diameter  at  breast  height 

Product  code  (2) 

Tree  height 

Number  of  trees  of  this 

description  (3) 
Percent  sawlog  (4) 

(1)  When  using  the  summary  option  where  trees  are  not  tallied  by  Individual 
plots,  only  one  plot  card  is  used  and  it  is  coded  as  follows:   POINT =-9, 
TREES  =  number  of  plots  sampled . 

(2)  Product  codes:   Sawlog  (S  or  blank);  Pulpwood  (P) ;  Percent  Sawlog  (%) . 

(3)  If  IREP  is  coded  as  0  or  blank,  the  program  assumes  IREP  =  1. 

(4)  PCTSAW  is  ignored  unless  PROD  =  %.   PCTSAW  is  coded  as  a  percent  (0-100) 
and  the  F3.2  format  reads  it  as  a  decimal  (0-1.00). 


1-2 

12 

POINT 

3-5 

13 

TREES 

1-2 

A2 

SPEC 

4-5 

12 

DBH 

6 

Al 

PROD 

7-9 

F3. 

0 

HEIGHT 

10-11 

12 

IREP 

12-14 

F3. 

2 

PCTSAW 
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Examples  of  various  tree  card  codings  may  be  found  in  Appendix  A, 
page  ^^.   The  species  code  list  is  found  in  Appendix  A,  page  45. 

The  data  file  is  constructed  by  combining  the  title,  plot,  and  tree 
cards  as  follows.   The  first  card  in  the  file  is  the  title  card.   This  is 
followed  by  the  plot  card  for  plot  number  1.   The  plot  card  is  then  followed 
by  all  the  tree  cards  for  that  plot.   Next,  the  plot  card  for  plot  number  2 
is  included  followed  by  its  tree  cards;  and  so  on  for  the  remaining  plots. 
The  actual  value  of  POINT  (plot  number)  need  not  be  consecutive. 


An  example  of  an  INVENT  tally  sheet  and  the  resulting  data  file  listing  may 
be  found  in  Appendix  A,  pages  46  and  47. 

Alternate  Form  Class  File  Coding 
The  user  has  the  option  to  specify  a  set  of  alternate  form  classes. 
Form  class  only  affects  volume  when  tree  heights  are  measured  to  a  merchant- 
able top  (see  page  3) .  A  file  containing  these  alternate  form  classes  must 
reside  on  a  disk  file.   This  file  consists  of  four  rows  and  eight  columns 


organized  as  follows: 


Row  1 
Row  2 
Row  3 
Row  4 


Form  classes  for  species  1-8   (see  species  list  in 
"      "     "     "     9-16  Appendix  A,  page  45) 

17_24 
It       II      II      II     2S— '^2 


Each  row  is  formatted  8(12, X). 
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As  an  example,   assume  the  following  form  classes  were  desired: 

All  pines  80 

Hemlock,  Spruces,  Ash,  78 

Black  Cherry 

Balsam  Fir,  Tamarack,  Cedar 

Other  Softwoods ,  Sugar  Maple       76 

White  Birch,  Yellow  Birch 

All  other  species  75 

The  data  file  that  would  implement  these  form  classes  is: 

80  80  80  78  76  76  78  78 

80  76  76  76  75  78  75  75 

75  78  76  76  75  75  75  75 

75  75  75  75  75  75  75  75 

Running  the  Program 

Running  INVENT  consists  of  interactively  supplying  the  program  all  job 
control  information.  This  information  includes:  alternate  form  class 
specification  (if  desired),  diameter  class  size,  number  of  strata,  confidence 
level,  input  file  name  for  each  strata,  sampling  method  (point  or  plot), 
basal  area  factor  or  plot  size,  method  of  height  measurement,  acreage  of 
each  stratum,  and  title  of  combined  strata  output  (if  more  than  one  stratum 
sampled) . 

This  interactive  portion  of  the  program  is  self-explanatory.  The  program 
asks  the  user  certain  questions  while  listing  the  possible  answers.  The 
number  of  questions  asked  depends  on  the  options  used  and  the  number  of  strata 
sampled.   The  questions  are  asked  as  each  data  file  is  read  and  processed. 
Under  conditions  of  heavy  computer  usage  there  may  be  a  time  lapse  of  as  much 
as  two  minutes  between  certain  questions.   This  process  is  best  explained  by 
example. 

Tables  1-16  are  part  of  a  59  page  inventory  report  based  on  a  3  compart- 
ment cruise  of  Whitaker  Woods  in  North  Conway,  N.H.   Listed  below  is  the 
entire  user-program  dialogue  that  produced  that  output.   Each  of  the  questions 
asked  has  been  numbered  for  reference  and  user  responses  are  underlined. 
Most  questions  are  totally  self-explanatory  and  therefore  only  a  few  will  be 
discussed. 
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♦RUN  INUENT 


1,      DO  YOU  WANT  TO  READ  IN  ALTERNATE  FORM  CLASSES  FROM  DISK? 
NO 


2. 
3. 


6. 

7. 
8. 

9. 

10. 
11. 
12. 


ENTER  DIAMETER  CLASS  SIZE 
2 


ENTER  NUMBER  OF  STRATA 
3 


INVENT  USES  THE  95%  CONFIDENCE  LEVEL  AS  A  DEFAULT. 
DO  YOU  WISH  TO  USE  A  DIFFERENT  LEVEL? 


5      ENTER  INPUT  FILE  NAME  FOR  STRATUM  *    1 
CONWIT.Cl 


WHITAKER  WOODS   1978    COMP.  1 


SPECIFY  SAMPLING  METHOD 

TYPE  POINT  OR  PLOT 
POINT 


ENTER  B.A.F. 
10 


ARE  HEIGHTS  IN  FEET  OR  SIXTEEN  FT.  LOGS? 

TYPE  FEET  OR  LOGS 
LOGS 


ARE  HEIGHTS  TOTAL  OR  MERCHANTABLE? 

TYPE  TOTAL f  MERCH»  OR  MIXED 
MERCH 


WERE  THE  HEIGHTS  OF  ALL  VOLUME  TREES  MEASURED? 


ENTER  ACREAGE  OF  TRACT. 
42 


ENTER  INPUT  FILE  NAME  FOR  STRATUM  ♦    2 
C0NWIT.C2 


WHITAKER  WOODS   1978    COMP.  2 


13. 

14. 
15. 


16. 


-21- 


SPECIFY  SAMPLING  METHOD 

TYPE  POINT  OR  PLOT  OR  SAME 
SAME 


ENTER  ACREAGE  OF  TRACT. 
42 


ENTER  INPUT  FILE  NAME  FOR  STRATUM  *    3 
C0NUIT.C3 


UHITAKER  WOODS   1978    COMP.  3 


SPECIFY  SAMPLING  METHOD 

TYPE  POINT  OR  PLOT  OR  SAME 
SAME 


17  ENTER  ACREAGE  OF  TRACT, 

18  WHAT  DO  YOU  WANT  TO  TITLE  THE  COMBINED  STRATA  OUTPUT? 
WHITAKER  WOODS   1978    STRATIFIED  TOTAL    CQMPS.  1  -  3 

END  OF  EXECUTION 

CPU  time:  4.32   ELAPSED  TIME.*  5:31.23 

EXIT 

Answering  NO  to  questions  1  and  4  indicates  the  user  will  use  the  default 
form  classes  and  the  default  confidence  level  of  95  percent. 

Question  9  asks  the  user  if  heights  are  total  or  merchantable  and  offers 
three  responses:   TOTAL,  MERCH,  or  MIXED.   The  response  MIXED  would  have 
indicated  that  softwoods  were  measured  by  total  height  while  hardwoods  were 
measured  by  merchantable  height.   In  this  case,  however,  MERCH  was  answered 
indicating  merchantable  height  was  used  on  all  trees. 

There  would  be  a  pause  of  about  thirty  seconds  between  questions  11  and 
12  as  the  program  processes  stratum  1. 

Question  13  asks  the  user  to  specify  the  sampling  method  as  did  question 
6.  However,  question  13  offers  an  additional  response:   SAME.  Answering  SAME 
indicates  that  in  stratum  2  the  sampling  method,  basal  area  factor,  and  method 
of  height  measurement  are  identical  to  those  in  the  previous  stratum.   There- 
fore questions  7-10  are  not  repeated. 
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Once  again  there  would  be  a  pause  of  about  thirty  seconds  between 
questions  lA  and  15  and  between  questions  16  and  17. 

Question  18  asks  for  a  title  for  the  combined  strata  output.  This 
title  may  be  up  to  70  characters  in  length  and  is  printed  at  the  top  of  each 
table  in  the  combined  strata  output  section. 

After  an  additional  pause  of  about  thirty  seconds  the  end  of  execution 
message  appears  and  the  output  is  sent  to  the  line  printer  queue. 

The  following  example  is  Included  to  Illustrate  the  alternate  form 
class  option  and  the  alternate  confidence  level  option. 

♦RUN  INVENT 


1. 

lA. 
iB. 


2. 

3. 

4. 

4a. 


DO  YOU  WANT  TO  READ  IN  ALTERNATE  FORM  CLASSES  FROM  DISK? 
YES 


ENTER  FORM  CLASS  FILE  NAME 
UHIT.FRM 


DO  YOU  WANT  A  LIST  OF  THE  FORM  CLASSES  JUST  READ? 
YES 


UP  =80  RP  =  80  PP  =80  HM  =  78 

BF  =  76  TA  =  76  SP  =  78  NS  =  78 

SC  =80  CE  =  76  OS  =76  SM  =  76 

RM  =75  UA  =  78  WI  =75  AS  =  75 

YP  =  75  BC  =  78  WB  =  76  YB  =  76 

SB  =  75  GB  =  75  BE  =  75  BW  =  75 

RO  =  75  WO  =  75  BO  =75  EL  =  75 

(3U  =  75  HI  =  75  HH  =  75  OH  =  75 

ENTER  DIAMETER  CLASS  SIZE 


ENTER  NUMBER  OF  STRATA 
3 


INVENT  USES  THE  95%  CONFIDENCE  LEVEL  AS  A  DEFAULT. 

DO  YOU  WISH  TO  USE  A  DIFFERENT  LEVEL? 

YES 

ENTER  CONFIDENCE  LEVEL  AS  A  DECIMAL  IN  THE  RANGE  .65  TO  .99 

.90 
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Answering  YES  to  question  1  envokes  the  alternate  form  class  option. 
Question  la  asks  for  the  file  name  while  question  lb  gives  the  user  the 
option  to  list  the  new  form  classes.   These  form  classes  are  the  ones  de- 
scribed on  page  19.   Answering  YES  to  question  4  allows  the  user  to  specify 
an  alternate  confidence  level. 
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SECTION  3.   PROGRAM  DOCUMENTATION 

The  FORTRAN  program  consists  of  the  main  program:   INVENT;  nine  sub- 
programs:  FRMCLS,  LINREG,  GUTS,  VOLMER,  VOLTOT,  STRAT ,  OUTPl,  MDSTI,  and 
MDNRIS;  and  five  named  common  storage  areas:   COMUNE,  FORMCL,  REG,  OUTl, 
and  STOUT.  A  description  of  each  along  with  a  variable  dictionary  are 
presented  below.   The  program  listings  are  found  in  Appendix  E. 


Common  Area  COMUNE 

COMUNE  transfers  job  control  information  from  the  main  program  to  the 

subroutines. 

Description 

Acreage  of  compartment 

Prism  factor  or  plot  size 

Diameter  class  size 

Sampling  method 

Input  file  name 

Height  measurement  length  technique 

Conversion  factor  to  transform  logs 

to  feet 
Two-tailed  area  under  t-curve 
Regression  option  indicator 
Output  title  indicating  sampling 

method  used 
Compartment  title 
Height  measurement  top  technique 


Variabl( 

2  Type      D 
Real 

imension 

ACRES 

BAF 

Real 

- 

CLASIZ 

Int 

- 

FRAME 

Real(C*5)l 

- 

IFILE 

D.P.(C*10) 

- 

LOGANS 

Int  (C*4) 

- 

MULT 

Int 

- 

PROBLV 

Real 

_ 

REGANS 

Int  (c*3) 

- 

SAMMTD 

D.P.(C*8) 

- 

TITLE 

Int  (C*5) 

14 

TOPANS 

Int  (C*5 

- 

Common  Area  FORMCL 
FORMCL  transfers  the  form  class  array. 
Variable  Type      Dimension  Description 

FC      Int  32       Form  class  array 

Common  Area  REG 
REG  transfers  information  to  and  from  the  regression  subroutine. 


Description 

Average  height 

Minimum  diameter  sampled  for  height 
Number  of  height  trees  sampled 
Regression  slope 
Regression  Intercept 
Coefficient  of  determination 
Standard  error  of  the  regression 
Species  code 

Number  of  non-cull  trees  above 
minimum  d.b.h.  sampled 


■(C*5)  indicates  the  variable  is  used  to  store  5  characters. 


Variable 

I     Type 

Dimension 

MEANY 

Real 

35 

MINDBH 

Int 

32 

NTR 

Int 

35 

REGCOF 

Real 

33 

REGINT 

Real 

33 

RSQR 

Real 

33 

SEREG 

Real 

35 

SPECIE 

Int  (C*2) 

32 

TOTNTR 

Int 

35 
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Common  Area  OUTl 


OUTl  transfers  processed  data  to  the  output  subroutine. 


Variable 

.  Type 

Dimension 

BASPAC 

Real 

40,35 

BASTOT 

Real 

35 

BDTOT 

Real 

35 

BDVLSM 

Real 

40,35 

CBTOT 

Real 

35 

CBVLSM 

Real 

40,35 

FINPOP 

Real 

— 

NPNTS 

Int 

— 

PLTOT 

Real 

35 

PLVLSM 

Real 

35 

SEBRD 

Real 

35 

SECUB 

Real 

35 

SEPUL 

Real 

35 

TREPAC 

Real 

40,35 

TRETOT 

Real 

35 

VARBRD 

Real 

35 

VARPUL 

Real 

35 

Description 

Basal  area/acre  by  d.b.h. 

Basal  area/acre 

Sawlog  bd.ft.  volume 

Sawlog  bd.ft.  volume  by  d.b.h. 

Sawlog  cub.ft.  volume 

Sawlog  cub.ft.  volume  by  d.b.h. 

Finite  population  correction 

Number  of  sample  units  sampled 

Pulpwood  cub . ft .  volume 

Pulpwood  cub.ft.  volume  by  d.b.h. 

Sawlog  bd.ft.  volume  standard  error 

Sawlog  cub.  ft.  volume  standard  error 

Pulpwood  cub.ft.  volume  standard  error 

Trees/acre  by  d.b.h. 

Trees/acre 

Sawlog  bd.ft.  volume  variance 

Pulpwood  cub . ft .  volume  variance 

Note;  All  volumes  are  per  acre 


Common  Area  STOUT 
STOUT  transfers  combined  strata  information  to  the  output  subroutine. 


Description 

Sawlog  bd.ft.  voltime 

Sawlog  bd.ft.  volume  by  d.b.h. 

Sawlog  cub.ft.  volume 

Sawlog  cub.ft.  volume  by  d.b.h. 

Sawlog  bd.ft.  volume  effective  d.f. 

Sawlog  cub.ft.  volume  effective  d.f. 

Pulpwood  cub.ft.  volume  effective  d.f. 

Pulpwood  cub.ft.  volume 

Pulpwood  cub.ft.  volume  by  d.b.h. 

Sawlog  bd.ft.  volume  standard  error 

Sawlog  cub.ft.  volume  standard  error 

Pulpwood  cub.ft.  volume  standard  error 

Total  acres  sampled 

Note :  All  volumes  are  total 


Variable 

Type 

Dimension 

BDTST 

Real 

35 

BDVLST 

Real 

40,35 

CBTST 

Real 

35 

CBVLST 

Real 

40,35 

EFDFBD 

Real 

35 

EFDFCB 

Real 

35 

EFDFPL 

Real 

35 

PLTST 

Real 

35 

PLVLST 

Real 

40,35 

STSEBD 

Real 

35 

STSECB 

Real 

35 

STSEPL 

Real 

35 

TACRES 

Real 

— 

Main  Program  INVENT 
INVENT  provides  interactive  user  communication  and  subroutine  control. 
Input  and  output  are  to  unit  5  which  is  the  user's  TTY  on  the  DEC-10.   The 
following  job  control  information  is  acquired  interactively  from  the  user: 

(1)  Alternate  form  class  specification. 

(2)  Diameter  class  size. 

(3)  Nvnnber  of  strata. 
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(4)  Confidence  level. 

(5)  Input  file  name. 

(6)  Sampling  method  and  sampling  frame. 

(7)  Method  of  height  measurement. 

(8)  Acreage  of  strata. 

(9)  Title  of  combined  strata  output  (if  more  than  one  strata) . 
INVENT  opens  the  output  file  on  unit  3  (LPT  on  DEC-10)  naming  it 

INVENT. LPT.   INVENT  also  sets  the  default  confidence  level  to  95  percent. 

The  following  subroutine  calls  are  made  in  INVENT:   FRMCLS,  LINREG, 
GUTS,  OUTPl,  0UTP2,  STRAT,  and  STRAT2. 

Common  Area  Used:   COMUNE. 

Other 

Variables  Type     Dimension  Description 

ALPHA  Real  —  Confidence  level  as  a  decimal 

ALPHCH  Real(C*3)  —  Confidence  level  change  flag 

FRMANS  Real(C*3)  —  Form  class  change  flag 

III  Int  —  Loop  counter  for  strata  number 

ISTRAT  Int  —  Number  of  strata 

REPEAT  D.P.(C*8)  ~  User  message 

SVFRAM  Real(C*5)  —  Saves  the  contents  of  FRAME 

Subroutine  FRMCLS 

FRMCLS  is  called  by  INVENT  if  the  user  specifies  that  he  wishes  to  use 
alternate  form  classes.   FRMCLS  interactively  asks  the  user  to  specify  the 
name  of  the  file  where  the  alternate  form  classes  are  stored.   This  file  is 
read  from  unit  1  (disk  on  DEC-10).   As  in  INVENT,  unit  5  is  used  to  communicate 
with  the  user.   The  user  is  given  the  option  to  list  the  new  form  classes  after 
they  are  read  from  the  file. 

Common  Areas  Used:   FORMCL  and  REG. 

Other 

Variables  Type     Dimension  Description 

FCFILE    D.P.(C*10)   —        Input  file  name 

I         Int         —        Implied  do  loop  index 

LSTANS    Int         —        Form  class  list  option  flag 

Subroutine  LINREG 
LINREG  is  called  by  INVENT  if  the  user  specifies  that  only  the  heights 
of  a  subsample  of  the  volume  trees  were  measured.   LINREG  then  reads  through 
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the  data  file  specified  by  the  variable  IFILE  and  located  on  unit  21  (disk 
on  DEC-10) .   For  every  species  with  three  or  more  height  trees  measured, 
LINREG  fits  the  following  regression  model:  HEIGHT  =  B  +  B  *  1/DBH.   The 
following  statistics  are  calculated  and  stored  in  common  area  REG: 

(1)  Regression  intercept. 

(2)  Regression  slope. 

(3)  Standard  error  of  the  regression. 

(4)  Number  of  height  trees  measured. 

(5)  Total  number  of  merchantable  trees  sampled. 

(6)  Average  height. 

2 

(7)  Coefficient  of  determination  (R  ) . 

(8)  Mlnimtim  diameter  of  all  height  trees  sampled  if  less  than  8  inches. 
LINREG  ignores  cull  trees  (HEIGHT  =  999) .   It  also  checks  that  input  parameters 
are  within  the  following  ranges : 

1  <_  DBH  <_  40     0  <_  HEIGHT  £  7  if  height  measured  by  logs. 

0  <_  HEIGHT  <^  200  if  height  measured  in  feet. 
If  data  is  found  outside  these  ranges,  the  value  is  set  to  the  minimum  or 
maximum  of  the  range.   If  an  illegal  species  code  appears,  the  tree  is  ignored 
and  an  error  message  is  written  to  unit  5. 

LINREG  also  calculates  certain  statistics  for  the  Softwood,  Hardwood  and 
All  Species  accumulations.  These  statistics  include  (3),  (4),  (5)  and  (6) 


above . 

Common 

Area  Used: 

REG 

, 

Other 

Variables 

Type 

Dimension 

CORMNY 

Real 

_. 

CORNTR 

Real 

— 

COVXY 

Real 

33 

DBH 

Int 

~ 

HEIGHT 

Real 

— - 

IDUMMY 

Int 

— - 

IREP 

Int 

— 

ISPEC 

Int 

— 

ITYPE 

Int 

— 

JTYPE 

Int 

— 

LOOP 

Int 

— 

MEANX 

Real 

35 

MEANY 

Real 

35 

MNYCOR 

Real 

3 

Description 

Corrected  mean  height  for  species 

Corrected  number  of  trees  sampled 

Covariance  of  X  and  Y 

Diameter  at  breast  height 

Tree  height 

DO  loop  index  for  reps  on  input 

Number  of  reps  on  input 

Species  loop  counter 

Forest  type  code 

Forest  type  code 

Loop  counter 

Mean  of  X 

Mean  of  Y 

Corrected  mean  height  for  forest  types 
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Other 


Variable 

'.s     Type 

Dimension 

Description 

MSE 

Real 

— 

Mean  square  error 

NTRCOR 

Int 

3 

Corrected  number  of  trees  for  forest  type 

NUMREG 

Int 

3 

Number  of  regressions  run 

POINT 

Int 

— 

Point  number 

SPEC 

Int(C*2) 

— 

Species  input  code 

SSE 

Real 

— 

Sum  of  the  squares  of  error 

TREES 

Int 

— 

Number  of  trees  on  a  sample  unit 

TYPSSE 

Real 

3 

Forest  type  error  sum  of  squares 

VARX 

Real 

33 

Variance  of  X 

VARY 

Real 

33 

Variance  of  Y 

X 

Real 

33 

1/D.B.H. 

XSUM 

Real 

33 

Sum  of  X 

XSUMSQ 

Real 

33 

Sum  of  X^ 

Y 

Real 

33 

Height 

YSUM 

Real 

35 

Sum  of  Y 

YSUMSQ 

Real 

33 

Sum  of  y2 

Subroutine  GUTS 

GUTS  is  called  by  INVENT  to  read  the  input  data  and  calculate  the  basal 
area  per  acre,  trees  per  acre,  and  all  of  the  volumes  along  with  their  stan- 
dard errors.   GUTS  also  calculates  the  number  of  units  sampled  and  the  percent 
cruise.   Either  subroutine  VOLMER  or  VOLTOT  is  called  by  GUTS  repeatedly  to 
calculate  tree  volume. 

GUTS  reads  the  input  file  specified  by  the  variable  IFILE  and  located  on 
unit  21  (disk  on  DEC-10) .   It  checks  the  range  of  input  data  in  the  same  way 
LINREG  does  except  that  the  minimum  allowed  DBH  is  set  at  the  value  of  CLASIZ 
Instead  of  1. 

GUTS  checks  the  value  of  POINT  on  the  first  plot  to  see  if  POINT  =  -9. 
Such  a  value  indicates  the  data  is  presented  in  summary  mode  (not  tallied  by 
individual  plot)  and  therefore  no  statistical  analysis  of  the  sampling  error 
can  be  made.   It  also  indicates  that  the  value  of  TREES  is  the  number  of  plots 
sampled . 

GUTS  reads  all  the  data  accumulating  sums  of  basal  areas,  number  of  trees, 
sawlog  board  foot  volume,  sawlog  cubic  foot  volume,  and  pulpwood  cubic  foot 
volume.   After  all  the  data  is  read,  the  accumulated  sums  are  expanded  into 
per  acre  estimates  based  on  the  sampling  technique  used  (point  or  plot)  .   The 
variances  and  standard  errors  of  these  volumes  are  also  calculated.   Before 
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attempting  to  calculate  variances,  GUTS  checks  to  see  if  more  than  one  plot 
has  been  sampled.   If  not,  GUTS  sets  itself  into  summary  mode  and  does  not 
attempt  to  calculate  any  variances  or  standard  errors. 


Common  Areas  Used:   COMUNE,  REG,  and  OUTl. 
Other 


Variables 

L  Type 

Dimension 

DBVOL 

Real 

^_ 

BRDPT 

Real 

35 

BRDSQ 

Real 

35 

CUBPT 

Real 

35 

CUBSQ 

Real 

35 

CUBVOL 

Real 

— 

DBH 

Int 

— 

FACTOR 

Real 

— 

HEIGHT 

Real 

— 

I 

Int 

— 

IREP 

Int 

— 

I  SPEC 

Int 

— 

I  START 

Int 

2 

ISTOP 

Int 

2 

LOOP 

Int 

— 

NPLOTS 

Int 

— 

PCTCRZ 

Real 

— 

PCTSAW 

Real 

— 

PLTSIZ 

Real 

— 

POINT 

Int 

— 

PROD 

Real(C*l) 

— 

PULPT 

Real 

35 

PULSQ 

Real 

35 

PULVOL 

Real 

— 

REGWGT 

Real 

— 

SPEC 

Int(C*2) 

— 

TRECNT 

Int 

40,32 

TREES 

Int 

— 

VARCUB 

Real 

— 

Description 


Board  foot  volume  of  a  tree 

Sawlog  bd.ft.  volume  of  a  sample  unit 

Sum  of  BRDPT^ 

Sawlog  cub.ft.  volume  of  a  sample  unit 

Sum  of  CUBPt2 

Cubic  foot  volume  of  a  tree 

Diameter  at  breast  height 

Sample  method  expansion  factor 

Tree  height 

DO  loop  index 

Number  of  reps  on  input 

Species  loop  counter 

First  of  a  forest  type 

Last  of  a  forest  type 

Loop  counter 

Number  of  plots  sampled 

Percent  cruise 

Percent  sawlog 

Fixed  radius  plot  size 

Point  number 

Product  code 

Pulpwood  cub.ft,  volume  of  a  sample  unit 

Sum  of  PULPt2 

Pulpwood  cub.ft.  volume  of  a  tree 

Weight  for  regression  standard  error 

Species  input  code 

Number  of  trees  tallied  by  d.b.h. 

Number  of  trees  on  a  sample  unit 

Sawlog  cub.ft.  volume  variance 


Subroutine  VOLMER 

VOLMER  is  called  by  GUTS  to  calculate  the  merchantable  cubic  and  board 
foot  volume  of  individual  trees  when  tree  height  is  specified  to  a  merchant- 
able top.  The  input  parameters  are  d.b.h.,  merchantable  height,  and  species. 

The  cubic  foot  volume  equations  are  from  Barnard  et  al.  (1969).   To 
use  these  equations  all  species  must  be  placed  in  one  of  seventeen  groups. 
The  species  grouping  is  as  follows: 
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Spec 

ies 

in 

Group 

WP, 

RP, 

OH 

SP 

BF 

HM 

PP, 

TA, 

NS,  SC 

CE 

SM 

RM, 

YP 

WA, 

AS 

BC 

WB, 

YB, 

SB,  GB 

BE 

B 

RO, 

BO, 

GU 

WO, 

EL 

HI 

WI, 

HH, 

OH 

Species  Group 
Number 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 


_   The  board  foot  volume  equation  is  from  Wiant  and  Castenada  (1977), 

This  equation  is  a  formulation  of  the  Mesavage  and  Girard  (1946)  form  class 

tables.  A  form  class  must  be  specified  for  each  of  the  32  species.  VOLMER 

sets  the  form  classes  to  the  following  defaults : 

Hemlock  70 

All  other  softwoods  76 
All  hardwoods       74 

The  user  can  specify  alternate  form  classes  at  execution  time. 

The  cubic  foot  equations  are  based  on  a  4"  merchantable  top  while  the 
board  foot  equations  are  based  on  a  10"  merchantable  top.  The  equations 
are  presented  in  full  in  Appendix  C,  pages  56  and  57. 

Common  Area  Used:   FORMCL. 


Other 


Variabl 

es  Type 

D 

imension 

BDHT 

Real 



BDVOL 

Real 

— 

CBHT 

Real 

— 

CUBINT 

Real 

17 

CUBSLP 

Real 

17 

CUBVOL 

Real 

— 

DBH 

Int 

— 

GRP 

Int 

32 

HEIGHT 

Real 

— 

I 

Int 

— 

ISPEC 

Int 

— 

WIANT 

Real  : 

Functii 

on 

Description 

Board  ft.  merchantable  height  in  logs 

Board  ft.  volume  of  tree 

Cubic  ft.  merchantable  height  in  feet 

Cubic  ft.  volume  equation  intercept 

Cubic  ft.  volume  equation  slope 

Cubic  ft.  volume  of  tree 

Diameter  at  breast  height 

Species  group  for  cub.ft.  volume  equation 

Tree  height 

Implied  DO  loop  index  for  data  statements 

Species  loop  counter 

Computes  board  foot  volume 
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Sub routine  VOLTOT 

VOLTOT  is  called  by  GUTS  to  calculate  the  merchantable  cubic  and  board 
foot  volume  of  individual  trees  when  total  tree  height  is  measured.  The 
input  parameters  are  d.b.h.,  total  height,  and  species. 

Both  the  cubic  and  board  foot  equations  are  from  Honer  (1967).   In 

applying  these  formulas,  the  gross  cubic  foot  volume  of  the  tree  is  first 

calculated.  The  merchantable  cubic  and  board  foot  volumes  are  then  derived 

from  the  gross  volume.  To  use  these  equations  it  is  necessary  to  place  all 

species  into  one  of  16  groups.   The  species  grouping  is  as  follows: 

Species  Group  Species  in 

Number  Group 

1  WP 

2  RP,  PP.  TA,  SC,  OS 

3  SP,  NS 

4  BF 

5  CE 

6  HM 

7  AS 

8  WB,  SB,  GB 

9  YB 

10  SM,  RM,  YP 

11  BW 

12  BE 

13  BC 

14  EL 

15  WI,  HH,  OH 

16  RO,  WO,  BO  GU,  HI,  WA 

Stump  height  and  merchantable  top  are  variables  in  the  equations. 

VOLTOT  sets  these  to  the  following  defaults: 

Stump  Height  0.5  ft. 

Cubic  Foot  Top    4  in. 

Softwood  Board  Foot  Top    6  in. 

Hardwood  Board  Foot  Top    8  in. 

The  equations  are  presented  in  full  in  Appendix  C,  page  58. 


No  common  area  used. 
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Variables 

.  Type 

Dimension 

ASDRBF 

Real 

3 

ASDRCF 

Real 

3 

BDVOL 

Real 

— 

BRDTOP 

Real 

2 

CUB INT 

Real 

16 

CUBSLP 

Real 

16 

CUBTOP 

Real 

2 

CUBVOL 

Real 

— 

DBH 

Int 

— 

GRP 

Int 

32 

HEIGHT 

Real 

— - 

I 

Int 

— 

ISPEC 

Int 

— 

ITYPE 

Int 

— 

STUMP 

Real 

— 

TCBVOL 

Real 

— 

XBRD 

Real 

_ 

XCUB 


Real 


Description 

Board  foot  volume  equation  coefficients 

Cubic  foot  volume  equation  coefficients 

Sawlog  bd.ft.  volume  of  tree 

Bd.ft.  merchantable  top  diameter 

Cub.ft.  volume  equation  intercept 

Cub.ft.  volume  equation  slope 

Cub.ft.  merchantable  top  diameter 

Cubic  ft.  volume  of  tree 

Diameter  at  breast  height 

Species  group  for  volume  equations 

Tree  Height 

Implied  DO  loop  index  for  data  statements 

Species  loop  counter 

Forest  type  indicator 

Stump  height 

Total  cubic  foot  volume 

Intermediate  step  in  bd.ft.  volume 

calculation 
Intermediate  step  in  cub.ft.  volume 

calculation 


Subroutine  STRAT 

STRAT  is  called  by  INVENT  if  more  than  one  stratum  are  being  processed. 
STRAT  calculates  combined  strata  statistics  for  sawlog  cubic  foot  volumes, 
sawlog  board  foot  volumes  and  pulpwood  cubic  foot  volumes.  These  statistics 
include  total  volume  over  all  strata  along  with  their  standard  error  and 
effective  degrees  of  freedom.   STRAT  also  calculates  the  total  acreage  of 
all  the  strata  combined. 

STRAT  has  two  entry  points,  STRAT  and  STRAT2.   STRAT  is  used  as  an 
entry  point  after  each  individual  stratum  is  processed.   Sums  are  accumulated 
for  volumes,  variances,  and  squared  variances.   STRAT2  is  called  after  all 
of  the  strata  have  been  processed.  The  previously  accumulated  sums  are  used 
to  calculate  total  volumes,  standard  errors,  and  effective  degrees  of  freedom. 


Common  Areas  Used:   OUTl  and  STOUT. 


Other 

Variables 

Type 

Dimension 

ACRES 

Real 

— 

CLASIZ 

Int 

— 

DBH 

Int 

— 

DFBD2 

Real 

32 

DFCB2 

Real 

32 

DFPL2 

Real 

32 

ISPEC 

Int 

— 

STVMBD 

Real 

35 

STVMCB 

Real 

35 

STVMPL 

Real 

35 

I 

Int 

— 

J 

Int 

_ 

Description 

Acreage  of  compartment 
Diameter  class  size 
Diameter  at  breast  height 
Intermediate  calculation  in  EFDFBD 
Intermediate  calculation  in  EFDFCB 
Intermediate  calculation  in  EFDFPL 
Species  loop  counter 
Sawlog  bd.ft.  variance 
Sawlog  cub.ft.  variance 
Pulpwood  cub . ft .  variance 
Implied  DO  loop  index 
Implied  DO  loop  index 
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Subroutine  OUTPl 

OUTPl  is  called  by  INVENT  to  produce  all  of  the  program's  output.   The 
subroutine  has  two  entry  points,  OUTPl  and  0UTP2.   OUTPl  is  used  to  produce 
the  individual  stratum  output  while  0UTP2  is  used  to  produce  the  stratified 
total  output.   All  output  is  written  to  unit  3  (LPT  on  DEC-10)  creating  a 
file  named  INVENT. LPT. 

OUTPl  produces  the  following  output  tables  for  each  stratum: 

1.  Summary  of  job  control  related  information  (Table  1). 

2.  Analysis  of  the  height-dbh  relationship  if  regression  option 
used  (Table  17) . 

3.  Individual  stand  and  stock  table  for  each  species  sampled 
(Tables  2,  3). 

4.  Stand  and  stock  tables  for  Softwood  summary.  Hardwood  summary, 
and  All  Species  summary  (Tables  4,  5,  6). 

5.  Table  of  species  composition  by  percent  (Table  7). 

6.  Table  of  volume  totals  for  all  species  (Table  8). 

7.  Table  of  volume  totals  for  all  species  expanded  by  acreage 
(Table  9). 

The  stand  and  stock  tables  include  basal  area/acre,  trees/acre,  sawlog  cubic 

and  board  foot  volumes  per  acre,  and  pulpwood  cubic  foot  volume  per  acre, 

all  presented  by  diameter  class  as  well  as  a  total  over  all  diameters.  A 

confidence  interval  expressed  as  a  percent  is  given  for  each  volume  total. 

Also  listed  are  mean  stand  diameter,  merchantable  m.s.d.  (mean  stand  diameter 

of  trees  with  dbh  greater  than  or  equal  to  8  inches),  percent  cruise,  sawlog 

board  foot  coefficient  of  variation,  and  pulpwood  cubic  foot  coefficient  of 

variation.   Furthermore,  the  title  of  each  compartment  is  printed  on  the  top 

of  each  page  and  the  pages  are  consecutively  numbered. 

0UTP2  is  used  as  an  entry  point  to  produce  the  following  stratified  total 

output  if  more  than  one  stratum  were  sampled: 

1.  Summary  of  job  control  related  information  (Table  10). 

2.  Volume  tables  for  each  species  sampled  (Table  11,  12). 

3.  Volume  tables  for  Softwood  summary.  Hardwood  summary,  and  All  Species 
summary  (Tables  13-15). 

4.  Table  of  volume  totals  for  all  species  (Table  16). 
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All  volume  figures  from  0UTP2  are  total  volumes  expanded  over  the  entire 
acreage  of  all  the  strata  combined.   The  individual  species  tables  and  the 
three  summary  tables  list  volumes  by  diameter  class  as  well  as  giving  a 
total  over  all  diameters.   A  confidence  interval  expressed  as  a  percent  is 
given  for  each  volume  total.   Each  page  of  the  output  is  consecutively 
numbered  and  a  specified  title  is  printed  at  the  top  of  each  page. 

Both  OUTPl  and  0UTP2,  when  producing  individual  species  tables,  will 
print  the  form  class  used  to  calculate  board  foot  volume  next  to  the  species 
title  if  the  trees  were  measured  to  a  merchantable  top.   Otherwise,  no  form 
class  information  is  printed. 

Common  Areas  Used:   COMUNE,  FORMCL,  OUTl,  REG,  and  STOUT. 


Other 


Variable 

es  Type 

Dimension 

AFCNUM 

ReaKC 

*2) 



AFCTIT 

ReaKC 

*5) 

— 

ALPHA 

Real 

— 

BASCOM 

Real 

— 

BDCOM 

Real 

— 

CUBCOM 

Real 

— 

DBH 

Int 

— 

DP 

Real 

— 

I 

Int 

— 

IBDCV 

Int 

— 

IBDTST 

Int 

— 

IBDVST 

Int 

— 

ICBTST 

Int 



ICBVST 

Int 

— 

IDATE 

int(c*: 

5) 

2 

lERBRD 

Int 

— 

lERCUB 

Int 

— 

lERPUL 

Int 

— 

IFLAG 

Int 

— 

IPAGE 

Int 

— 

IPLCV 

Int 

— 

IPLTST 

Int 

— 

IPLVST 

Int 

— 

ISPEC 

Int 



ISTRA 

Int 

— 

J 

Int 

— 

MSD 

Real 

— 

OFRAME 

ReaKC^ 

*5) 

— 

PCTCRZ 

Real 

— 

PULCOM 

Real 

— 

RBSTOT 

Real 

— 

RECUSE 

ReaKC^ 

*A) 

— 

RMSD 

Real 

— 

RTRTOT 

Real 

— 

SPETIT 

Int(C* 

5) 

35,3 

Description 

Form  class  number  as  a  literal 

Output  title 

Confidence  level  as  a  decimal 

Percentage  of  basal  area 

Percentage  of  sawlog  bd.ft.  volume 

Percentage  of  sawlog  cub.ft.  volume 

Diameter  at  breast  height 

Degrees  of  freedom 

Implied  DO  loop  index  for  data  statements 

Board  ft.  coefficient  of  variation 

Total  sawlog  bd.ft.  volume  all  strata 

Total  sawlog  bd.ft.  volume  all  strata  by 

dbh 
Total  sawlog  cb.ft.  volume  all  strata 
Total  sawlog  bd.ft.  volume  all  strata  by 

dbh 
Date 

Percent  error  in  sawlog  bd.ft.  volume 
Percent  error  in  sawlog  cb.ft.  volume 
Percent  error  in  pulpwood  cb.ft.  volume 
Warning  message  flag 
Page  number 

Pulpwood  coefficient  of  variation 
Total  pulpwood  cb.ft.  volume  all  strata 
Total  pulpwood  cb.ft.  volume  all  strata 

by  dbh 
Species  loop  counter 
Number  of  strata  sampled 
Implied  DO  loop  index 
Mean  stand  diameter 
Saves  the  contents  of  FRAME 
Estimated  percent  cruise  for  point  sampling 
Percentage  of  pulpwood  volume 
Basal  area  of  trees  with  dbh  <^  8" 
Output  message 

Merchantable  mean  stand  diameter 
Trees/acre  of  trees  with  dbh  <^  8" 
Species  titles  for  output 
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Other 


Variabl 

.es 

Type 

Dimension 

TOTABD 

Int 

__ 

TOTACB 

Int 

— 

TOTAPL 

Int 

— 

TOTPTS 

Int 

— 

TRECOM 

Real 

— 

TVAL 

Real 

— 

Description 

Total  sawlog  bd.ft.  volume 
Total  sawlog  cb.ft.  volume 
Total  pulpwood  volume 
Total  number  of  sample  units 
Percentage  of  trees/acre 
t -value 


Subroutines  MDSTI  and  MDNRIS 

MDSTI  and  MDNRIS  are  adapted  from  the  International  Mathematical  and 
Statistical  Libraries,  Inc.  (IMSL) .   Both  subroutines  have  been  edited 
removing  capabilities  not  needed  in  this  application.  The  remaining  code 
is  reproduced  in  Appendix  E  in  accordance  with  the  IMSL  policy  on  research 
work  as  stated  on  page  INTRO-22  of  the  July,  1977  IMSL  library  manual. 

MDSTI  is  called  by  OUTPl  to  calculate  the  percentage  point  of  the 
Student's  t  distribution.  MDNRIS  is  called  by  MDSTI  and  is  used  to  evaluate 
the  inverse  normal  probability  distribution  function.  MDSTI  gives  5  signif- 
icant digit  accuracy  (IMSL  1977).   Furthermore,  table  26.10,  p.  999  of  the 
Handbook  of  Mathematical  Functions  edited  by  Abramowitz  and  Stegun  was 
duplicated  by  MDSTI  (IMSL  1977) . 

No  variable  dictionary  is  presented.  The  interested  reader  may  refer 
to  pages  MDSTI-1,  MERFI/MERFCI/MDNRIS-1,  and  MERFI/MERFCI/MDNRIS-2  of  the 
1977  IMSL  Library  manual. 
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SECTION  4.   STATISTICAL  ANALYSIS 

Overview 

The  purpose  of  INVENT 's  statistical  analysis  is  to  indicate  the  reli- 
ability of  the  volume  estimates.   This  is  accomplished  by  computing  a 
confidence  interval  for  all  volume  estimates  that  are  summed  across  the 
diameter  classes.   This  confidence  interval  is  specified  as  a  percent  of 
the  respective  volume  estimate.  The  user  specifies  the  confidence  level 
as  a  decimal  in  the  range  .65  to  .99.   All  statistical  analysis  assumes  a 
stratified  random  sample  where  each  sample  unit  is  a  cluster.   There  is  an 
option  to  double  sample  for  height  within  the  clusters. 

INVENT  calculates  confidence  intervals  for  sawlog  cubic  foot  volume, 
sawlog  board  foot  volume,  and  pulpwood  cubic  foot  volume.   The  analysis 
for  each  volume  is  carried  out  independently  of  the  other  volumes.   Since 
the  method  is  exactly  the  same  for  each  volume,  the  calculation  of  error 
for  an  arbitrary  volume  V  is  presented.   Throughout  the  discussion  we  will 
assume  the  reader  is  familiar  with  the  techniques  for  calculating  the  mean 
volume  per  acre  when  the  heights  of  all  volume  trees  are  measured.   (Cf. 
Barrett  and  Nutt,  1975;  Avery,  1975;  Husch,  Miller,  and  Beers,  1972; 
Dillworth  and  Bell,  1976.) 

The  statistical  analysis  is  presented  in  the  same  sequence  as  the 
program  performs  it.   For  each  individual  stratum  the  standard  error  of  the 
volume  must  be  calculated.   If  the  height  subsampling  option  is  used,  the 
program  first  calculates  the  standard  error  due  to  height  estimation.  Next 
the  percent  cruise  of  the  stratum  is  calculated.  This  is  used  to  refine 
the  standard  error  due  to  the  variation  between  clusters  which  is  calculated 
next.   This  standard  error  is  then  combined  with  the  height  estimation 
standard  error  to  give  the  volume  standard  error  of  the  stratum.   Finally, 
the  volume  standard  errors  of  all  the  strata  are  combined  to  yield  the 
standard  error  of  the  stratified  volume  estimate  along  with  the  effective 
degrees  of  freedom  of  that  estimate. 
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Standard  Error  of  Height  Estimation 
INVENT  allows  the  user  to  double  sample  within  each  cluster  for  tree 
height.   For  each  species  where  three  or  more  height  trees  are  measured 
within  the  stratum,  a  regression  relationship  between  height  and  the  inverse 
of  d.b.h.  is  analysed.   When  calculating  volume  for  a  tree  whose  height  was 
not  measured,  the  regression  equation  is  used  to  generate  a  height  for  that 
tree.   If  the  user  were  required  to  measure  the  heights  of  all  trees  on 
some  of  the  clusters  while  measuring  no  heights  on  the  remaining  clusters, 
the  volume  standard  error  could  be  calculated  using  techniques  developed  by 
Donald  Bruce  (1961)  and  Floyd  A.  Johnson  (1958).   Bruce 's  method  is  outlined 
in  Appendix  B,  page  53.   Johnson's  method  is  an  application  of  double 
sampling  with  ratio  of  means  estimation  and  this  method  is  presented  in 
Appendix  B,  page  52.  Neither  method  Is  applicable  if  the  user  is  free  to 
subsample  a  portion  of  the  trees  within  each  cluster.  A  new  statistical 
technique  was  developed  to  approximate  the  standard  error  of  the  resulting 
volume  estimate. 

The  authors  believe  that  this  new  technique  is  far  more  efficient  in 
field  application.   It  allows  the  user  to  sample  heights  more  intensely  on 
the  more  valuable  species.   It  provides  far  more  flexibility  than  the 
traditional  approach. 

H  =  height    D  =  d.b.h.    T  =  //  trees  sampled    t  =  //  height  trees  sampled 
Regression  Model:   H  =  B^  +  B^  ^"n"^  "*"  ^ 

„   ^      /-u^     c2     MSE 
Variance  (H)  =  S„  =  — — - 
u     t 


where  MSE  =  mean  square  error  of  the  regression  (for 

computational  formula  see  Neter  and  Wasserman, 
1974,  p.  45). 

^H 
Standard  Error  (H)   =  SE„  =  — 2_  (iqO) 

"      H 

SE  is  calculated  for  every  species  where  three  or  more  height  trees 
H 

are  measured.   It  is  also  necessary  to  calculate  SE^  for  the  Softwood,  Hard- 
wood and  All  Species  summary  volumes.  The  method  used  is  from  Neter  and 
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Wasserman  pp.  163-4  (1974).   It  is  equivalent  to  fitting  a  regression  for 

each  of  the  three  groups  using  pooled  data  from  each  species  in  the  group. 

However  the  method  used  is  far  more  efficient.   We  simply  define  T  ,  t  ,  and 

MSE  for  each  summary  group  as  a  function  of  T,  t,  and  MSE  for  every  species 

in  that  group.   For  instance,  suppose  we  want  SE  for  hardwoods  where  the 

n 

hardwood  group  in  this  stratum  consists  of  three  species.  Then, 

* 

T   =  T  +  T  +  T 
1    2   ^3 

.            MSE, (t,-2)  +MSE_(t,-2)  +  MSE-(t_-2) 
MSE   =  ^^         2  2         3  3 


(t^-2)  +  (t2-2)  +  (t3-2) 


Using  these  new  values,  SE  is  calculated  as  before. 

H 

Percent  Cruise 
Percent  cruise  is  simply  the  percent  of  the  population  actually  sampled. 
When  sampling  with  fixed-radius  plots,  percent  cruise  is  easily  calculated 
as  follows: 

Let       P  =  percent  cruise 

PLT  =  plot  size  in  acres 

n  =  number  of  plots  sampled 
AC  =  tract  size  in  acres 
Then.      P  =  "  '^J^'^        (100) 

When  sampling  with  variable-radius  plots  (point  sampling)  the  true  value 
of  P  is  unknown.   However,  based  on  our  sample  estimate  of  MSD  we  can  estimate 
the  value  of  P  (Cunia,  1959). 

Let      P  =  percent  cruise 

n  =  number  of  sample  points 
AC  =  tract  size  in  acres 
BAF  -  basal  area  factor  of  angle  used 

B  =  total  basal  area  of  all  trees  on  tract 
T  =  total  number  of  trees  on  tract 

MSD  =  mean  stand  diameter  (quadratic  mean) 

100  •  n  •  B    _  0.5454  •  n  •  MSD^ 
'  T  •  AC  •  BAF  ~        AC  •  BAF 
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It  should  be  noted  that  in  a  variable-radius  cruise,  the  percent  cruise  may 
be  different  for  each  species.   See  tables  2  and  3. 

In  forest  sampling,  we  deal  with  finite  populations.  The  formula  we 
use  to  estimate  the  variance  of  the  mean  volume  per  acre  for  each  stratum 
includes  a  term  known  as  the  finite  population  correction  (FPC) .  The  FPC 
reduces  the  variance  of  the  mean  when  a  significant  (usually  >  5  percent) 
portion  of  the  population  has  been  sampled. 

Let  N  =  number  of  sample  units  in  population 
n  =  number  of  units  actually  sampled 

N  -  n 
Then,  FPC  =       "       "^ 


N 

AC 
In  plot  sampling  N  =         and  n  =  number  of  plots  sampled.   However, 

if  Lt  L 

in  point  sampling  there  is  no  easy  was  to  calculate  N.  Therefore,  mensura- 
tionists  have  traditionally  ignored  the  FPC  when  analysing  point  cruises. 
This  can  be  quite  inefficient  when  tract  sizes  are  small  as  is  often  the 
case  in  New  Hampshire.   This  problem  can  be  easily  overcome  using  the 
previous  result  permitting  the  calculation  of  percent  cruise  when  point 
sampling.   It  can  easily  be  shown  that  the  FPC  is  simply  a  function  of  P 
(percent  cruise) . 

FPC  =   N  -  n   ^  _N D_  =  1  _  ^  =  1     P 


N        N      N  N         100 

INVENT  uses  this  relationship  to  refine  the  variance  estimates  when  point 

p 
sampling  as  well  as  plot  sampling.   INVENT  calculates  FPC  as  1  -  y^  . 

The  standard  error  due  to  the  variation  in  volume  between  clusters,  SE^, 

is  calculated  as  follows: 

Y  =  mean  volume /acre  of  stratum 

2 
S  =  sample  variance  of  volume 

n  =  number  of  sample  units 

1 


then,   SE^  = 


(100)   (FPC) 


If  the  height  of  all  volume  trees  was  measured,  then  SE—  is  the  standard 
error  of  the  volume  estimate  for  the  stratum,   SE   .   If,  however,  the 
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height  double  sampling  option  was  used,  SEL   must  reflect  both  SE—  and  SE„, 
the  standard  error  due  to  height  estimation.   In  this  case,  SE^   is 
calculated  as  follows : 


^^voL  =  ^/  4  ^  -H^  ^'4 ) 


Combined  Strata  Statistics 
We  have  seen  how  INVENT  calculates  the  mean  volume/acre  with  its  standard 
error  for  each  stratum  in  the  sample.  The  final  step  in  the  statistical 
analysis  is  to  use  the  individual  stratum  statistics  to  calculate  the  total 
volume  across  all  strata  along  with  the  standard  error. 
Let       L  =  number  of  strata  sampled 

AC.  =  size  of  stratum  i  in  acres 

N  =  number  of  sample  units  in  stratum  i 

/ 

n  =  number  of  units  sampled  in  stratum  i 

Y  =  mean  volume  per  acre  of  stratum  i 


VOL  =  total  volume  of  stratum  i  =  AC  "  Y 


^% 


OL  =  volume  standard  error  of  stratum  i  as  a  percent 

2 

2        2       2 
VOL  =   (SE^Q,  •  VOL  )  /  100   =  variance  of  total  volume 


1 
then  we  calculate  the  stratified  total  statistics  as  follows: 


AC  =  Z     AC   =  total  acreage  sampled 
i=l   ^ 


L 

VOL  =   Z     VOL   =  stratified  total  volume 
i=l    ^ 


J I  'K 


»^0L "°« 

Having  calculated  VOL  and  SE^  ,  all  that  remains  in  calculating  the 
confidence  interval  is  determining  the  appropriate  percentage  point  of  the 

Student's  t  distribution.   In  order  to  do  this  we  must  calculate  the 

2 
degrees  of  freedom  associated  with  S    .   Satterthwaite  (1946)  states  that 

2 
the  exact  distribution  of  the  complex  estimate  S    is  too  involved  for 

VULi 
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everyday  use.   He  suggests  approximating  the  true  distribution  with  a  Chi- 
square  distribution  of  equal  variance.   Cochran  (1977  p.  96)  uses  this 
approach  and  defines  the  degrees  of  freedom  of  the  approximating  Chi-square 
as^ 


DF 


eff 


i=l   ^ 

i=i  ^i"^ 


where  g . 


"i  ^^i  "  ^i^ 


S.   =  sample  variance 


A  more  appropriate  form  for  our  purposes  is  given  by 
L 


DF 


eff 


ifl  '^°^i 


L   VOL . 
i=l  ^i-  ^ 


The  equivalence  of  the  two  formulas  is  shown  in  Appendix  B,  page  54. 

Using  this  result  the  appropriate  value  t  is  derived  from  the  Student's 
t  distribution  and  the  confidence  interval  is  given  by: 


VOL  ±  fSE. 


VOL 


-43- 


APPENDIX  A 


DATA  CODING  INFORMATION 
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EXAMPLES  OF  TREE  CODING  OPTIONS 
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INVENT  Species  List 


White  Pine 

WP 

Red  Pine 

RP 

Pitch  Pine 

PP 

Hemlock 

HM 

Balsam  Fir 

BF 

Tamarack 

TA 

Spruce 

SP 

Norway  Spruce 

NS 

Scotch  Pine 

SC 

Cedar 

CE 

Other  Softwoods 

OS 

Sugar  Maple 

SM 

Red  Maple 

RM 

Ash 

WA 

Willow 

WI 

Aspen 

AS 

Yellow  Poplar 

YP 

Black  Cherry 

BC 

White  Birch 

WB 

Yellow  Birch 

YB 

Sweet  Birch 

SB 

Grey  Birch 

GB 

Beech 

BE 

Basswood 

BW 

Red  Oak 

RO 

White  Oak 

WO 

Black  Oak 

BO 

Elm 

EL 

Gum 

GU 

Hickory 

HI 

Hophornbeam 

HH 

Other  Hardwoods 

OH 
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INVENT  TALLY  SHEET 


COMR 

^ 

POINT 

± 

ff 

TREES 

1  1 

SPECIES 

DBH 

HEIGHT 

PRODUCT 

i=t: 

f^P 

tW 

l.O 

wP 

ao 

H.O 

iRP 

/^ 

z.^ 

WjP 

/9r 

H.O 

wP 

o^ 

s.o 

wP 

I"? 

l.iT 

wp 

l(o 

z.o 

wp 

0.0 

c/.O 

wp 

XX 

v.r 

U^\P 

-x^ 

V'^T 

RP 

l(o 

3.0 

COMR 

/ 

L 

POINT 

^ 

^           /^ 

TREES 

^ 

SPECIES 

DBH 

HEIGHT 

FRODUC  T 

i=tr 

luP 

207 

-i.O 

^uP 

^2. 

^.  s 

RP 

2^ 

3.r 

RP 

l^ 

2>.S" 

({o 

N 

LO 

(lo 

:i^ 

\.o 

Ro 

1? 

1.  o 

wp 

Q.L, 

?'C 

^ 

COMR 
POINT 

a 

^       /^ 

TREES 

/^ 

SPECIES 

DBH 

HEIGHT 

PRODUCT 

:«r 

vue 

/z 

a  r 

LV-'P 

9.0 

3.6~ 

^'P 

3a 

v.o 

(ip 

/6 

3.0 

a^ 

/  o 

/.o 

Ce 

OG 

•  r 

P 

A 

ujp 

/(:? 

3>-S 

UjP 

/G 

^.o 

wP 

O 

ao 

Rp 

/<^ 

a.o 

wp 

^(o 

'V.O 

COMR 
POINT 

i 

V           ^ 

■TREES 

IV 

SPECIES 

DBH 

HEIGHT 

PRODUCT 

:fct 

u^P 

2-:i 

"i.O 

wP 

2  o 

H-O 

u^f 

3  0 

"i-o 

( 

wP 

:^o 

Z.'T 

ucP 

ic 

H^O 

3. 

U-P 

a.(^ 

9.S- 

Cc/P 

/t- 

"i.O 

U.P 

^? 

H.  r 

u/P 

/§■ 

l.o 

U/  P 

:ikc 

w.  s~ 

wP 

?.? 

^.  s~ 

wP 

5:i, 

w.o 

^p 

3LC 

^.s~ 
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INVENT  DATA  FILE 


WHITAKER  WOODS 

01  11 

RP  14 

3.0 

WP  20 

4,0 

RP  14 

3.5 

WP  18 

4.0 

WP  12 

3.0 

WP  18 

3.5 

WP  16 

3.0 

WP  20 

4.0 

WP  22 

4.5 

WP  22 

4.5 

RP  16 

3.0 

02  12 

WP  12 

2.5 

WP  20 

3.5 

WP  32 

4.0 

RP  16 

3.0 

BE  10 

1.0 

BE  06P0.5  2 

WP  16 

3.5 

WP  16 

4.0 

WP  12 

2.0 

RP  16 

3.0 

WP  26 

4.0 

03  09 

WP  26 

4.0 

WP  32 

4.5 

RP  22 

3.5 

RP  18 

3.5 

RO  14 

1.0 

RO  22 

1.0 

RO  18 

1.0 

WP  26 

5.0  2 

04  14 

WP  22 

4,0 

WP  30 

4.0 

WP  30 

4.0 

WP  20 

3.5 

WP  20 

4.0  2 

WP  26 

4.5 

WP  16 

1 . 0 

UP  18 

4.5 

WP  18 

3.0 

WP  20 

4.5 

WP  28 

4.5 

WP  32 

4.0 

WP  20 

4.5 

05  13 

WP  24 

3.0 

RP  16 

2.0 

WP  32 

4.0 

WP  26 

3.5 

RO  12 

1.0 

HM  18 

1.5 

HM  .1.6 

1.0 

RP  12 

1.0 

RP  10 

1.0 

RO  16 

1.0  2 

UP  14 

2.5 

HM  18 

2.0 

1978         COMP.     1 
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INVENT    Processing  Request  Form 


Forester  Information: 


Name: 


Firm  Name;_ 
Address: 


Tract  Information: 


Number  of  compartments: 


Do  you  wish  to  use  alternate  form  classes?  (  yes  or  no  ) 

If  yes  enter  fona  classes  below. 


WP 

RP 

PP 

HM 

BF 

TA 

SP 

NS 

sc 

CE 

OS 

SM 

RM 

WA 

WI 

AS 

YP 

BC 

WB 

YB 

SB 

GB 

BE 

BW 

RO 

WO 

BO 

EL 

GU 

HI 

HH 

OH 

Diameter  class  size 


INVENT  uses  the  9^%   confidence  level  as  a  default.  If  you  wish 
to  use  a  different  level,  please  specify  (65  -  99) 

Tract  Description  (up  to  70  characters) 
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IKVENT    Compartment  Information  Form 


Compartment  Name  and/or  Description  (up  to  TO  characters) 


Sampling  Method  (  Point,  Plot,  Strip,  or  10055) 

If  Point:   BAT  = 


If  Plot  :  Plot  size  =  acres 

If  Strip:  Total  acreage  of  all  strips  =  acres 

Tree  Height  Measurement  Units  (  l6  ft.  logs  or  feet  ) 

Tree  top  specification  (  merchantable,  total,  or  mixed  ) 

Were  the  heights  of  all  volume  trees  measured?  (  yes  or  no  ) 

Acreage  of  compartment :  acres 
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APPENDIX  B 


STATISTICAL  FORMULAS  AND  DERIVATIONS 
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DOUBLE  SAMPLING  WITH  RATIO  OF  MEANS  EXPANSION 


X.      =   Basal  area  per  acre  at  cluster  1 


Y      =   Volume  per  acre  at  cluster  i, 
N       =   Number  of  possible  sample  units 
n      =   Total  number  of  units  sampled 


n      =   Number  of  volume  units  sampled 
o 


X      =   Average  basal  area  of  all  units  sampled 


X      =   Average  basal  area  of  the  volume  units  sampled 


Y      =   Average  volume  of  the  volume  units  sampled 
o 


Y      =   Ratio  estimate  of  the  mean  volume  per  acre 


Y 

B  -  -T- 

O 


Y  =  B  •  X 

r 


t  -  C-^X^l  (-ii^^^  4(-^ 
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STATISTICAL  ANALYSIS  OF  A  VARIABLE  PLOT  CRUISE  USING  BRUCE 'S  METHOD 

(Bruce  1961) 

Notation 
y  -       volume  points 

a      -   all  points 
n      -   number  of  points 


TC 


tree  count 


VBAR 


volume  to  basal  area  ratio 


VOL 


volume  per  acre 


SE 


standard  error  expressed  as  a  percent 


CF 


correction  factor 


Example  of  Notation; 

SE    -  standard  error  of  the  tree  count  on  all  points  expressed  as  a 
a   percent 


Estimated  volume  =  VOL 


TC  /n 

g  g 

TC  /n 
Y  Y  J 


=   (VOL  •  TC   •  n  )/(TC  •  n  ) 
Y    a    Y     Y   a 


/   2         2 

Uncorrected  standard  error  =  /  SE,„.„  +  SE 

.  /   VBAR      TC 
V       Y       a 


Corrected  standard  error  =  CF- (uncorrected  standard  error) 

/OL 


^hc 


where  CF  = 


<AR  +< 
Y       Y 
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FORMULA  FOR  EFFECTIVE  DEGREES  OF  FREEDOM 


DF 


eff 


i=l  "i"^ 


where  g . 


^i  ^^i-\^ 


L  N.  (N.-n^)   ^^y 
LV     n,    /  ^1 


1=1 


n^-1 


L  S^  _  (N,-n,)  n2  Y 
1=1  "l      ^1      J 


1=1 


n^-1 


a  \-  »o ' 

a  ^vo.)^ 

4         4 

J.          1 

i-l  "l"^ 

1=1    ^r^ 
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APPENDIX  C 


VOLUME  FORMULAS 
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CUBIC  FOOT  VOLUME  EQUATIONS 

BASED  ON  MERCHANTABLE  HEIGHT 

(Barnard  et  al.  1973) 


2 

General  form  of  equation:   VOLUME  =  a  +  b  (DBH  -HEIGHT) 


where  height  is  measured  to  a  A"  top  in  ft. 


Species  Group 

a 

b 

1 

3.5142 

0.00236 

2 

2.1998 

.00257 

3 

1.4793 

.00272 

4 

2.4784 

.00242 

5 

-0.0496 

.00303 

6 

1.5817 

.00259 

7 

1.6823 

.00310 

8 

1.1763 

.00310 

9 

1.0067 

.00293 

10 

1.1809 

.00292 

11 

1.1339 

.00281 

12 

1.2851 

.00326 

13 

0.8976 

.00349 

14 

1.2027 

.00301 

15 

1.0009 

.00300 

16 

1.4438 

.00316 

17 

0.8591 

.00309 

For  the  species  composition  of  each  group  see  page  31. 
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BOARD  FOOT  VOLUME  EQUATION 

BASED  ON  MERCHANTABLE  HEIGHT 

(Wlant  and  Castenada  1978) 
(INT.  1/4-lnch  Rule) 

D   =   diameter  at  breast  height 

H   =   merchantable  height  to  a  10"  top  in  16  foot  logs 

FC  =   Girard  form  class 


VOLUME  =  [ (aQ+a^H+a2H^)  +  (bQ+b^H+b2H^)D  +  (Cq+c^H+C2H^)D^]   [ (FC-78) (.03)+!] 


^0 

ss 

-13.35212 

^ 

= 

9.58615 

^2 

= 

1.52968 

^0 

= 

1.79620 

^1 

= 

-2.59995 

^2 

= 

-0.27465 

•=0 

= 

0.04482 

'^l 

= 

0.45997 

<=2 

= 

-0.00961 

Equations: 
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CUBIC  AND  BOARD  FOOT  VOLUME  EQUATIONS 

BASED  ON  TOTAL  HEIGHT 

(Honer  1967) 
(INT.  1/4-inch  Rule) 

D  =  diameter  at  breast  height 

H  =  total  height  in  feet 

GCVL  =  gross  cubic  foot  volume 

MCVL  =  merchantable  cubic  foot  volume 

MBVL  =  merchantable  board  foot  volume 

STH  =  stump  height 

CTOP  =  cubic  foot  merchantable  top  diameter 

BTOP  =  board  foot  merchantable  top  diameter 

GCVL  =     ° 


H 


X  = 


,_CTOPJ_, 


H  +  STH 


H 


MCVL  = 


[b  +  b,  X  +  b.  X  J  •  GCVL 
o    i      i. 


Y  = 


P  BTOP   ,   p  H  +  STH  , 


MBVL 


[c  +  c-  Y  +  c„  Y^]  •  GCVL 
o    i      / 


where 


STH 

= 

0. 

,5  ft. 

CTOP 

= 

4, 

.0  in. 

for  all  species 

BTOP 

= 

6, 

,0  in. 

for  softwoods 

8, 

.0  in. 

for  hardwoods 

t'l 

= 

0, 

.9604 

= 

-0, 

.1660 

b2 

= 

-0, 

.7868 

<^0 
'^l 
^2 

= 

5, 

.4332 

= 

-1, 

.6281 

= 

-4, 

.4710 

Species 

Group 

^0 

*1 

1 

0.691 

363.676 

2 

0.710 

355.623 

3 

1.226 

315.832 

4 

2.139 

301.634 

5 

4.167 

244.906 

6 

1.112 

350.092 

7 

-0.312 

436.683 

8 

2.222 

300.373 

9 

1.449 

344.754 

10 

1.046 

383.972 

11 

0.948 

401.456 

12 

0.959 

334.829 

13 

0.033 

393.336 

14 

0.634 

440.496 

15 

1.877 

332.585 

16 

1.512 

336.509 

For  the  species  composition  of  each  group  see  page  32. 
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APPENDIX  D 


INVENTORY  DESIGN  AND  SAMPLE  SIZE  ESTIMATION 
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Inventory  Design 

INVENT  analyzes  data  based  on  two  sampling  designs: 

(1)  simple  random  sampling, 

(2)  stratified  random  sampling  with  simple  random  sampling  within 
each  stratum. 

Sampling  Is  assumed  without  replacement — that  Is,  sample  units  (plots 
or  points)  are  far  enough  apart  so  that  a  tree  is  sampled  on  only  one  sample 
unit. 

INVENT  allows  for  the  possibility  of  double  sampling  regression  estima- 
tion of  volume.  The  error  associated  with  regression  estimation  is  small 
relative  to  the  error  among  sample  units.  Thus,  it  is  ignored  in  the 
discussion  given  below  on  sample  size  approximation. 

Foresters  commonly  use  systematic  sampling  rather  than  simple  random 
sampling.   Generally,  estimates  from  systematic  sampling  are  somewhat  more 
precise  than  from  simple  random  sampling  provided  the  systematic  sample 
units  are  spread  well  over  the  tract  sampled.  Therefore,  confidence 
intervals  based  on  simple  random  equations  provide  somewhat  conservative 
estimates  for  well-representative  systematically  collected  data.   For  the 
same  reason,  sample  size  is  usually  overestimated. 

Sample  Size 
Sample  size  is  treated  briefly.   More  detailed  explanations  are  given 
in  Barrett  and  Nutt  (1975) ,  Cochran  (1977) ,  Freese  (1962) ,  and  Dllworth  and 
Bell  (1976). 


If  a  tree  is  sampled  on  more  than  one  sample  unit,  the  estimated  means  are 
still  unbiased  but  the  variance  is  underestimated. 


-61- 

Slmple  Random  Sample 

Sample  size  from  simple  random  sampling  is  estimated  by  the  following 

equation: ^ 

r  t  CV  ,^ 
n  =  I — ]  Equation  1 

where  t  is  a  t-value,  CV  is  the  coefficient  of  variation,  and  E  is  the 
allowable  error  in  percent. 

The  t-value  depends  on  the  confidence  level  desired.  The  approximate 
t-values  for  the  most  widely  used  confidence  levels  are: 


Confidence  level 

(percent)  t-value 

80  1.3 

90  1.7 

95  2.0 

99  2.7 


The  coefficient  of  variation  is  a  measure  of  the  variability  among 
volume  estimates  in  a  forest  based  on  a  particular  technique  of  sampling. 

Assume  that  from  past  experience  a  forester  estimates  that  the 
coefficient  of  variation  is  50  percent  when  he  samples  with  a  10  B.A.F. 
prism.  He  wants  to  estimate  the  volume  on  a  tract  with  an  allowable  error 
of  10  percent  at  the  95  percent  confidence  level.  That  is,  he  wants  the 

chances  to  be  95  in  100  that  his  estimate  after  completing  the  inventory 

3 
is  within  10  percent  of  the  actual  volume. 

The  sample  size  by  Equation  1  is  estimated  as  follows: 

,  =  [_E_C^L]'  =   [^iP^]'  =  100  points 


Thus,  the  forester  must  select  100  points. 


2 

The  estimate  of  sample  size  using  simple  random  and  the  subsequent  estimates 
using  stratified  random  sampling  can  be  adjusted  by  the  equation 


n  = 


a   1  +  P/100 

where  P  is  the  percent  cruise.  This  additional  refinement  in  estimating 
sample  size  is  generally  not  necessary  except  on  small  tracts. 

3 
To  be  within  10  percent  of  the  actual  volume  implies  that  errors  such  as  tree 

measurement  errors,  errors  in  estimating  form  class,  and  errors  in  volume 

equations  are  negligible. 
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Stratified  Random  Sampling 

Often  a  forester  divides  the  area  into  compartments  or  stands  before 
beginning  the  inventory.   In  other  words  he  stratifies  the  area. 

If  the  forester's  primary  interest  is  in  the  estimates  within  each 
stratum,  then  he  can  specify  the  allowable  error  for  each  stratum  and  apply 
Equation  1  to  estimate  sample  size. 

If  his  primary  interest  is  in  the  mean  volume  per  acre  or  total  volume 
over  all  strata,  then  he  can  estimate  sample  size  by  one  of  two  techniques. 
Each  technique  depends  on  how  sample  units  are  allocated  to  strata.   The 
two  techniques  will  be  illustrated  with  an  example. 

First  consider  proportional  allocation.   By  proportional  allocation, 

the  forester  determines  sample  size  and  allocates  sample  units  to  strata  in 

proportion  to  the  acreage  of  each  stratum.   The  equation  for  sample  size  is 

t   2   L 
n  =  [-^]    I     W.CV7         Equation  2 

1=1 

where  t  is  the  t-value,  E  is  the  allowable  error,  W.  is  the  proportion  of 
acreage  in  stratum  i,  CV.  is  the  coefficient  of  variation  in  stratum  i, 
and  L  is  the  total  number  of  strata. 

Assume  a  forester  has  the  following  information  on  an  80-acre  tract  he 
wishes  to  sample  with  a  10-B.A.F.  prism. 


Stratum 

Acreage 

\ 

CV^ 

W^CV2 

1. 

2. 
3. 

small  pines 
large  pines 
mixed  pine-hardwood 

36 
18 
26 

.450 
.225 
.325 

30 
50 
80 

405.00 

562,50 

2080.00 

Totals 

80 

1.000 

— 

3047.50 

Suppose  he  wants  to  estimate  the  mean  volume  per  acre  within  10  percent 

with  95  percent  confidence.   Then,  by  applying  equation  2,  he  estimates  the 

number  of  points  required. 

2   L  2 

n  =  [-^]  Z     W^CV^  =  [j|-]   (3047.50)  =  122  points 

He  would  allocate  points  to  strata  as  follows: 
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Stratum         W  n'W 


1 

.450 

55 

2 

.225 

27 

3 

.325 

40 

Totals        1.000  122 


By  the  second  technique  of  allocation,  called  Neyman  allocation,  the  forester 
takes  into  account  the  coefficient  of  variation  as  well  as  the  relative 
acreage  in  each  stratum  in  assigning  sample  units  to  strata.  Neyman  allocation 
minimized  the  number  of  sample  units  required  to  achieve  a  specified  precision. 
The  equation  for  sample  size  is 

2 


n  =  i-^]' 


T.      W  CV 
1=1 


-i2 

Equation  3 


To  apply  this  equation  to  the  data  given  above  we  require  the  following 
computations : 


Stratum 

w. 

1 

CV^ 

W^CV^ 

1 
2 
3 

.450 
.225 
.325 

30 
50 
80 

13.50 
11.25 
26.00 

Totals 

1.000 

— 

50.75 

Assuming  an  allowable  error  of  10  percent  with  95  percent  confidence,  the 

forester  estimates  the  number  of  points  required. 

-i2 


n  =  [^]' 


1     W  CV 
i=l  ^ 


=  [j|-]   (50.75)^  =  103  points 


The  points  are  allocated  to  strata  as  follows ; 


L 

Stratum 

W.CV^ 

n-W  CV  /  l   W.CV. 
i  1  i=l  ^ 

1 

13 

.50 

27 

2 

11 

.25 

23 

3 

26 

.00 

53 

Totals 

50 

75 

103 

Proportional  and  Neyman  allocation  are  two  of  the  most  common  allocations. 
Other  techniques  of  allocation  are  given  in  the  references  cited. 
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APPENDIX  E 


PROGRAM  LISTING 
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c 
c* 

c 

c 
c 
c 
c 
c 
c 
c* 


PROGBAH   INVERT 
INVENT        VEB.    3  1/1/79         DAVID   S.    LINDEN      I.  N.  E.  R.      O.N.  H. 


THIS    IS    THE   HAIN    PBOGRAH    OF    INVENT 

IT   PROVIDES    THO    HAT    OSER   COHflDNICATION    AND  SOBBOOTINE   CONTROL 

DOOBLE   PBECISION   IFILE, SAHHTD, REPEAT 
INTEGER    LOG ANS, TOP ANS.CLASIZ, REGANS, TITLE (14) 

COHHON/COHONE/CLASI Z, TITLE, HO LT,PBOBLV,BAF, TOP ANS.BEGANS, ACRE 
1,  FRAllE,SAnHTD,  LOGANS,  IFILE 
DATA    PHOBLV/.05/ 
OPEN   OOTPOT   FILE 

OPEN (0NIT=3 , ACCESS= • SEQOOT* , FILE=« INVENT. LPT' ) 

OSER    COHBONICATION    PORTION 
****** ********* •#««****«•****** 

ASK    OSER    IF    HE    WANTS   TO   OSE   ALTERNATE   F.C.  •$ 
10    »RITE(5,20) 

20    FOBHAT(*0DO    ¥00    WANT   TO   READ  IN   ALTERNATE   FORM   CLASSES    •, 
1'FROa    DISK?*) 
READ(5,270)     FBHANS 
IF(FRHANS.BQ. 'YES')       CALL  FRllCLS 
IF<FBHANS.NE. 'NO*. AND.FRBANS.NE.'IES')       GOTO    10 
HAVE  OSER    SPECIFY   DBH   CLASS    SIZE 
WRITE(5,30) 
30    FORHAT(«0 ENTER    DIAHETER   CLASS   SIZE') 
READ(5,*)    CLASIZ 
HAVE  OSER   SPECIFY    NOBBEB   OF   STRATA 
NRITE(5,i(0) 
40    FORHATCO  ENTER    NOHBER   OF  STRATA') 
READ(5,*)    ISTRAT 
ASK    OSER    FOR    NEW   VALDE   OF    ALPHA    IF    DESIRED 
50    HBITE(5,60) 

60    FOHMATCOINVENT    OSES   THE   95X  CONFIDENCE   LEVEL   AS    A    DEPAOLT.  '  / 
1'    DO  YOO    WISH   TO   OSE   A    DIFFERENT   LEVEL?') 
READ(5,270)    ALPHCH 

IF(ALPHCH.NE.  ' YES*. AND. ALPHCH.  NE. 'NO')       GOTO   50 
IF (ALPHCH. Eft. 'NO')       GOTO    90 
70    HRITE(5,80) 

eo  FORHATC  ENTER  CONFIDENCE  LEVEL  AS  A  DECIHAL  IN  THE  RANGE  ', 
1'.65  TO  .99') 
READ(5,*)  ALPHA 

IF (ALPHA.GT. 0.99. OB.ALPHA.lt. 0.65)   GOTO  70 
PROBLV=1.0- ALPHA 
90  CONTINOE 


BEGIN  PROGBAB  CONTROL  LOOP 

DO  310  111=1, ISTRAT 

HAVE  OSER   SPECIFY    INPOT    FILE 

WRITE(5,100)    III 
100    FORBATCO  ENTER    INPOT    FILE   NABE   FOB    STRATOB    •    *,I3) 

READ(5,110)    IFILE 
110    FORBAT(AIO) 
OPEN   INPOT    FILE 

0PEN(0NIT=21,ACCESS='SEQIN',FILE=IFILE) 
BEAD  INVENTORY    TITLE    FROB    FILE 

HEAD(21,120)    TITLE 
120    FORBAT(14A5) 
WRITE  TITLE   ON   TEBBINAL 

WRITE (5, 130)    TITLE 
130    FORBATCO', 14A5/) 

SVFRAH=FRAaE 
HAVE  OSER    SPECIFY    SAHPLING    HETHOD 

IF(III.GT.  1)       BEPEAT='    OR   SABE' 
140    WRITE(5,150)    REPEAT 
150    FOBBATCOSPECIFY   SABPLING   HETHOD'/'         TYPE   POINT  OR    PL0T',A8) 

READ(5,240)  FBABE 

IF(FRAHE. EQ.' SABE'. AND. III. GT.1)     160  ,  170 


BAINOOIO 
BAINa02a 
IIAIN0030 

♦♦♦♦•BAINOOHO 
HAIN0050 
HAIN0060 
BAIN0070 
HAIH0080 
HAIN0090 
HAINOIOO 
HAIH0110 

♦•***BAIN0120 
BAIN0130 
HAIN0149 
HAlei0150 
BAIN0160 
HAIH0170 
HAIH0180 
BAIN0190 
HAIN0230 
HAIN0210 
BAIN9220 
HAIN0230 
BAIH0240 
HAIN0250 
HAIN0260 
NAIN0270 
HAIN0280 
BAIN0230 
aAIN0300 
BAIN0310 
BAtN032Q 
HAIN0330 
HAIN03II0 
HAIN0350 
BAIN0360 
HAIN0370 
aAIN0380 
aAIN0390 
HAIN0400 
BAIN3419 
BAIN0420 
BAIN0430 
BAIN044Q 
HAIN0450 
BAIN04&0 
BAIN0470 
BAINa43Q 
aAIN0490 
BAIN0500 
BAINQ510 
HAIN0520 
BAIN0530 
HAIN0540 
BAIN0550 
aAIN0560 
HAIN0570 
aAIN0580 
aAIN05)Q 
BAIN0600 
BAIN0610 
HAIN0620 
NAIN0630 
NAINOenO 
aAIH0650 
HAIN0660 
BAIN0670 
BAIN0680 
BAIN063O 
HAIN0700 
BAIN0710 
HAIH0720 
HAIN0730 
aAIN0740 
BAIN0750 
BAIN0760 
BAINa770 
HAIN0780 


-67- 


169    FEAHE=SVFEAH 

GOTO    260 
173    IF (FRiHE.NE-« POINT* -AND. FBAflE. HE. 'PLOT')       GOTO    1U0 

IF{FBAHE. EQ. 'POINT')       S AHHTD= 'B. A. F. • 

IF(FEAaE-EQ. 'PLOT')       SABJITD=' PLOT   SIZE' 
HAVE   USER   SPECIFY    BAF  OB    PLOT    SIZE 

HBITE(5,180)     SAHHTD 
180    FOBHATCOENTEB    '  ,  A9) 

READ(5,*)    BAF 
ASK    USEE   IF    HEIGHTS    ABE   IN    LOGS 
190    HBITE  (5,200) 

200    FOaaAT('0ARE   HEIGHTS    IN    FEET    OB    SIXTEEN    FT.     LOGS?'/ 
1'  TYPE   FEET    OE   LOGS') 

BEAD(5,210)     LOSANS 
210    FOBHAT(A(t) 

IF(LOGANS.NE. 'FEET'.  AND.LOGANS.NE. 'LOGS')   GOTO  190 
ASK  OSES  IF  HEIGHTS  ABE  HEBCHANTABLE  OE  TOTAL 
220  WBITE(5,230) 

230    FOBHAT('0ARE   HEIGHTS   TOTAL   OE    MEfiCH ANTABLE?'/ 
1'  TYPE  TOTAL,    HEBCH,    OB    flIXED') 

EEAD(5,240) TOPANS 
2U0    F0BBAT(A5) 

IF (LOGANS. EQ. 'FEET')       H0LT=1 

IF(LOGAHS.EQ. 'LOGS')       JJ0LT=16 

IF (TOPANS- NE. ' TOT AL' . AND. TOPANS- NE. • HEECH' . AND. TOPANS. NE.  'MIXED') 
1      GOTO    220 
ASK    OSES    IF   BEGBESSION    BOOTIBE   IS    NEEDED 
250    HEITE(5,260) 
260    FOBBATCOHEBE    THE    HEIGHTS    OF    ALL   TOLOBE   TBEES    BEASOEED?') 

BEAD(5,270)     BEGANS 
270   F0EHAT(A3) 

IF (BEGANS. NE. 'YES' .AND. EEGANS- NE. • NO')       GOTO    250 
ASK    USEE   FOB    ACBEAGE 
280    HBITE(5,290) 
290    FOBHATCOENTEB    ACBEAGE  OF  TEACT. ' ) 

BEAD(5,*)    ACEES 

IF  NEEDED  BON  THE  HEIGHT  BEGBESSION  BOOTINE 

***  **4>4<  ****«**•«««****  ******««4i**4i  41 4<***  ******  *4>*4<4i****4>*4i*4i« 


IF(BEGANS. E0.'HO')   CALL  LIHEEG  (IFILE, TOPANS) 

****************************************** 

BEAD  IN  TBEE  DATA  AND  PBOCESS 
****************************************** 

300  CALL  GUTS 

********************************** 
CALL  THE  OOTPOT  EOOTINE 
********************************** 

CALL  00TP1 

**************************************** 

ACCOHOLATE   TOTALS    FBOB    EACH    STBATA 
**************************************** 

IFIISTBAT.GT. 1)       CALL    STBAT (ACEES, CLASIZ) 

310    CONTINUE 

END    PEOGRAB   CONTROL   LOOP 


************************************* 
CALCULATE   COHBINED    STEATA    STATISTICS 
************************************* 

IF(ISTEAT. GT. 1)       CALL    STBAT2 


******************************** 
OOTPOT    STBATA    COBBINED   TOTALS 
******************************** 

IF(ISTEAr.GT. 1)       CALL   O0TP2 

EHD 


BAIN0790 
HAIN0800 
BAINOBIO 
HAIN0820 
flAIN0830 
BAINOSltO 
HAIN0850 
BAINOSSO 
aAIN0870 
HAIN0880 
flAIN0890 
HAIN0900 
BAIN0910 
BAIN0920 
BAIN0930 
BAIN09ltO 
HAIN0950 
BAIN09Sa 
HAIN0970 
HAIN0980 
BAIN099a 
HAIN1000 
BAIHIOIO 
HAIN1020 
HAIN1030 
BAIN10II0 
BAINIOSO 
HAINIOSO 
HAIN1070 
BAIN1080 
HAIN1090 
HAIN1100 
BAIN1110 
HAIN1120 
BAIN1130 
HAINIIUO 
HAIN1150 
BAIN1160 
BAIN1170 
HAIN11S0 
flAINl190 
aAIHl200 
flAIN1210 
HAIN1220 
BAIN123Q 
BAIN1240 
BAIN1250 
BAIN125Q 
HAIN1270 
HAIN1280 
aAIN1290 
BAIN1300 
BAIN1310 
aAINl320 
aAIN1330 
aAIN13l»0 
BAINISSO 
BAIN1360 
aAIN1370 
aRIH1380 
aAIN1390 
BAIN1430 
BAINKtlO 
HAIN1420 
BAXNIItSO 
BAINItlO 
BAIN1450 
flAIN1460 
aAIN1<«70 
BAIN11130 
BAIN1490 
HAIN1500 
BAINISIO 
aAIN1520 
UAIN1530 
UAIN1540 
BAIN1550 
HAIN1560 
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C  FBHC0013 

SOBROOTIBE    FRUCLS  FBBC0020 

C  FFMC0030 

C***********************************************************  *******  *****FB.tlCQOHO 

C  PEBC0050 

C         THIS    SOBPEOGRAH    READS    IN    AN    AREA?    OF    FOEH    CLASSES    FEOH    A  FEnC0050 

C         OSER    SPECIFIED    DISK    FILE.       THIS    DISK    FILE    IS    FREE    FOBHATTED  FRHC0070 

C         OSING    LIST    DIRECT  ED    INPUT.       THE   USEE    IS    GIVEN    THE    OPTION  FEHCOOSO 

C         TO    LIST    THE    FOBB    CLASSES    AFTER    THEY    AEE    READ.  FEMC0090 

C  FEMC0  190 

C** ***«******«>******•***********«********«****«************«***** *******FRHCO 110 

C  FRHC0120 

INTEGER    FC(32) ,SPECIE(32)  PBMC3130 

DOOBLE    PRECISION    FCFILE  FRBC0140 

REAL    fiEGCOF(33) ,BEGINT(33) ,SEBEG(35) .BSQB (33) ,HERNY(35)  FBHCOISO 

INTEGEB    aiNDBH(32) ,NTR(35) ,TOT»TR(35)  FRHC0150 

COMaON/SEG/REGCOF,BEGINT,SPECIE,HIHDBH,SEREG,NTR,  RSQR.HEANY  FRI1C0170 

I.TOTNTB  FBHC0130 

COHBON/FORMCL/FC  PRaC0190 

C    HAVE   OSER    SPECIFY    FOBH    CLASS    FILE   NABE  FBHC0230 

«EITE(5,10)  FBHC0210 

10    FOEBAT('0  ENTER    FOEH    CLASS    FILE    NAME')  FEI1C0220 

Z    BEAD  FOBM  CLASSES  FBOH  FILE  FEMC0230 

BEAD(5,20)  FCFILE  FRMC02aO 

20    FORHAT(AIO)  FRHC0250 

OPEN(FILE=FCFILE,ACCESS=«SEQIN',0HIT=1)  PRaC0  260 

READ(1,30)     FC  FRWC0270 

30    FOBHAT(8{I2,X))  FBHC0230 

C    ASK    USEE    IF    HE    WANTS    TO   LIST    FOEH    CLASSES    JUST    BEAD  FBBC0290 

UO    BEITE{5,50)  FEHC0300 

50    FORHAT('0DO    YOO    WANT    A    LIST    OF    THE    FOEH    CLASSES    JOST    READ?')  FBaC0310 

READ(5,60)     LSTANS  FBaC0320 

60    F0RHAT(A3)  FBBC0330 

IF(LSTANS.NE. ■ YES' .AND.LSTANS. NE. 'NO')            GOTO    HO  FBHC0340 

IF (LSTANS. EQ. 'NO')     EETOEN  FEHC0350 

WRITE(5,70)  (SPECI2(I) ,FC(I) ,1=1,32)  FRaC0350 

70  FOBBAT('0',8(/'  •,4(A2,'  =  ',I2,5X)))  FRaC0370 

EETOBN  FEHC0330 

END  FBHC0390 
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SOBROOTIHE    LIHREG (IPILE.TOPfiNS) 
C 

C 

THIS  SOBPBOGRAH  READS  IB  DATA  FE3H  ALL  THE  SAHPLE  TREES 

IT  FINDS  THE  TREES  HITH  HEIGHT  HEASOREMENTS 

IT  THEN  FITS  A  REGRESSION  HODEL  FOR  EACH  SPECIE  SAHPLED 


HODEL:   HEIGHT  =  BO  +  B 1  *  1/DBH 


C 

c 
c 
c 

c 
c 

DOUBLE    PRECISION    IFILE 

REAL    SEREG(35) ,RSQR(33) ,MSE,TYPSSE (3) 

REAL    VARX(33) , V ARY (33) , COVXY (33) , REGCOF  (33)  ,REGINT{33) 
REAL    X(33),Y(33) ,XS0a(33) ,YS0B(35) ,XS0HSQ(33) ,YS0HSQ(33) 
REAL   XYS0H(33)  ,  BEANX  (33)  ,  tlEANY  (35)  ,  HEIGHT,  HNYCOR  (3) 
INTEGER    SPECIE(32),NTR(35) , POINT, TREES, DBH, SPEC .HINDBH  (32) 
INTEGER    TOTNTR(35)  ,N0IiREG(3)  ,NTRC0E(3)  ,TOPANS 
COHHON/REG/REG CO F, REG INT,  SPECIE, BIN DBH,SE8EG,NTR,BSQB,HEA NY 
1,T0TNTR 
Z    SET    HINiaOB    DBH'S    TO    8 
DO    10    ISPEC=1,32 
10         aiNDBH(ISPEC) =8 


ZERO    ALL    ACCOHOLATOES 

DO  20  JTyPE=1,3 

NTRCOR(JrYPE)=0 

HNYCOR (JTYPE)=0 

NOHREG(jrYPE)=0 
20    TYPSSE(JTYPE)=0 
DO  30  ISPEC=1,35 

NTR (ISPEC) =0 ; YSOM (ISPEC) =0 

SEREG(ISPEC) =0; BEANY (ISPEC) =0 

IF(ISPEC.GT.32)   GOTO  30 

XSOa(ISPEC) =0;XSOHSQ(ISPEC)=0 

¥SOBSQ(ISPEC)=0:XYSUB(ISPEC) =0 

R EGCOF( ISPEC) =0 ;EEGINT (ISPEC) =0 ; RSQR (ISPEC)  =0 
30  CONTINUE 


INPUT  POINT  INFOEBATION 

«0  READ(21,50,END=130)  POINT, TREES 
50  FORBAT(I2,X,I2) 
CHECK  TO  SEE  IF  DATA  IS  IN  SOHHAEY  FORB 
IF(POINT. NE.-9)   GOTO  60 
TREES=1000000 
INPUT  INDIVIDUAL  TREE  DATA 
60  LOOP=0 
70  IF (LOOP. GE. TREES)   GOTO  «0 

READ{21,80, END=«0)  SPEC, DBH, HEIGHT, lEEP 
80  FORMAT(A2,X,I2,1X,F3.0,I2) 
IF(IREP.EQ.O)   IEEP=1 
LOOP=LOOP+IREP 
CHECK  TO  SEE  IP  TREE  IS  A  CULL 
IF  (HEIGHT. EQ. 999)   GOTO  70 
CHECK  DBH 

IF(DBH. LT. 1)   DBH=1 
IF(DBH.GT.«0)   DBH=<*0 
CHECK  HEIGHT 

IP(HEIGHT-LT.O.)   HEIGHT=0. 

IF (TOPAHS„ EQ. • FEET' . AND. HEIGHT.  GT. 200.  )   HEIGHT=200. 
IF (TOPANS. EQ. ' LOGS •• AND. HEIGHT. GT.  7.)   HEIGHT=7. 
CHECK  TO  SEE  IF  HEIGHT  OF  THIS  TREE  WAS  BEASURED 

IF(HEIGHr. EQ. OJ   GOTO  70 
DECODE  SPECIE  CODES 
DO  90  ISPEC=1,32 
90    IF(SPEC.EQ. SPECIE (ISPEC))   GOTO  110 
PRINT  EEEOE  BESSAGE 

WBITE(5,100)  SPEC, POINT 
100  FOEBAT(«0***ERROR***  ILLEGAL  SPECIE  CODE:  ',A2,'  AT  POINT: 
l/'ODATA  FOR  THIS  TREE  IGNORED!  EXECUTION  CONTINUES.') 
GOTO  70 
DETERBINE  HINIBOB  BEECHANTABLE  DBH 
110  IF(DBH,LT.BINDBH(ISPEC)  )   flIHDBH (IS  PEC) =DBH 


LNBG0010 
LNBG0020 
LNRG0030 

♦»»»»LBBGOO<»0 
LNRG0050 
LNRSOOSO 
LNBG0070 
LNBG0030 
LNSG0090 
LNRGOIOO 
LNBG0110 

♦•♦♦♦LB8G0120 
LKRG0130 
LNBG01U0 
LNRG0150 
LNBG01S0 
LNBG0170 
LNRG0130 
LNBG0190 
LNRG0230 
LNBG0210 
LNSG0220 
LNRG023a 
LNBG02'tO 
LKBG0259 
LNBG0260 
LNBG0270 
LNBG0280 
LNBG02J0 
LNBG0300 
LNBG0319 
LNBG0320 
LNBG0330 
LNBG0  3ltO 
LNEG0350 
LNBG0360 
LNBG0370 
LNBG0330 
LNBG0390 
LNBGOtOO 
LNBG0410 
LNBG0(»20 
LNBSOU30 
LNBG0440 
LNBG01I50 
LNRG0t60 
LNBG0470 
LNRGOaSO 
LNRG0490 
LNB60530 
LNBG0510 
LNBG0523 
LNRG0530 
LBBGOSUO 
LNBGOSSO 
LNBG0560 
LHBG0570 
LNRG0580 
LNBG0590 
LNBG0600 
LNBG0610 
LNBG0620 
LHBG0630 
LNRG06t»0 
LNRG0650 
LNRG0660 
LNRG0670 
LNBG06aO 
LNRG0690 
LNBG3700 
LNBG0710 
,12,  LNBG0720 
LNRG0730 
LBBG07ltO 
LNEG0750 
LNBG0760 
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C    COnPOTE    ACCOHOLATED   SOHS    FOB    INDIVIDOAL    SPECIE 

ITYPE=33 

IF(ISPEC.  GE.  12)        ITYPE=3IJ 

DO    120    ID0HMY=1 ,IREP 
X(ISPEC)  =  1.0/DBH 
T(ISPEC) =HEIGHT 
NTB(ISPEC)=NTB(ISPEq  +1 
BTB(ITYPE) =NTR(ITyPE) +1 
NTB{    35      )=NTE(    35      ) +1 
XSOJ!(ISPEC)  =XSOH(ISPEC)  *t  (ISPEC) 
TSO  H  (IS  PEC)  =  Y  SO  H  ( IS  PEC)  ♦Y  ( ISP  EC) 
YSDH(ITyPE) =YSOH{ITYPE)*Y  (ISPEC) 
¥SOa(    35      )=YSOH(    35      )+Y(ISPEC) 
XSOMSQ (ISPEC)  =XSOMSQ (ISPEC) +X (ISPEC) **2 
YSOHSQ(ISPEC)  =YSOaSQ(ISPEC) ^Y (ISPEC)**2 
XYSOM(ISPEC)=XYSOa(ISPEC) +X(ISPEC) *Y(ISPEC) 


120 


130 


COHTINOE 
GOTO  70 
CONTINUE 


COUPOTE  REGRESSION  PABAHETEBS  FOR  ALL  SPECIES 


DO    170   ISPEC=1,32 
JTYPE=1 

IF(ISPEC-GE.12)       JTYPE=2 
IF(NTR(ISPEC)  .LT.  1)        GOTO    170 
flEaNY(ISPEC) =YSOH (ISPEC) /HTR(ISPEC) 
C    CHECK    TO   SEE    IF    ENOOGH   TREES    EXIST    FDR    REGRESSION 

IF(NTR(ISPEC) .GE.3)       GOTO    140 
C    CALOLATE   CORRECTION    SOHS    FOR    SPECIES    WHERE    NO    REGRESSION    HAS    BUN 
HHYCOR (jrYPE) =BNYCOB (JTYPE) +YSOM (ISPEC) 
NTHCOR(JTYPE) =NTRCOR (JTYPE) +NTR (ISPEC) 
GOTO    170 
1«0         NOaREG (JTYPE) =N0aBEG  (JTYPE) +1 

aEANX (ISPEC) =XS0a (ISPEC) /NTB(ISPEC) 
VARX (ISPEC)  =  ( XSOaSQ  (ISPEC) - 

(XSOa{ISPEC)  **2/NTB (ISPEC)  )) /(NTB (ISPEC) -1) 
VARY (ISPEC)=(YSOBSQ(ISPEC) - 

(YSOa(ISPEC) ♦♦2/HTB (ISPEC) ))/(NTB (ISPEC) -1) 
COVXY (ISPEC)  =  (XYSOB (ISPEC)  - 

XS OH (ISPEC)  *YS0a (ISPEC) /NTH (ISPEC) ) /(NT R (IS  PEC) -1) 
TO   SEE    IF    VARIANCE    IS    ZERO 

IF  (VARX  (ISPEC)  -GT.  0.00000  LAND.       VARY  (ISPEC).  GT.O.OOOOOI) 
GOTO    150 

REG INT (ISPEC) =HEANY (ISPEC) 
GOTO    160 
150  REGCOF(ISPEC)=COVXY (ISPEC) /VARX (ISPEC) 

REGINT(ISPEC) =MEANY (ISPEC) -REGCOF  (ISPEC) *H£ANX (ISPEC) 
160    IF(VABY (ISPEC). LT. 0.0000001)   GOTO  170 

SSE  =  YSOHSe (ISPEC) -REGINT (ISPEC)  ♦YSOa(ISPEq - 
1      BEGCOF(ISPEC) *XYS0H (ISPEC) 

TYPSSE(jrYPE) =TYPSSE(JTYPE) ♦SSE 
aSE=SSE/(NTB(ISPEC)-2) 

SEBEG(ISPEC)=SQBT(HSE/NTB(ISPEC)  )  * 100/HEANY (I  SPEC) 
RSQR  (ISPEC) =100* (COVXY (ISPEC)  *»2) / (VA RX (ISPEC) ♦VABY (ISPEC)) 
CONTINUE 


1 


1 

CHECK 

1 


170 


COHPOTE    STANDARD   ERROR    FOB    SOFTWOOD, HAR DHOOD,    ALL    SPECIES 

NTRC0B(3) =NTRC0B(1) +NTBC0B(2) 
HNYC0E{3)  =aNYC0B(1) +aHYC0E(2) 
TYPSSE(3)  =TYPSSE(1)  ■••TYPSSE(2) 
NUHREG(3)=N0HREG(1)  +N0HBEG(2) 
DO    180    JTYPE=1,3 

ISPEC=JTYPE+32 

IF (NTR (ISPEC) .LE.O)       GOTO    180 

HEANY  (ISPEC) =YSOa(ISPEC) /NTR(ISPEC) 

CORNTB=NTR(ISPEC) -NTBCOE (JTYPE) 

IF(COBNTB.LE.  0)        GOTO    180 

COBHN¥= (YSOH(ISPEC) -HNYCOB (JTYPE) ) /CORNTB 

HSE=rYPSSE(JrYPE) /(CORNTB- 2*N0aEEG(JTYPE)  ) 

SEREG (ISPEC) =SQRT(aSE/CORNTB) ♦lOO/CORBNY 
180  CONTINUE 

CLOSE  (0BIT=21) 

0PEN(UNIT=21,ACCESS=«  SEQIH* ,FILE=IFILE) 

SKIP    RECORD    21 

RETURN 

END 


LNRG0770 
LNBG0780 
LNBG0790 
LNBG0800 
LKRGOSIO 
LNRG0820 
LNBG0B3a 
LNBG08U0 
LNBG0850 
LNRG0860 
LNRG0870 
LNRG0880 
LNRG0890 
LNRG0900 
LNRG0910 
LNRG0920 
LNR60930 
LKBG09l«0 
LNRG0950 
LNRG0960 
LNRG0970 
LNRG0990 
LNRG0990 
LNRG1030 
LNRG1010 
LNBG1020 
LNRG1030 
LNRGIOHO 
LNRG10S0 
LNRG1060 
LNRG1070 
LHRG1080 
LHRG1090 
LNRG1100 
LNRG1110 
LNRG1120 
LKRG1130 
LNRGIIitO 
LNRG1150 
LNRG1160 
LNRG1170 
LNRG1130 
LtlRG1190 
LNRG1200 
LKRG1210 
LUeG1220 
LNBG1230 
LNRG12U0 
LNEG1250 
LNBG1260 
LNRG1270 
LNRG1230 
LNRG1290 
LNRG1300 
LNRG1310 
LNBG1320 
LNBG1330 
LNRG1340 
LNBG1350 
LNRG1360 
LNfiG1370 
LNRG1380 
LNBG1390 
LNRGItOO 
LNRGIltIO 
LNBG1420 
LNRG1(t30 
LNRGIUaO 
LNRG1450 
LNBG1460 
LNRG1470 
LNRGIdSO 
LNRG1<(90 
LNRGISOO 
LNRG1510 
LNRal520 
LNRG1530 
LNRG1540 
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c 

c* 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c* 


SOBBOOTIHE   GOTS 
THIS    SOBPBOGBAH    IS    THE    'GOTS'    OF    INVEST. 


GOTS    PEBFOBaS    THE    FOLLOWING: 

1.  INPUTS    PLOT    AND    TREE    DATi 

2.  CALCOLATES    BASAL    AREA,    TREES    PER    ACRE    AND   ALL    VOLOHES 

3.  CALCULATES    VABIAHCES    AND    STAHDABD    ERRORS 

**«««***««**«4<  **•*«*«*********«*«**  4i**4i***4i*********4<«*  ********** 

DOUBLE    PRECISION    SAatlTD.IPILE 

INTEGER    REGANS, ISTART (2), IST0P(2) , DBH.SPEC, TOTNTB (35) 

INTEGER    SPECIE (32) , POINT,  TREES, TBECNT (UO, 32) 

INTEGER    CLASIZ,TITLE(14) .MOLT.TOPANS.HINDBH (32) .NrR(35) 

REAL    VABBRD(35)  ,VABP0L(35)  ,SEREG(35)  ,BSQR  (33) ,HEAN¥(35) 

REAL  BEGINT(33) ,REGCOF(33) ,SECDB(35) ,SEBRD(35) 

REAL  C0BPT(35)  ,BRDPT(35) ,P0LPT(35) ,C0BSQ(35)  ,  BRDSQ(35) 

REAL    BASPAC(U0,35)  ,  BASTOT(35)  ,  rREPAC(«»0,35)  ,TBETOT(35) 

REAL  CBVLSa(l»0,35)  ,CBTOT(35)  ,BDVLSM  (40  ,  35)  ,  BDTOT  (35)  ,  HEIGHT 

REAL  SEP0L(35) , P0LSQ(35)  ,  PLVLSH (40, 35) ,PLT0T(35) 

CO HHON/COaUNE/CLASIZ, TITLE, HOLT, PR0BLV,BAF,T3PANS, BEG ANS,A:RE 

1,FEAHE, SAHHTD, LOGANS, IPILE 
C0HH0N/REG/REGC0F,BEGINT, SPECIE,HINDBH,SEREG,NTR,RSCB,HEAHf 

I.TOTNTR 
COBHON/OOTVBASPAC,  BASrOT,BDTOT,BDVLSa,CBTOT,CBVLSa,TREPAC,TR 

1,SECnB,SEBBD,lJPHTS,PLVLSB,PLTOT,SEP0L,FINP0P,  VARBRD,VABPIJL 
DATA    (SPECIE(I) , 1=1, 32) /• HP' , • BP* , ' PP* » ' Ha« , ' BF« , 'TA* , 'SP' , • 

1,«SC«,«CE','0S« ,«SM«,«Ba«, 'HA', 'HI'.'AS'j'YP', 'BC,  'HB' 

2,'lfB','SB',  'GB'  ,'BE«,  'BH*  ,'R0«,  'HO'  ,'BO','EL',  'GO' 

3,'HI','HH', 'OH'/ 
DATA    ISTART(I) /1/,ISTART(2) /12/,IST0P(1) /1 1/,IST0P(2|  /32/ 

ZERO   ALL   ACCOaULATORS 

DO    10    ISPEC=1,35 

BASTOT(ISPEC)=0:TBET0T(ISPEC)  =0; CBTOT (ISPEC)  =0;  BDTOT (ISPEC) 

PLTOT(ISPEC) =0:TOTHTR(ISPEC)  =0 

COBSQ(ISPEC) =0    ;     BRDSQ (ISPEC) =0    ;    POLSQ  (ISPEC) =0 

SECUB(ISPEC)=0:SEBED(ISPEC) =0 ; SEP OL (ISPEC) =0 

DO    10    DBH=CLASIZ,40,CLASIZ 

TREPAC(DBH, ISPEC)  =0; BASP AC ( DBH,  ISPEC)  =0 
CBVLSB(DBH,  ISPEC)  =0;BDVLSH(DBH,  IS  PEq=0 
PLVLSB(DBH,  ISPEC)  =0 

IF (IS  PEC. LT- 33)   TRECNT (DBH.ISPEC) =0 
10  CONTINUE 
NP»TS=0 

O***** ******* *****«**************#******4>* 
BEAD  IN  ALL  TREE  DATA  AND  ACCUBULATE  SUBS 
****************************************** 

INPUT  POINT  OR  PLOT  INFOBMATION 


NS 


20    READ(21,30,END=140)     POINT, TREES 
30    FORHAT (12,13) 
CHECK    TO    SEE    IF    DATA    IS    IN    SOHBARY    FORB 
IF (POINT- NE- -9)        GOTO    40 
NPNTS=TREES 
TREES=1000000 
INPUT    INDIVIDUAL   TREE    DATA 
40    LOOP=0 
50    IF  (LOOP.  GE.  TREES)       GOTO    120 

BEAD(21,60,END=120)     S PEC, DEH, PROD, HEIGHT, IREP, PCTS AW 
60         FORHAT(A2,X,I2,A1,F3.  0,  I2,F3-2) 
IF(IREP.EQ.O)       IREP=1 
LOOP=LOOP*IEEP 
DECODE    SPECIE   CODES 
DO  70    1=1,32 
70  IP(SPEC.EQ.SPECIE(I))       GOTO   90 

PRINT    EBBOB    HESSAGE 

WRITE(5,80)     SPEC, POINT 
80    FORBAT('0***EBROB«**    ILLEGAL    SPECIE   CODE:     •,A2,'    AT    POINT: 
V'ODATA    FOR    THIS    THEE    IGNORED!    EXECUTION   CONTINUES.') 
GOTO    50 


GUTS0010 
SUTS0020 
GUTS0030 

***»*aOIS0040 
G0TS0050 
GOTS0060 
GOTS0070 
G0TS0080 
GUTS0090 
GUTS0100 
GUIS0110 
GnTS0120 
G0TS0130 

*****GUTS0140 
GOrS0150 
GUTSOI^O 
GOTS0170 
GOTSOISO 
GUTS0  190 
GUTS0200 
GnTS0210 
GUTS0220 
G0TS0230 
GaTS0240 

S  GOTS0250 
GUTS0260 
GUTS0270 
GnTS0280 

ETOI  GUTS0290 
GUIS0300 
GOTS0310 
GUTS0320 
GOTS0330 
GOTS0340 
G0TS0350 
GaTS03S0 
GUTS0  37O 
GUTS0330 
GUTS0390 
GOIS0400 
GnTS0410 
G0TS0420 
SnTS0430 
GOTS0440 
GUTS0450 
GUTS0460 
GUTS0470 
SnTS0480 
G0TS0490 
SUTS0500 
GUTS0510 
GUIS0520 
GOrS0530 
G0rS0540 
GUTS0550 
GnTS0560 
30IS0570 
G0TS0580 
GOTS0590 
GUIS0600 
SUTSOeiO 
GnTS0620 
GUTS0630 

soTsoeuo 

GOrS0650 
SOTS0660 
GUTS0670 
GnTS0680 
GUTS069a 
SUTS07aO 
G0TS0710 
GOTS0720 
GOTS0730 
GUrS07l»0 
G0TS0750 
GUTS0760 
GOrS0770 
GUTS0780 


=0 


.12, 
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C  IF 


C 

c 
c 

c 


CHECK  DBH  TO  SEE  IF  HITHIH  LEGAL  LIHITS 
90  IF(DBH-LT-CLASIZ)   DBH=CLRSIZ 

IF (DBH- GT. 40)   DBB=aO 
T&LLT  TBEES  FOB  EACH  SPECIES 

ISPEC=I 

T8ECNT (DBH,ISPEC)=TEECNT(DBH, ISPEC)  ♦IBEP 
CHECK  TO  SEE  IF  DIAHETER  IS  ABOVE  aiHIMOM  LIBIT 

IF(DBH.LT.HIND6H(ISPEq . AND. PEGANS. EQ. • NO • )       GO 
aCCOHOLATE    HOHBEE    OF    HOH-COLL    TBEES    ABOVE    BIN    DBH    S 

TOTHrE(ISPEq=TOTNTB(ISPEC)  +IBEP 

IF  (ISPEC.  LE- 11)       TOTNTE(33)  =TOTHTE(33)  -HEEP 

IF  (ISP EC.  GE.  12)        TOTNTE(3a)=TOTlITE(3a)+IBEP 

T0TNTE(35)=T0rNTE(35)  +IBEP 
HEIGHT  HAS  NOT  HEASUBED  GENEBATE  HEIGHT  FEOH  BEG 

IF(HEIGHT.EQ.O. AND. REGANS. EQ.« NO' )   HEIGHr=BEGI 

IF (HEIGHT. EQ.O)   GOTO  50 
COHPOTE  VOLOHES 
CHECK  TO  SEE  IF  TREE  A  COLL 

IF(HEIGHT.EQ.999)   GOTO  50 

HEIGHT= HEIGHT *HULT 
CHECK  HEIGHT  TO  SEE  IF  BITHIN  LEGAL  LIHITS 

IF(HEIGHT. LT. 0)   HEIGHT=0 

IF  (BOLT.  EO.  1.  AND.HEIGHT.GT.  200)   HEIGHT=200. 

IF(flOLT.EQ. 16. AND. HEIGHT. GT. 128)   HEIGHT=128 

IF  (TOPANS.EQ.  TOTAL')   CALL  VOLTOT  (DBH,  HEIGHT, I 

IF(TOPANS-EQ. 'BERCH')   CALL  VOLHER(DBH, HEIGHT, I 

IF(TOPANS„EQ.'BIXED'. AND. ISPEC- GT. 11)  CALL 

1ISPEC,C0BVOL, BDVOL) 

IF (TOPANS. EQ. ' flIXED • . AND. ISPEC.  L£. 1 1 ) 
1 1S  PEC, CO  BVOL, BDVOL) 

POLVOL=0 
CHECK  TO  SEE  IF  %    SAHLOG  OPTION  OSED 

IF(PBOD.  NE. 'X')   GOTO  100 
CHECK  PCTSAW  TO  SEE  IF  WITHIN  LEGAL  LIHITS 

IF(PCTSAH.LT-O)        PCTSAH=0.0 

IF (PCTSAW. GT. 1.0)       PCTSAW=1-0 

BDVOL  =PCTSAW*BDVOL 

P0LVOL= (1. -PCTSAW) ♦COBVOL 

C0BV0L=PCTSAW*C0BVOL 

GOTO  110 
CHECK  TO  SEE  IF  TBEE  IS  POLPWOOD 
100  IF(PBOD.  NE. 'P')   GOTO  110 

P0LVOL=C0BVOL;C0BVOL=0:BDVOL=0 
110  CBVLSHiDBH,ISPEC)=CBVLSfl(DBH,ISPEC)+C0BVOL*IEEP 

BDVLSH(DBH,  ISPEC) =BDVLSH(DBH, ISPEC) +BDVOL*IEEP 

PLVLSa(DBH, ISPEC)  =PLVLSB(DBH, ISPEC)  +PDLVOL*IBEP 


TO    50 
ABPLED 


BESSION 

NI  (I)>BEGCOF 


SPEC, COB VOL, 

SPEC, COBVOL, 

VOLHER(DBH, 


CALL    VOLTOT (DBH, 


ACCOBOLATE    SOUS    FOB    VOLOBE    STANDAED    ERRORS 


CHECK    TO    SEE    IF    DATA    IS 
IF  (POINT. EQ. -9)        GOTO    50 


IN    SOHHAEY    FOEB 


IF(FBABE-EQ. 'PLOT')       FACTOB= 1.0/BAF 

IF (FRABE. EQ. ' POINT' )       FACTOE=BAF/(. 005a5a*DBH** 

COBPT (ISPEC) =C0BPT (ISPEC)  +C0BVOL*IEEP*FACTOE 

BRDPT(ISPEC)=BBDPT(ISPEC) ♦BDVOL ♦IEEP*FACTOB 

POLPT(ISPEC)=P0LPT(ISPEC)+PULVOL»IEEP*FACTOB 

IF(ISPEC„LE.11)        ITYPE=33 

IF(ISPEC.GT.II)        IT?PE=Jl» 

COBPT (ITrPE)  =  COBPT(ITTPE)  +C0BVOL*IBEP*FACT0R 

BBDPT(ITTPE)=BBDPT(ITYPE) +BDVOL«IEEP*FACTOE 

PDLPT(ITJPE)=POLPT(ITYPE) ♦P0LVOL*IEEP*FACTOE 

COBPT (35) =C0BPr (35) ♦C0BVOL*IEEP*FACTOR 

BBDPT(35)  =BRDPT(35)  ♦BDVOL*IEEP»FACTOE 

POLPT (35) =POLPr (35) ♦P0LVOL*IEEP*FACTOR 


2) 


GOTO    50 
120    CONTINOE 
IF    IN    SOHHART    BODE    DON'T    ACCOHOLATE 

IF (POINT. EQ. -9)        GOTO    1U0 
TALLY    NOHBER    OF    POINTS    SABPLED 

NPNTS=HPNrS+1 


VOLOHES    OB   INCREBENT    NPNTS 


G0TS0790 

Gorsosoo 

GOTS081Q 
GDTS0820 
GDTS0830 
GOTSOBItO 
G0TS0850 
GOTSQBSa 
GOTS0870 
G0TS0880 
60TS0890 
60TS0900 
GnTS0910 
GOTS0920 
(I)/DBHG0rS0930 
SOTS09li0 
G0rS0950 
GOIS0960 
G0rS0970 
GOTS0980 
GOrS0990 
G0TS1000 
GOTS1010 
G0TS1020 
BDVOL)  aOTS1030 
BDVOL)  GOTS10H0 
HEIGHT, GOTS1050 
GOTS1060 
HEIGHT, GOTS 1070 
G0TS1080 
GOTS1090 
G0TS1100 
GUTS1110 
GOTS1120 
G0TS1130 

GOTsma 

GOTS1150 
GQTSIIBO 
G0rS1170 
GOTS1130 
G0TS1190 
S0TS1200 
GOTS1210 
G0TS1220 
GOrS1230 
GOrS12U0 
G0TS1250 
GDTS1260 
G0rS1270 
GaTS1280 
G0TS1290 

sorsi300 

GOTS1310 
G0TS1320 
GOTS1330 
G0TS1340 
G0TS1350 

GOTsneo 

G0TS1370 
GDTS1330 
G0TS1390 
5UTS1400 
GOTSiaiO 
G0TS1U20 
GOTS1U30 
GOIS1440 
GOTS1450 
50151460 

Gorsia70 

G0TS1430 
G0TS1490 
GOTS1500 
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ACCOHOLATE    SOHS    OP    SQOAEES    FOB    VOLOHES 
DO    130    ISPEC=1,35 

COBSQ(ISPEC)  =COBSQ(ISPEC)  -fCOBPT  (ISPEC)  **2 
BEDSQ(ISPEC) =BRDSQ{ISPEC)  +BEDPT (ISPEC)  **2 
POLSQ(ISPEC) =POLSQ(ISPEC) +POLPT (ISPEC)**2 
COBPT  (ISPEC) =0 
POI.PT(ISPEC)  =0 
130         BBDPT (ISPEC)  =0 

GOTO  20 
1U0  CONTINOE 

IF  ONLY  ONE  POI HT  SAHPLED,  SET  IN  SOHHARY  BODE 
IP(HPNTS. EQ. 1)   P0INT=-9 

DETEBMINE  SAHPLING  METHOD 
***««******«***«*««  ««>*•  ******* 


IF(FRAHE-EQ.'PLOT«)   GOTO   170 

*************************************************** 
POINT    SAHPLING 

CALCULATE    PER    ACRE    ESTIMATES    FOR    IHDIVIDOAL    SPECIES 
*************************************************** 


DO    160    ISPEC=1,32 

DO    150    DBH=CLASIZ,l»0,CLASIZ 
C    CHECK    TO    SEE    IF    ANY    TREES    TALLIED 

IF  (TEECNICDBH,  ISPEC)  .EQ.O)       GOTO    150 
C    BASAL    AREA    CALCULATIONS 

BASPAC(DBH, ISPEC) =TRECNT (DBH, ISPEC) 'BAF/NPNTS 
BASrOT (ISPEC) =BASTOT( ISPEC)  +BASP AC (DB H, ISPEC) 
C    TREES    PEB    ACRE   CALCULATIONS 

TREPAC(DBH,  ISPEC)  =BASPAC  (DBH,  ISPEC)  /(.  005l*5a*DBH**2) 
TBETOT (I SPEC) =TRETOT (ISPEC)  +TREPAC (DBH, ISPEC) 
C   CUBIC    VOLUME    CALCULATIONS 

CBVLSH (DBH, ISPEC) =CBVLSM (DBH, ISPEC)  *BAF/ 
1  (NPNTS*.005«5<»*DBH**2) 

CBTOT (ISPEC) =CBTOT( ISPEC)  +CBVLSM (DBH, ISPEC) 
C    BOARD    FOOT    VOLUME    CALCULATIONS 

BDVLSM(DBH, ISPEC)  =BDVLSa (DBH,  ISPEC)  ♦BAF/ 
1  (NPNTS*.005«54*DBH**2) 

BDTOT (ISPEC) =BDTOT (ISPEC) +BDVLSH (DBH, ISPEC) 
C    PULPWOOD    CALCULATIONS 

PLVLSM(DBH, ISPEC) =PLVLSfl(DBH, ISPEC)  *BAF/ 
1  (NPNTS*.005a5tt*DBH**2) 

PLTOT  (I  SPEC) =PLTOT (ISPEC)  ♦PLVLSM (DBH, ISPEC) 
150         CONTINOE 
160    CONTINOE 
GOTO    200 


*************************************************** 
PLOT    SAMPLING 

CALCULATE    PER    ACRE    ESTIMATES    FOR    INDIVIDUAL   SPECIES 
*************************************************** 


170    PLTSIZ=BAF 

NPLOTS=NPNTS 

DO    190    ISPEC=1,32 

DO    180    DBH=CLASIZ,l»0,CLASIZ 

CHECK    TO    SEE    IF    ANY    TBEES    TALLIED 
IF(TRECNr  (DBH,  ISPEC).  EQ-0) 

TREES    PEB    ACRE   CALCULATIONS 

TBEPAC(DBH, ISPEC) =TRECNT(DB 
TRETOT (ISPEC) =TRETOT(ISPEC) 

BASAL    AREA    CALCULATIONS 

BASPAC(DBH, ISPEQ  =TEEPAC(DB 
BASTOT (ISPEC) =BASTOT(ISPEC) 

CUBIC    VOLUME    CALCULATIONS 

CBVLSM(DBH, ISPEC)  =CBVLSH(DB 
CBTOT (I  SPEC) =CBTOT (ISPEC) +C 

BOARD    FOOT    VOLUME    CALCULATIONS 

BDVLSH(DBH, ISPEC) =BDVLSM(DB 
BDTOT  (ISPEC) =BDTOT( ISPEC)  +B 

PULPHOOD  CALCULATIONS 

PLVLSM (DBH, ISPEC)  =PLVLSM (DB 
PLTOT (I  SPEC) =PLTOT (ISPEC) +P 

183    CONTINOE 

190  CONTINUE 


GOTO  180 

H,ISPEC) /PLTSIZ/NPLOTS 
+TREPAC (DBH, ISPEC) 

H,  ISPEC)  *.0054  5'»*DBH*»2 
+BASPAC(DBH, ISPEC) 

H, ISPEC) /PLTSIZ/NPLOTS 
BVLSH(DBH,  ISPEC) 

H, ISPEC) /PLTSIZ/NPLOTS 
DVLSM(DBH,ISPEC) 

H, ISPEC) /PLTSIZ/NPLOTS 
LVLSM(DBH, ISPEC) 


GOTS1510 
GUTS1520 
GUTS1530 

Gursisuo 

GUTS1550 
G0TS1560 
G0TS1570 
GUTSISBO 
G0TS1590 
GOTS1600 

Gursieio 

G0rS1620 
GOTS1630 
G0rS16l»0 
GOrS1650 

Guisieeo 

G0TS1670 
GUrsi680 
G0TS1690 
GOrS1700 
G0rS1710 
GOrS1720 
G0TS1730 
GUrS17l»0 
GUTS1750 
GUrS1760 
GUrS1770 
G0rS1780 
GOrS1790 
GOTS1800 
G0TS1810 
GOTS1820 
GUrS1830 
GOTS18U0 
GOTS1850 
G0rS1860 
G0TS1870 
G0TS1880 
GUTS1890 
G0TS1900 
G0TS1910 
GUTS1920 
GUTS1930 
GOrS19U0 
GUTS1950 
GUrS1960 
G0rS1970 
GUrS1980 
GUTS1990 
GUTS2000 
GUTS2010 
6UrS2020 
GUTS2030 
GniS20lt0 
G0rS2050 
G0TS2060 
GUTS2070 
G0IS2080 
GUTS2090 
GUrS2100 
GUTS2110 
GUTS2120 
GOTS2130 

Gors2no 

G0TS2150 
GUTS2160 
GUrS2170 
GUTS2180 
GOTS2190 
GOTS2200 
GUTS2210 
G0TS2220 
GOrS2230 
G0TS22a0 
GUTS2250 
G0rS2260 
GUTS2270 
GUrS2280 
3UTS2290 
GUIS2300 
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CALCULATE    TOTALS    FOB    SOFTWOODS    ARD    HARDWOODS 
^^^^^^*m*************************** ********** 


200    DO   230    ITIPE=33,34 

DO    220    DBH=CLASIZ,UO,CLASIZ 

DO    210    ISPEC=ISTABT(ITYPE-32|  ,ISTOP 

BASPAC(DBH,ITyPE) =BASPAC( DBH, ITYP 

TBEPAC(DBB,  ITYPE)=TEEPAC(DBH,ITIP 

CBVLSH(DBH, ITSPE)  =CBVLSH ( DBH, ITYP 

BDVLSH(DBH,ITyPE)=BDVLSIl(DBH,ITYP 

210  PLVLSH(DBH, ITtPE) =PLVLSH (DBH, ITYP 

BASTOT(ITYPE)=BASTOT{ITYPE)  ♦BASPACC 

TEETOT(ITYPE)=TRETOT(ITYPE)  +TEEPAC( 

CBTOT  (ITYPE)  =CBTOT  (ITYPE)  -t-CBVLSH  (DB 

BDTOT(ITYPE)=BDTOT(ITYPE)+BDVLSH(DB 

223  PLTOT (ITYPE) =PLTOT( ITYPE)  ♦PLVLSa{DB 


(ITYPE-32) 

E)  +BASPAC(DBH,  ISPEC) 

E)  ♦TREPaC(DBH, ISPEC) 

E)  ♦CBVLSH(DBH, ISPEC) 

E)  +BDVLSa(DBH,  ISPEC) 

E)+PLVLSa(DBH,ISPEC) 

DBH, ITYPE) 

DBH,  ITYPE) 

H, ITYPE) 

H, ITYPE) 

H, ITYPE) 


230  CONTINOE 

*4t««4i*  «*«*««*««*«««*»****«*  ****** 

CALCULATE    TOTALS    FOE    ALL    SPECIES 
********************************* 

DO    2«0    DBH=CLASIZ,40,CLASIZ 

BASPAC(DBH,35)=BASPAC (DBH, 33)  ♦BASPiC (DBH, 34) 
TBEPAC(DBH,  35)=TBEPAC(DBH,  33)  +TEEPAC  (DB  H,  34) 
CBVLSa(DBH,35)=CBVLSH(DBH,33) +CBVLSH (DBH, 34) 
BDVLSH(DBH, 35)=BDVLSH(DBH, 33)  +BDVLSH (DBH, 34) 
PLVLSa(DBH,35)=PLVLSH(DBH,33)*PLVLSB(DBH,34) 
240    CONTINUE 

BASTOT(35) =BASTOT (33) ♦BASTOT(34 ) 

TBETOT (35) =TEErOT (33) +TBETOT (34) 

CBTOT(35) =CBT0T(33)  +CBTOT(34) 

BDTOT(35)=BDT0T  (33)  ♦BDT0T(34) 

PLT0T(35) =PLT0T{33) +PLT0T (34) 


*************************** 

CALCULATE    STANDARD    ERRORS 
*************************** 

FINPOP=1.0 

IF(FRAHE.EQ.' PLOT')       FINPOP= ( (ACRES/BAF) -NPNTS) /( ACEES/BAF) 
IF(FINPOP. LT. 0.0)        FINPOP=0.0 
CHECK    TO    SEE    IF    DATA    IS    IN    SUMHAEY    FOEB 
IF(POINT.  EQ. -9)        GOTO    270 
DO    260    ISPEC=1,35 

IF(TRETOT(ISPEC) .Eg.O)       GOTO    260 
IF(FBAHE. EQ. • PLOT')        GOTO    250 

PCTCEZ=100*NPNTS*BASTOT(ISPEC) /TRET OT (ISPEC) /ACEES/BAF 
IF(PCTCBZ-GT.  100.0)       PCTCRZ  =  100.0 
FINPOP= (100-PCTCRZ)  /100 
250         VAECnB= (CUBSQ (ISPEC) -NPNTS* (CBTOT (ISPEC) **2) ) /(NPNTS-I) 

VAEBBD(ISPEC)  =  (BEDSQ(ISPEC)-NPNTS*  (BDTOT  (ISPEC)  **2)  )  /(HPNTS-1) 
VAEPUL (ISPEC)  =  (PULSQ (IS PEC) -NPNTS *1PLT0T (ISPEC) •*2)  )/(NPNTS-1) 
SECOB(ISPEC) =SQBT((VAECUB/BPNTS)»FINPOP) 
SEBED(ISPEC)=SQET((VAEBBD(ISPEC)  /NPNTS)  *FINPOP) 
260         SEPUL(ISPEC) =SQRT ( (VARPOL(ISPEC) /NPNTS) *PINPOP) 
CHECK    TO   SEE    IF    REGRESSION    ROUTINE    USED 
270    IF(EEGANS.EQ. 'YES')        RETURN 
DO    280   ISPEC=1,35 

IF(rOTNTE(ISPEC)  .  EQ.O)        GOTO    280 

EEGHGT  =  PLOAT(TOTNTR (ISPEC) -NTR(IS  PEC) ) /FLOAT (TOTNTR (ISPEC) ) 
SECUB(ISPEC)=SQRT(SEC0B(ISPEC)»*2    * 
1  REGHGT»(SEBEG (ISPEC)  »CBTOT(ISPEC)  /100.)  **2) 

SEBRD(ISPEC)  =SQRT(SEBRD(ISPEC)**2    ♦ 
1  REGWGT* (SEEEG(ISPEC) *BDTOT( ISPEC)  /I  00.)  **2) 

SEPOL (ISPEC) =SQRT(SEPUL(ISPEC) **2    * 
1       EEGWGT*(SEREG(ISPEC)*PLTOT(ISPEC)/100.) ♦♦2) 
280  CONTINUE 
RETURN 
END 


GUTS2310 
GBTS2320 
GUrS2330 
GUTS2340 
G0rS2350 
GaTS2360 
GUr32370 
GnTS2330 
GUrS2390 
GUTS2400 
GOTS2410 
GUrS2420 
GUTS2430 
GnTS2440 
G0rS2450 
GUTS24S0 
GnTS2470 
GUTS24B0 
GUTS2490 
GUTS2530 
GUrS2510 
GUTS2520 
GUTS2530 
G0TS2540 
GOrS2553 
GOTS2560 
GUTS2570 
G0rS2580 
GnTS2590 
GUTS2600 
GUTS2610 
GUTS2620 
GOrS2630 
GUTS2640 
GUrS2650 
GUrS2660 
GUrS2670 
G0TS2680 
GUTS2690 
Gnrs2700 
GUTS2710 
GUrS2720 
GUTS2730 
GnTS2740 
G0TS2750 
GUTS2760 
GnTS2770 
GUrS2780 
GUTS2790 
GUT528a0 
GUTS2810 
GOTS2820 
G0rS2830 
G0TS2840 
GOTS2850 
GUTS2860 
SOrS2870 
GnTS2880 
GnTS2890 
GUTS2900 
GOTS2910 
GUTS2920 
G0TS2930 
G0TS2940 
G0TS2950 
G0rS2960 
GUrS2970 
GOrS2980 
GUTS2990 
GOTS3030 
GUTS3010 
GUTS3020 


-75- 


C  VLMROOIO 

SOBBOOTINE  VOLHEB (DBH, HEIGHT, ISPEC, COBVOL.BDVOL)  VLHR0020 

C  VLHR0030 

C  VLBB9050 

C      THIS    SOBPEOGBAM    COHPOTES    THE    HERCHANTABLE    COBIC    AND    BOARD    F30r  VLHR0060 

C      VOLOHES    OF    INDIVIDOAL    TREES.       THE    INPOT    PARAHETEES    ABE    DBH,  VLME0070 

C      MERCHANTABLE    HEIGHT,    AND    SPECIES.  VLMBOOSO 

C  VLHB0090 
C      THE   COBIC    FOOT    EQUATIONS    ARE    FROM    BAENABD,    JOSEPH    E.     ET    AL.     1973.            VLHE0100 

C      FOBEST   SOEVET    COBIC-FOOT    VOLOHE    EQUATIONS.    OSDA    FOBEST    SEBV.  VLHB0110 

C      NE.    FOREST    EXP.    STA.    RESEARCH    NOTE    NE-66.  VLnR0120 

C  VL8R0130 
C  THE  BOARD  FOOT  EQUATIONS  ABE  FBOH  WIANT,  H.V.  AND  CASTENADA,  P.  1977  VLHEOIUO 
C      MESAVAGE    AND    GIRAED'S    VOLUME    TABLES    FORMULATED.       BESOURCE    INVENTOBJ      VLME0150 

C      NOTES,    USDI,    BLB,     DENVEE,    COLO.  VLMR0160 

C  VLHB0170 
C«*«* ♦♦*♦♦»*♦♦♦♦*****♦***♦****♦***♦*»*****♦♦****♦***♦*♦♦*♦♦****♦♦♦♦**♦♦» VLB EO 180 

C  VLHE0190 

REAL    HEIGHT, C0BIBT(17) ,C0BSLP(17)  ,BDHT  VLMEO2D0 

INTEGEB    DBH,GBP(32) ,D,FC(32)  VLMR0210 

COBBON/FOEBCL/FC  VLaR0220 

C  VLBR0230 

C    DEFINE    FORB    CLASS    FOE    THE   SPECIES  VLBE02U3 

C  VLMR0250 

DATA    (FC(I)  ,I=1,32)/3*76,70,7*76,21*7«»/  VLHE0260 

C  VLBB0270 

C    DEFINE    VOLUME    EQUATION    GBOUPING    FOB    THE    SPECIES  VLaB0280 

C  VLBB0290 

DATA    (GEP(I) , 1= 1,32) /I, 1, 5, U, 3, 5, 2, 5, 5, 6, 1,7, 8,  9,  17,9,8,10  VLBR0  300 

1,4*11,12,13,14, 15,14, 15,1  a, 16, 2*17/  VLMB0310 

C  VLaB0320 

C    DEFINE    THE    INTEBCEPTS    FOE    THE    COBIC    VOLOBE    EQUATIONS  VLSB0330 

C  VLBB0340 

DATA    (CUBINT(I)  ,1=1 , 1 7) /3. 5142,2. 1998 , 1 . 4793, 2. 4784,-. 0496, 1.5817    VLSR0  350 

1,1.6823,1. 1763,1.0067,1.1809, 1. 1339, 1-2  851 ,. 8976, 1 .2027,1  . 0009  VLBR0  35  3 

2, 1.4438.. 8591/  VLHR0370 

C  VLBR0380 

C  DEFINE  THE  SLOPES  FOR  THE  COBIC  VOLOBE  EQUATIONS  VLMR0390 

C  VLBR34DD 

DATA    (COBSLP(I) ,1=1 ,1 7) /236 , 257 , 272, 242,3 03, 259, 310, 310,293 ,2 92         VLBR0410 

1,281,326,349,3  01,300,316,309/  VLaR0420 

C  VLBB0433 

C    FONCTION    TO    COBPOTE    BOARD    FOOT    VOLOBES    FOR    G. F. C    78  VLaR0440 

C  VLBR0450 

HIANT(DBH,BDHT)=(-13.35212*9.58615»BDHT+1.52968*BDHT**2)+  VLaR0460 

1  (1. 79620-2. 59995*BDHT-.27465*BDHT**2)*DBH+  VLHR0470 

2  (.04482+.45997*BDHT-.00961*BDHT**2) *DBH**2  VLBB0480 
C  VLaE0490 

C    DETERBINE    BOABD    FOOT    HEIGHT    TO   THE    HALF-LOG  VLMB0530 

C  VLBEOSIO 

BDHT=HEIGHT/16-0  VLBE0520 

CBHT=HEIGHT  VLaE0530 

C  VLBB0540 

C    CALCULATE    CUBIC    AND    BOARD    FOOT    VOLUBE  VLBR0550 

C  VLMR0560 
C0BVOL=COBINT(GEP(ISPEC))  +.  0000 1 ♦COBS LP (GBP (ISPEC) ) ♦CBHr*DBR**2         VLMR0570 

BDVOL=HIANr(DBH, BDHT)*(((FC(ISPEC)-78)*.03)+1)  VLaR0580 

IF(BDHT.LT. 0.5)        BDVOL=0  VLMR0590 

IF{CUBVOL-LT.O)       C0BVOL=0  VLaR0600 

IF(BDVOL. LT.O)        BDVOL=0  VLBE0610 

HETOBN  VLBB0620 

EUO  VLaB0630 
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C  VLTT0013 

SOBEOOTINE    VOLTOT  (DBH,  HEIviHT.ISPEC.COBVOL,  BDVOL)  VLTT0020 

C  VLTT0030 

C ******************************************************** ************«**VLTT0O(tO 

C  VLTT0050 

C      THIS    SOBPBOGBAH    COHPOTES    THE    HEECHANTftBLE    COBIC    AND    BOABD    FOOT  VLTTOOeO 

C      VOLOHES    OF    INDIVIDUAL    TE EES.       THE    INPDT    PAEABETEES    ARE    DBH,  VLTTOOTO 

C      TDTAL    HEIGHT,    AND    SPECIES.  VLTTOOaO 

C  VLTT0090 

C      ALL    OF    THE    EQUATIONS    OSED    ABE    FEOM    HONEB,    T.  G.     1967.    STANDAED    VOLOBE    VLTT0100 

C      TABLES    AND    BEECHANTABLE    CONVEBSIOH    FACTOES    FOE    THE    COHBEECIAL    TREE  VLTTOIIO 

C      SPECIES    OF    CENTEAL    AND    EASTEBN    CANADA.    FHB-X-5,    FBI,    OTTAWA,    CANADA  VLTT0120 

C  VLrT0130 

C**** ******************************************************* ************VLrTO lltO 

C  VLTT0150 

INTEGEE   GBP(32) ,DBH  VLTTOieO 

REAL   C0BINr(16)  ,C0BSLP(16) ,ASDBCF(3) ,ASDBBF(3) ,BEDT0P(2) .COBTOP {2 ) VLTTO  170 

C  VLTT0180 

C    DEFINE    VOLOBE    EQOATION    GBOOPING    FOB    THE    TBEE    SPECIES  VLTT0  190 

C  VLTT0200 

DATA  (GBP (I), 1=1, 32) /I, 2, 2, 6, a, 2, 3, 3, 2, 5, 2, 10, 10, 16, 15,7,  10, 1 3  VLTT0210 

1,8,9,8,8, 12, 11, 3*16,  H», 16,16,15,15/  VLTT0220 

c  vLrro230 

C  DEFINE  THE  INTEBCEPTS  FOB  THE  TOTAL  COBIC  FT  VOL  EQOATIOH  VLTT0240 

C  VLTT0250 

DATA    {CO BINT (I) , 1= 1, 16) /- 691 , . 710, 1.226,2.  139 ,4. 167, 1.  112,-.  3 12  VLTT0260 

1,2.222,1.««9,1-0i»6,  .9H8,.959,-033,.63i»,1.877,1.512/  VLrT0270 

C  VLrT0280 

C    DEFINE   THE   SLOPES    FOB   THE   TOTAL    COBIC    FT    VOL    EQUATION  VLTT0290 

C  VLTT0300 

DATA    (COBSLP(I) , 1= 1, 16) /363- 676,355.623 ,315. 832,30 1. 63U, 2au. 9  06  VLTT0310 

1,350.09  2,iJ36.683,300-373, 3ai*.754,383-972,<»01.i»56,33i».829  VLTT0320 

2, 393. 366,  4i»0.<»96, 332.  585, 336-509/  VLTT0330 

C  VLrT03i»0 

C    DEFINE    STOMP    HEIGHT    AND    BEECH    TOPS  VLTT0350 

C  VLTT0  350 

DATA  STOBP/0. 5/  VLTT0370 

DATA  (CaBTOP(I) ,I=1,2)/a.0,U.0/  VLTT0380 

DATA  (BEDTOP(I) ,1=1,2) /6.0, 8.0/  VLrT0390 

C  VLTTOUOO 

C  DEFINE  ADJUSTED  SQED  DIA  EATIO  COBIC  FT  VOL  CONVERSION  COEFFICIENTS  VLTTOIIO 

C  VLTT0420 

DATA  (ASDECF(I) ,1=1, 3)/.960«,-. 1660, -.7868/  VLTT0H30 

C  VLTTOUUO 

C    DEFINE    ADJUSTED   SQBD    DIA    EATIO    BOAE D   FT    VOL    CONVERSION    COEFFICIENTS  VLTTOttSO 

C  VLTT0460 

DATA  (ASDSBF(I)  ,1=1, 3)/5.«332,- 1.6281,-4.4710/  VLTT0470 

C  VLTTOUSO 

C  DETERMINE  SB  OB  HW  TYPE  VLTT0490 

C  VLTT0500 

IF(ISPEC.LE.II)        ITIPE=1  VLTTOSIO 

IF(ISPEC.  GE.  12)       ITYPE=2  VLTr0520 

c  vLrro530 

C    COBPOTE   TOTAL    COBIC    FOOT    VOLOBE  VLrT0540 

C  VLTTOSSO 

TCBVOL=DBH**2/(C0BINT(GBP(ISPEC))  ♦COBSLP (GBP (ISPEC) ) /HEISHT)  VLTT0560 

C  VLrr0570 

C  COBPOTE  BEECH  COBIC  AND  BOABD  FT  VOLUMES  VLTT0580 

C  VLTr0590 

XCUB=( (CUBTOP(ITYPE) /DBH) **2) * (1. O+STOMP/HEIGHT)  VLrT0600 

XBE&=  ((BEDrOP(ITYPE)/DBH)  **2)  *(  1.  0  +  STOMP/HEI3HT)  VLTT0610 

C0BVOL=TCBVOL*(ASDRCF(1)+ASDECF(2) *XCOB+ASDECF (3)  *XC0B**2)  VLTT0620 

BDVOL  =TCBV0L*(ASDBBF(1) ♦ASDEBF (2) •XBED+ASDRBF (3)  »XBBD**2)  VLrT0630 

IF(COBVOL.LT. 0)   C0BVOL=0  VLTT0640 

IF(BDVOL.LT.O)   BDVOL=0  VLTT0650 

BETOBN  VLTT0660 

END  VLrT0670 
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SOBROUTINE  STRAT (ACRES, CLASIZ) 
C 

Q******1f^***********i***t  *************  ******************m**^t.t.ttlr*m****^ti 
C 

C    THIS  SOBPROGEAH  CALCULATES  COMBINED  STRATA  STATISTICS 
C 
C********************************************************************** 

c 

INTEGER    CLASIZ 

REAL    STSECB(35) ,STSEBD(35) .EFDFCB (35) ,EFDPBD(35) 

REAL   SrVHCB(35) ,STVHBD(35)  , DFCB2(35) ,DFBD2(35) 

REAL    STV11PL(35)  ,DFPL2(35)  ,PLTST(35)  ,STSEPL(35)  ,  EFDFPL  (35) 

REAL    PLVLST(l»0,35) 

REAL  CBTST(35)  ,  BDTST(35)  ,  CBVLST  (1*0,  35)  ,  BDVLST  (40,  35) 

COHMON/O0T1/BASPAC(«»0,35)  ,BASTOT(35)  ,BDTOT(35)  ,  BDVLSH  (40,  35)  , 
1CBTOT(35) ,CBVLSn(40,35) , TREP AC (40, 35) , TRETOT(35) 
2,SECOE(35) ,SEBRD(35) , NPNTS, PLVLSH (40, 35) ,PLT0T(35) ,SEP0L(35) 
3,FISPOP,VARBED(35) ,VARPOL(35) 

COHMON/STOOT/CBTST, BDTST,CBVLSr, BDVLST, STSECB ,STSEBD, EFDFCB 
1 , EFDFBD, TA:RES, PLTST, PL VLST ,STS  EPL, EFDFPL 
C    ZERO    ACCOBOLATORS    AND   OTHER    VARIABLES 

DATA    ((CBVLST(J,I) ,BDVLST(J,I) ,PLVLST(J,I)     ,    J=1,40,1) 
I.CBTST  (I)  ,BDTST (I) ,PLTST(I) ,STVHCB(I) ,STVHBD(I)  ,STVnPL(I) 
2,DFCB2  (I) ,DFBD2  (I) ,DFPL2(I)  , EFDFCB (I) ,EFDFBD(I)  , EFDFPL  (I) 
3,STSECB(I) ,STSEBD(I) ,STSEPL(I)     ,    1=1,35)    /    4725»0.0    / 


TACRES=TACRES+ ACRES 
C 

C  ACCOaOLATE  SOMS  FOR  VOLOHES 
C 

DO    20   ISPEC=1,35 

00    10    DBH=CLASIZ, 40, CLASIZ 

CBVLST (DBH, ISPEC) =CBVLST (DBH, ISPEC) +CBVLSfl (DBH,ISPEC)  ♦ACRES 
BDVLST(DBH,ISPEC)=BDVLST(DBH,ISPEC) +BDVLSa ( DBH,  ISPEC)  *ACRES 
10  PLVLST{DBH, ISPEC)  =PLVLST (DBH, ISPEC)  +PLVLSB (DBH, ISPEC)  ♦ACRES 

CBTST (ISPEC) =CBTST(ISPEC)  ♦CBTOT (ISPEC)  ♦ACRES 
BDTST (ISPEC) =BDTST(ISPEC) ♦BDTOT (ISPEC) ♦ACRES 
PLTST  (ISPEC)  =PLTST(ISPEC)  -fPLTOT  (ISPEC)  ♦ACRES 
C 

C    ACCOnOLATE   SOHS    FOE    VARIANCE    AND    D.F. 
C 

IF(NPNTS.  LT.2)       GOTO    20 

ST VHCB (ISPEC) =STVBCB (ISPEC) ♦ (ACRES ^^2) ♦ (SEC OB (I  SPEC) ^^2) 
STVHBD ( ISPEC) =STVBBD (ISPEC) +(ACRES^^2)^(S EBRD (ISP EC) ^^2) 
STVHPL (IS  PEC) =STVBPL( ISPEC) ♦(ACRES^^2)^ (SEP OL (I  SPEC) **2) 
DFCB2(ISPEC) =DFCB2(ISPEC)  + ACEES^^4^SEC0B(ISPEC) ♦♦4/(NPNTS-1 ) 
DFBD2(ISPEC) =DFBD2 (ISPEC) ♦ACRES^^4^SEBRD(ISPEC) ^^4/ (NPNTS- 1) 
DFPL2 (ISPEC) =DFPL2(ISPEq  ♦ACBES^^4^SEPaL (ISPEC) ♦♦4/(NPNTS-1 ) 
20    CONTINUE 
RETURN 
C 

C ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦♦♦♦♦♦♦♦♦♦♦♦• 
C    CALCULATE    STANDAED    EEROBS    AND    EFFECTIVE    D.F. 
C      BASED    OH    PREVIOUSLY    ACCOHOLATED    SUflS 
C*************** *********************************** 
C 

ENTRY    ST  EAT  2 

DO    50    ISPEC=1,35 

IF(PLTST(ISPEC) .EQ.O)       GOTO   30 
STSEPL(ISPEC)=SQET(SrVBPL(ISPEC)  /TACSES^^2) 
IF(DFPL2(ISPEC)-EQ.0)        GOTO    30 
EFDFPL(ISPEC) =STVBPL (ISPEC) ♦♦2/DFPL2 (ISPEC) 
30         IF(CBTST(ISPEC)  .EQ.O)        GOTO    50 

STSECB (ISPEC) =SQET(STVBCB (ISPEC) /TACRES^^2) 
STSEBD(ISPEC) =SQBT (STVBBD (ISP  EC) /TACRES^^2) 
IF(DFCB2(ISPEC) .EQ.O)       GOTO    40 
EFDFCB(ISPEC) =STVBCB (ISPEC) ♦♦2/DFCB2 (ISPEC) 
40         IF (DFBD2 (ISPEC). EQ.O)        GOTO    50 

EFDFBD(ISPEC) =STVaBD (ISPEC) ♦♦2/DFBD2 (IS  PEC) 
50    CONTINUE 
RETURN 
END 


STET0010 
STRr0020 
STRT0030 

♦STRT0040 
STRT0050 
STBT0050 
STRT0070 

♦STET0080 
STRT0090 
STRT0130 
STRT0110 
STRT0120 
STRT0130 
STRT0140 
STRT0150 
STRT0160 
STBT0170 
STBT0180 
STRT0190 
STRT0230 
STRT0210 
STRT0220 
STRT0230 
STRT0240 
STET0250 
STRT0260 
STRT0270 
STRT0280 
STRT0290 
STRT0300 
STRT0310 
STRT0320 
STRT0330 
STRT0340 
STBT0350 
STBr0360 
STRT0370 
STRT0380 
STRT0390 
STRT0400 
STRT0410 
STRT0420 
STRT0430 
STRT0440 
STBT0450 
STBT0460 
STRT0470 
STRT0480 
STRT0490 
STRT0500 
STBr0510 
STRT0520 
STRT0530 
STRT0540 
STRT0550 
STRT0560 
STRT0570 
STBT0580 
STRT0590 
STRT0600 
STRT06  10 
STRT0620 
STRT0630 
STRT0640 
STBT0650 
STET0660 
STRT0670 
STRT0680 
STRT0690 
STRT0700 
STRT0710 
STRT0720 
STBT0730 
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C  OTP10010 

SOBEOOTIHE   O0TP1  OTP  10020 

C  OTP10030 
C*****  **************************************************************  **$i**oip  10040 

C  OTP10050 

C         THIS    SOBPROGRAH    PRODOCES    fiLL    OF    INVENTS    TABLES    AND    SOHHAFIES  OTP10060 

C  0TP10070 

C         THE    SOBPROGRAH    HAS    TWO    ENTRY    POINTS:       OOTPI         AND         00TP2  OTP10080 

C  OTP  10090 

C                 OOTPI    PRODOCES    ALL    OF   THE    INDIVIDOAL    STRATA    OOTPOT  OTP10130 

C                 00TP2    PRODUCES    THE    STRATIFIED    TOTAL    OOTPOT  OTP10110 

C  0TP10120 

INTEGER    T0TACB,T0TABD,T0TAPL,IDATE(2)  ,TOTNTR(35)  .TOTPTS  0TP10U0 

INTEGER    DBH, TOPANS, HOLT, REGANS, FC(32J  OTP10150 

INTEGER    SPETIT(35,3) ,TITLE(ia),CLASIZ,SPECIE(32)  0TP10160 

INTEGER    HINDBH(32) ,NTR(35)  OTP10170 

REAL    VARBRD(35) ,VARP0L(35) ,REGC0F(33) ,REGINT(33) ,SEREG(35)  OTP10180 

REAL    RSQR(33)  ,HEANY(35)  OTP10190 

REAL    TREPAC(<»0,35),TRETOT(35)  ,CBT0T(35)  ,BDT0T(35)  ,HSD  OTP10200 

REAL    SEC0B(35)  ,SEBED(35)  ,SEPUL(35)  ,  PLVLSB(<»0  ,  35)  ,PLTOT(35)  OTP10210 

REAL    CBVLSH(40,35) ,BDVLSH (40 ,35) ,BASPAC  (40, 35) , BASTOT (35)  OTP  10220 

DOOBLE    PRECISION    IFILE,SAHHTD  OTP10230 

COMHON/REG/EEGCOF, REGIHT, SPECIE, MI NDBH,SEREG, NT R,RSQR,HEANr  OTP 10  240 

1,T0TNTR  OTP10250 

COHHON/FORflCL/FC  OTP  10260 

COHHON/STO0T/CBTST(35) ,BDTST(35) , CBVLST (40,35) , BD VLST(40,35(  OTP10270 

1,STSECB(35) ,STSEBD(35) ,EFDFCB(35) ,EFDFBD(35) ,TACRES  OTP10280 

2,PLTST(35) ,PLVLST(40, 35),STSEPL (35) ,EFDFPL(35)  OTP  10290 

COBHON/COHUNE/CL AS IZ, TITLE, nOLT,P80BLV,BAF,TOPA US, REGANS, ACRES  OTP10300 

1,FRAaE,SA[lMTD,L0GANS,IFILE  OTP  10310 

C0HH0N/0UT1/BASPAC, BASrOT,BDTOT,BDVLSH,CBT0T,CBVLSH,TREPAC,TRETor    OTP10  320 

1,SECUB,SEBRD,SPNTS,PLVLSn,PLTOT,SEP0L,FINPOP,VARBRD,VARP0L  OTP  10330 

DATA    ((SPETIT  (I,J)  ,J=1,3)  ,1=1,35) /'WHITE    PINE  RED    PINE  ',OTP10343 

1'PITCH    PINE              HEMLOCK                      BALSAM    FIR              «,  OTP10350 

2«TAHARACK                    •,  OTP10360 

3'BD-BL-HH  SPROCEHORHAY  SPROCE   SCOTCH  PINE     CEDAR  ',    OTP10370 

4'OTHER    SOFTWOODSSOGAR    HAPLE  RED    MAPLE  WHITE   ASH  «,  OTP  10380 

5'HILLOB                         ASPEN                            YELLOW    POPLAR       •,  OTP10390 

6«BLACK    CHERRY         WHITE    BIRCH  YELLOW    BIRCH         SWEEP    BIRCH  ',  0TP10400 

7«GREY    BIRCH               BEECH                           BASSWOOD                    ',  OTP10410 

8«RED    OAK                      WHITE    OAK                 BLACK    OAK                 •,  OTP  10420 

9«ELH  GOH  HICKORY  HOPfiOENBEAM  ',         OTP10430 

1 'OTHER    HARDWOODSALL   SOFTWOODS       ',  OTP  10440 

2'ALL    HARDWOODS       ALL    SPECIES            '/  OTP10450 

C    ACCOMULATE    NDMBER    OF    STRATA    SAMPLED  OTP10460 

ISTRA=ISTRA*1  OTP  10470 

C    ACCOMOLATE   TOTAL    NUBBER    OF    POINTS    SAHPLED    IN    ALL    STRATA  OTP10480 

TOTPTS=TOTPTS  +  BPNTS  OTP  10490 

C  OTPIOSOO 

C    CALCOLATE    T- VALUE    FOB   THIS    STRATA  OTP 10510 

C  OTP10520 

TVAL=0  OTP10530 

DF=NPNTS-1  OTP 10540 

IF{DF. GE.  1)       CALL    MDSTI (PROBLV, DF.TVAL)  OTP10550 

IPAGE=IPAGE+1  OTP 10560 

C  OTP10570 

C******************lf  **************  OTP  10580 

C    PRINT    HEADING   ON   OOTPOT    FILE  OTP  10590 

C*********************************  OTP 106 00 

C  OTP10610 

CALL    DATE(IDATE)  OTP10620 

REGOSE='ALL'  OTP10630 

IF(EEGAMS.EQ.'NO')       R EGOSE=' SOME'  OTP  10640 

WEITE(3,100)    TITLE, IPAGE  OTP10650 

OFRAHE=FRAME  OTP 10660 

IF (FRAME.  EQ. 'PLOT')       OFRAHE='    PLOT'  OTP10670 

ALPHA=1.0-PROBLV  OTP10680 

WRITE  (3, 10)     IDATE,SAHBTD, BAF,0 FRAME, NPNTS,CLASIZ, rOPANS,  LOS ANS,  OTP 10690 

1REG0SE,ISTRA,IFILE, ACRES,  ALPHA  OTP10700 

10    FORflAT(//'    INVENT    VEH.    3          11/1/78         I.  N.  E.  R.          O.N.H-'/  0TP10710 

l/'ODATE   OF    RON:    ' , 2 A5//' 0' , A9, '    =' , F6. 2 ,28X, ' NOBBER    OF    ',A5,'S    ',    OTP10720 

2'SAHPLED    =• ,I4/«0DIAMETER   CLASS    SIZE   =',I3,21X  OTP10730 

3, 'HEIGHTS    BEASORED    BY    ',A5,'     HEIGHT    IN    ',A4/  0TP10740 

4«0HEIGHTS    OF    •,A4,'    VOLOBE   TREES    BEASORED ' ,8X ,' STRATUM    #',13/  OTP10750 

5'OINPOT    FILE    NAME:     • ,A1 0,18X, 'ACREAGE    OF   TRACT    =',F7.  1/  OTP10760 

6'OCONFIDENCE   LEVEL   =    '  , F3.  2)  OTP10770 
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C«** ********** **«**«****•***•****«•«» ****««*«**««««  OTP10780 

C    PRINT    BEGBESSION    ANRLYSIS    SOMMAEY  OTP 10790 

Q:^^^^^**^^^t *************************  **************  OT PI  0800 

C  0TP10810 

IP(REGANS.EQ. 'lES')        GOTO    90  OTP10820 

IPiGE=IPAGE*1  OTP10830 

HBITE(3,100)     TITLE, IPAGE  0TP108U0 

BBITE(3,20)  OTP10850 

20  FORaAT('0«,99(«*')/'  *',97X,'*V'  ♦' ,28X,  •  AHALJSIS  OF  THE  HEIGHT' ,  OTP  10860 

1*  -  DBH  EELiriONSHIP',2eX,'*V'  *',97X,«*V'  *' ,35X, «  MODEL:   HEX  •  ,  OTP10870 

2«GHT  =  BO  ♦  B1/DBH',3l»X, '*•/'  *',97X,'*V'  ',  99  (•*•)/' 0«  ,99  ('-«  )  OTP108B0 

3/«  '.ISX,':  •,5X,«B0«  ,<«X,' :•  ,5X,'B1' ,i»X,«:«  ,5X,'R' ,5X,':  STABDAE',  OTP10890 

IJ'D       :       AVEEAGE      :       VOLOHE         :       HEIGHT'/'  SPECIES  :     I NTEBC  •  ,orP1 0930 

5'EPT  :    SLOPE    :   SQOAEED   :   ERROR  %      :    HEIGHT   :    TREES   ',OTP10910 

6'  :    TREES'/'  ',99('-')/'  ')  OTP10920 

IFLAG=0  OTP10930 

DO    50    ISPEC=1,32  OTP1091J0 

IF(TOTNTE(ISPEC)- EQ.O)        GOTO    50  OTP10950 

IF(STB(ISPEC) .LT. 3.AND. NTB(ISPEC)-LT. TOTNTR(ISPEC) )       IFLAG=1  OTP 10960 

IF(NrB(ISPEC) .GE.3)       HBITE(3,30)     (SPETIT (ISPEC, J) , J= 1 , 3)  OTP10970 

1  ,EEGINT(ISPEC) ,B2GC0F(ISPEC) ,aSttE(ISPEC) ,SEBEG(ISPEC)     ,  OTP10980 

2  IJEAHY(ISPEC)  ,TOTNTa(ISPEC),NTE(ISPEC)  OTP10990 
30              FOEMATC    '  ,  3  (A5)  ,  2( ' :  '  ,  F9.  3, 2X)  ,  3  (' :  '  ,F8.  1  ,3X)  ,2  (•  :  '  ,18  ,3  X)  >  OTP11000 

IF(NTR(ISPEC) .LT. 3.AND.BTE(ISPEC) . EQ.  TOrNrB (ISPEC) )  OTP11010 

1  HRITE(3,60)  (SPETIT(ISPEC,J) ,J=1,3) ,SEBEG(ISPEC) ,  0TP11020 

2  flEABY (ISPEC)  , TorNTR(ISPEC)  , NTE (ISPEC)  OTP11030 
IF (BTR( ISPEC) .LT-3.ABD. HTE(ISPEC).LT.  TorNrR(ISPEC)  )  OTP 11 040 

1  HRITE(3,«0)  (SPETIT(ISPEC,J) ,J=1,3) ,SEREG (ISPEC) ,  OTP11050 

2  MEAHY  (ISPEC)  ,TOTHTR(ISPEC)  ,BTB(ISPEC)  OTP11060 
»0      FOEMATC  ',3A5,':',10X,'***  HAENIBG  ***',10X,  0TP11070 

1     2{';',F8. 1,3X),2(  •: ',I8,3X) )  OTP11030 

50  COBTIBOE  OTP  11090 

WEITE(3,250)  OTP11100 

DO  70  ISPEC=33,35  OTP11110 

WRITE (3, 60)   (SPETIT(ISPEC,J),J=1,3) ,SEBEG (ISPEC) , MEAN J (IS PEC(  OTP11120 

1       ,TOTNTR(ISPEC) ,NTE(ISPEC)  OTP11130 

60      FOBHATC  '  , 3A5 ,' : ' ,35X,2  (' : ' ,F8. 1 , 3X) ,2 (' : ' ,18, 3X) )  0TP11140 

70      IF(ISPEC. EQ. 34.0E.ISPEC. EQ. 35)   BBITE(3,250)  OTP11150 

IF(IFLAG.EQ.I)   WEITE(3,80)  OTP11160 

80  FOEHAT('0'/'0***  WAENING  ♦*♦     INSOFFICIEBT  HDMBER  OF  HEIGHT  TB«,0TP11170 

1'EES  HEASOBED  FOE  THIS  SPECIES!')  OTP11180 

C  OTP11190 

C  *****************************************  OTP11200 

C  OOTPOT  STAND  AND  STOCK  TABLES  OTP11210 

C  *****************************************  OTP  11 220 

C  OTP11230 

90  DO  210  ISPEC=1,35  OTP11240 

IF(TEEror (ISPEC). EQ.O. 0)   GOTO  210  OTP11250 

IPAGE=IPAGE+1  0TP11260 

HBITE(3,100)  TITLE, IPAGE  OTP11270 

100    FORMAT('1',14A5,26X,I3/'0')  DTP11280 

C  OTP11290 

C  THE  FOLLOWING  SECTION  COBTBOLS  THE  OOTPOT  OF  THE  SPECIES  FORM  CLASS  OTP11300 

C  PEIBTED  IB  THE  STABD  TABLE  HEADING.   THE  FOEM  CLASS  IS  NOT  PRINTED  OTP  11310 

C  ONDEE  THE  FOLLOWING  CONDITIONS:  OTP11320 

C          1.  STAND  TABLES  FOB  SOFTWOOD, HARDWOOD,  OB  ALL  SPECIES  OTP  11330 

C              (  ISPEC  >=  33  )  OTP113U0 

C          2.  TREES  MEASORED  BY  TOTAL  HEIGHT  (  TOPANS  =  "TOTAL"  )  0TP11350 

C          3.  STAND  TABLES  FOE  SOFTWOODS  WHEN  "MIXED"  HEIGHTS  OSED  OTP11360 

C              (  TOPANS  =  "MIXED"  AND  ISPEC  <=  11  )  OTP11370 

C  0TP11380 

AFcriT='  •  ;  AFCNOn='  •  OTP11390 

IF(ISPEC.GE.33)   GOTO  120  OTP11400 

IF(TOPANS.EQ. 'TOTAL')        GOTO    120  OTP11410 

IF(T0PAHS.EQ. 'MIXED'. AND.ISPECLE.  11)        GOTO    120  OrP11420 

AFCTIT='FC   =    '  OTP  11 430 

C    ENCODE    WEITES    THE    INTEGEB    VALOE    OF    FC    INTO    AFCBOP    AS   ASCII    CHARACTERS    OTP11440 

EBCODE(2,110,AFCN01I)     FC(ISPEC)  0TP11450 

110          FORMAT (12)  OTP11460 

120         WEITE(3,130)      (SPETIT (ISPEC, J) , J= 1 ,3) , AFCTIT, AFCNOM  OTP11470 

130    FOEHAT (•0',99(' *')/'     ','*', 97X, '*' /  OTP11480 

1,'    *',39X,3(A5) ,35X,A5,A2,X,'*'/'    ♦',97X,'*'/  0TP11490 

2'  ',99 (•*')/'0',99('-')/'  ',10X,':   BASAL  ABBA   :',  0TP11500 

3'      TEEES     :       SAWLOG        :       SAWLOG        :',  OTP11510 

4'      PULPWOOD'/'   DIAMETER  : '6X , 'PEE' , 5X , ' : ' ,6X, • PER* , 5X,' : • ,  0TP11520 

5'  CUBIC  FOOT  VOLUME  :  BOARD  FOOT  VOLUME  :  CUBIC  FOOT  VOLUME'/  OTP11530 

6'  ','   CLASS', 3X, ': ',5X,'ACEE',5X,  OTP11540 

7':',5X,'ACRE' ,5X,3(':',5X,'PEE  ACEE',6X)/'  ',99('-')/'  ')  0TP11550 
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DO    150    DBH=CLASIZ,aO,CLRSIZ 

IF{TREPAC(DBH,ISPEC).EQ.O.O)       GOTO    150 

HRITE(3,140)     DBH, BASPAC{DBH, ISPEC) ,TEEPAC (DBH.ISPEC) 
1  ,CBVLSH(DBH,ISPEC) , BDVLSH( DBH, ISPEC) , PLVLSM (DBH, ISP 

1IJ0  FOEMATC     •,i»X,I2,aX,' :•  ,F10.1,'JX,':  ',?10.  1,«Z,' :  SFIO 

1  2(9X,':', F10. 1)) 

150    CONTINOE 

IF  (PLTOT  (ISPEC). EQ.O)        IEEPaL=0 

IF(PLTOT(ISPEC)  .NE.0)IEEPOL=(SEPUL (ISPEC)  ♦TVAL*100/PLTOT (IS 
1.5 
IF  (CBTOT  (ISPEC)  .EQ.O)        IEECtrB=0 

IF {CBTOT( ISPEC)  -HE.O)  IEECOB= (SECUB (ISPEC) *TVAL* 100/CBTOT ( IS 
1.5 
IF(BDTOT(ISPEC) -EQ.O)        IERBBD=0 

IF (BDTOT (ISPEC) -HE. 0) IERBED= (SEBED (ISPEC) ♦TV AL* 100/BDTOT ( IS 
1.5 

WRITE(3,160)     BASTOT (ISPEC)  ,TEETOr (ISPEC) .CBFOT (ISPEC)  .lEECB 
IBDTOT(ISPEC) ,IEEBBD,PLTOT(ISPEC) ,IEEPOL 
160    PO  RMAT('0',  99  ('-M/'O'rax, 'TOTAL', 2X,':',3X,F7.  1,  <»X  ,  •  :  •  ,  3X, 
1,<»X,':',F10.1,'    +',I3,'X',2(3X,':',F10.1,«    +',13,'%')/ 
2'  +  ',52X,'_',2(19X,'_')/'0',99('-')) 
C    CALCOLATE    AND    OOTPOT    HSD    AND   C.  V. 

BSD=SQET( BASTOT (ISPEC) /TEETOT (ISPEC)  /-005U54) 
IBDCV=0    :    IPLCV=0 

IF(BDTOr (ISPEC) -NE.O)        IBDC?=SQET (VAE BR D (ISPEC) )*100/BDTOT( 
1.5 

IF (PLTOT (ISPEC) 
1.5 
HRITE(3,170)     HSD.IBDCV 
170    FOEHAT('0'/'0',5X,«nEAN    STAND    DIAHETEE    =■. 
1,9X,'B0AED    FOOT    COEFFICIENT    OF   VAEIATIOH 
C    CALCULATE    AND    OUTPOT    MERCHANTABLE    USD 
RBSTOT=BASTOT (ISPEC) 
ETRTOT=TRETOT(ISPEC) 
DO    180    DBH=CLASIZ,7,CLASIZ 

RBSTOT=RBSTOT-BASPAC (DBH, ISPEC) 
180         ETflTOT=RrETOT-TREP  AC  (DBH,  ISPEC) 
IF(RTRT0T.LT„0-01)       RHSD=0 

IF(BTETOT.GE.0.01)       EMSD=SQET  (R  BSTOT/BTRTOT/.  0051»5«) 
WEITE(3,190)     RHSD,IPLCV 
190    FOEHAT('0'/'0' ,5X,' HERCHANTABLE    H.S.D.    =    '.FH.I 

1,9X,' POLP-WOOD      COEFFICIENT    OF    VARIATION    =',Ii»,'%'l 
C 
C    CALCOLATE    AND    OOTPOT     (ESTIHATED)     %   CROISE 

IF(FRRHE.EQ.'PLOT')        PCTCHZ=100 ♦( 1 .0-FINPOP) 

IF(  FRAME.  EQ.  'POINT')       PCTCRZ=- 5«5<»*NPNTS*MSD*HSD/ACRES/BAF 
WRITE  (3,200)     PCTCRZ 
200    FORMAT('0'/'0' ,5X, 'PERCENT   CROISE   =',F5-1,'X') 
210    CONTINOE 
C 
C     «*«**4i«i*«««**4i**4[**  «******•  **4i****4i  ************  **4i** 

C  CALCOLATE  AND  OOTPOT  SPECIES  COMPOSITION  PERCENTAGES 

C 

IPAGE=IPAGE+1 
HRITE(3,100)     riTLE,IPAGE 
WRITE(3,220) 
220    FORMAT('0',99(' *')/•    *',97X,'*'/ 

1,'    *',30X, 'SPECIES    COMPOSITION    B?    PERCENT' , 37X, •*• / 

2'    *',97X,'*'/'     ',99('*')/'0',99(«-')/ 

3*     ',3(15X,':'),  2(5X,'SAWLOG',6X,':'),l»X,'P0LP»OOD'/ 


.NE.O)        IPLCV=SQRT (VARPUL(ISPEC) ) *100/PLrOT( 


'  ,  F4  .  1 
=',ia,'X') 


1'     •,3(15X,':') ,3X,«C0BIC   FOOT', 4X , ' : ' , ' 
5,<»X,': ',3X,'CaBIC    FOOT'/ 
6'     ',I»X, 'SPECIES', «X,':','       BASAL    AREA 
73(': ',5X,'V0L0ME«,6X)/'     ',99('-')/'     ') 
IF(CBTOr (35). EQ. 0)       CBTOT (35)=-999 


BOARD    FOOT' 


:',5X, 'TREES 


IF(BDrOT (35)-EQ.O) 
IF(PLTOT(35)-EQ.O) 
DO    21*0    ISPEC=1,3a 
IF  (TEETOT (ISPEC) 


BDTOT  (35)  =-999 
PLTOT (35) =-999 


EQ.0.0)       GOTO    2U0 
BASCOH=BASTOT (ISPEC) ♦lOO/BASTOT (35) 
TRECOM=TRET0r (ISPEC)  • 100/TRETOT (35) 
C0BCOM=CBTOT(ISPEC) •100/CBTOT  (35) 
BDCOM=BDror (ISPEC) ♦100/BDTOT (35) 
P0LCOM=PLTOT (ISPEC) ♦100/PLTOT (35) 


OTP11560 

0TP11570 

OTP11580 

EC)            OTP11590 

.1,           OTP11600 

OTP11610 

OT?11620 

OTP11530 

PEC) )+    OTPIieUO 

OTP 11650 

OTP11660 

PEQ  )  +    0TP11670 

OTP11630 

OTP11690 

PEC))  ♦    OTP  11700 

OTP11710 

B,               DTP  11720 

OTP11730 

F7.1          0TP117U0 

OTP11750 

OTP11760 

0TP11770 

OTP11780 

0TP11790 

ISPEC)  +OTP11800 

OTP11810 

ISPEC) +OTP 11820 

OTP11830 

0TP118«0 

OTP11850 

OTP11860 

0TP11870 

OTP11880 

0TP11890 

OTP11930 

0TP11910 

OTP11920 

DTP11930 

OTP119U0 

OTP11950 

0TP11960 

OTP11970 

0TP11980 

OTP11990 

OTP120D0 

OTP12010 

OTP12020 

OTP  12030 

OTP120IJ0 

OTP12050 

OTP 12060 

OTP12070 

OTP 12080 

OTP12090 

0TP12100 

OTP12110 

OTP12120 

0TP12130 

OTP121«0 

OTP12150 

0TP12160 

OTP12170 

0TP12180 

OTP12190 

OTP12200 

0TP12210 

OTP12220 

0TP12230 

OTP122U0 

OTP12250 

OTP12260 

OTP12270 

0TP12280 

OTP12290 

OTP12300 
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1 


230 
240 
250 


HBITE(3,230)      (SPETIT (ISPEC,J) , J= 1 ,3) , BASCOM         .TRECOH, 

CO  BCOM, BDCOH, POLCOM 
FORHATC     ',3  (AS)  , ':  '  ,2  (3X,F7.  2,  5X,  •  ;  • )  .  2  (UX,  F7.  2,  6X,  •  :  '  ) 
,UX,F7.  2) 

IF{ISPEC.EQ.32)        HRITE(3,250) 
WRITE(3,250) 

FOBOATCO', 99  (•-•)/'     ') 
IF(CBTOT(35) .EQ.-999)       CBTOT(35) =0 
IF(BDTOT(35). EQ.-999)        BDTOT(35)=0 
IF(PLTOT{35).EQ--999)        PI.TOT(35)  =0 

C 

C  ««***««***«'*•  *****4'<)>«**  *************************  ********** 

C  OUTPUT  VOLUHE  TOTALS  FOR  ALL  SPECIES 

C*** ***************************************** ************** 

C 

IPAGE=IPAGE+1 
WRITE  (3,100)     TITLE.IPAGE 
«RITE(3,260J 
260    FORHAT(«0', 99  {•*•)/•    *',97X,'*V 

1'    *«,30X, 'VOLUHE    TOTALS    FOR    ALL   S PECIES • , 38X, • *•/ 
2'    *',97X,'*V'     '  ,99('*«)/'0' ,99  ('-•)/ 
3«    • ,15X,2(' :«,11X,'SAHL0G',11X)  ,•: ', 10X, • PULP  WOOD'/ 
Ui  SPECIES  : sex, 'CUBIC    FOOT    VOLUHE • ,5X, •:', 6X, • BOARD' , 

5'    FOOT    V0L[JBE',5X,':' ,6X, 'CUBIC    FOOT    VOLUHE'/ 
6'     ',15X,2(': •,10X,'PER    ACRE', 10X) ,'C',10X, 'PER    ACRE'/ 
7'    ',99('-')/'    ') 
DO    290    ISPEC=1,35 


EQ. 0. AMD. ISPEC. LT.  33) 


IF  (CBTOT (ISPEC) .EQ.O.AND. PLTOT (ISPEC) 
1         GOTO  280 

HRITE(3,270)      (SPETIT (ISPEC, J(  , J= 1 , J)  ,CBTOT{ISPEC)  ,BDTOT (ISPEC) , 
1  PLTOT(ISPEC) 

270    FORHATC     '  ,3  (A5)  ,  3( ': ',  F17.  1,  1 1X)  ) 

280         IF(ISPEC.EQ.32.0R.ISPEC.EQ-31».OR.ISPEC.  EQ.35)       HBITE(3,250) 
290    CONTINUE 
C 
C    ***************************************************** 

C  OUTPUT  VOLUME  TOTALS  EXPANDED  BY  ACREAGE 

C 

IPAGE=IPAGE+1 
WBITE(3,100)    TITLE.IPAGE 
HRITE(3,300) 
300    FOBHAT('0',99('*')/'    *',97X,'*'/ 

1'    *' ,30X,'VOLOHE   TOTALS    EXPANDED    BY    ACREAGE' ,3«X,  •  *•/ 
2'    *',97X,'*'/«     ',99('*«)/'0', 99 ('-')/ 
3*     ',15X,2(«:',11X,'SAHLOG',11X) ,': ' , 10X , ' PULPWOOD  '  / 
«•  SPECIES  :',6X, 'CUBIC    FOOT    VOLUME ', 5X ,':' ,6X, • BOARD' , 

5'    FOOT    VOLUME', 5X, •:' ,6X, 'CUBIC    FOOT    VOLUHE'/ 
6'     ',99('-')/'    ') 
DO    330   ISPEC=1,35 

IF  (CBTOT (ISPEC) . EQ.O.AND. PLTOT (ISPEC)  ,  EQ. 0. AND. ISPEC. LT. 33) 
1         GOTO  320 

TOTACB=CBTOT(ISPEC) *ACRES+.5 
TOTABD=BDTOT (ISPEQ  *ACRES*.5 
TOTAPL=PLT0T (ISPEC) *ACRES+.5 

WRITE(3,310)      (SPETIT (ISPEC, J)  ,J=1,3) , TOTACB ,T0TABD,T0TAPL 
310         FORHATC    '  ,3  ( A5)  ,3  ( «:  '  ,1 17, 1 1  X)  ) 

320         IF(ISPEC.  EQ.32.OR.ISPEC.EQ.3U.0R.ISPEC.EQ.35)        BBITE(3,250| 
330    CONTINUE 
RETURN 


OTP12310 
OTP12320 
OTP  12330 
OTP123U0 
OTP12350 
OTP12360 
OTP12370 
OTP123B0 
0TP12390 
OTP12'»00 
OTP12mO 
OTP12II20 
OTP12430 
0TP12HH0 
OTP12l*50 
OTP  121160 
OTP12U70 
0TP12tt80 
OTP12l»90 
OTP12500 
0TP12510 
OTP12520 
OTP  12530 
OTP125a0 
OTP12550 
OTP  12560 
OTP12570 
OTP 12580 
OTP12590 
OTP12600 
0TP12610 
OTP12620 
OTP 12630 
OTP12640 
OTP12650 
OTP  12660 
OTP12670 
OTP 12680 
OTP12690 
0TP12700 
OTP 127 10 
OTP12720 
OTP  12730 
OTP127aO 
3TP12750 
OTP  12760 
Orpi2770 
OTP 12780 
OTP12790 
OTP12800 
0TP12810 
OTP12820 
OTP12830 
OTP12840 
OTP 12850 
OTP12860 
OTP12870 
OTP12880 
OTP12890 
OTP 12900 
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Q**^t* ************************  ******************************************* 

c 

C         OOTPOT    TABLES    FOE    STRATIFIED   TOTAL 

C 

C************ *******************  **************************************** 

C 

ENTRY    O0rP2 
WRITE(5,3«0) 
3U0    FOEHAT(«0HHAT    DO    100    WANT    TO    TITLE   THE   COHBINED    STRATA    OOTPOT?') 

BEAD(5,350)     TITLE 
350    F0EHAT(1<JA5) 
IPAGE=IPAGE+1 
WRITE(3,100)     TITLE, IPAGE 
C    PRINT    HEADING    ON   OOTPOT    FILE 
OFRAME=FRAaE 

IF(FRAHE. EQ. 'PLOT')       OFBAHE=»    PLOT* 

WRITE (3,360)    IDATE.CLASIZ, ISTEA,TACRES,OFRAHE ,TOTPTS, ALPS  A 
360    FORaAT(//'    INVENT    VER.    3         11/1/78         I. N. E. R.  O.N.H. •// 

VODATE   OF    RON:    '.aAS// 

2«0DIAMETER    CLASS    SIZE    = • ,13 ,21X , • NOHBEB   OF    STRATA    SAMPLED    =• 
3,I3/«0TOTAL    ACEEAGE   SAHPLED   =• ,F9. 1 , 13X,« TOTAL    •,A5 
l»,'S    SAMPLED    =«,I5/«0CONFIDENCE    LEVEL    =    ',F3.2) 

C 

C**** ********************************** 

C    OOTPOT    DIAMETER    CLASS    VOLOHE    TABLES 

c 

DO    U50   ISPEC=1,35 

IF(CBTST{ISPEC)  .EQ.O.  AHD.PLTST(ISPEC)  .EQ.O)       GOTO    1*50 
IPAGE=IPAGE+1 
WRITE(3,100)     TITLE, IPAGE 
C 

C    THE    FOLLOWING   SECTION    CONTROLS    FOEB    CLASS    PRINTING    AS    PREVIOOSLY 
C    DESCRIBED.       THE    INPOTS    OF    THE    LAST    STRATOH    ANALYSED   CONTROLS    WHAT 
C    IS    PRINTED    IN    THE   STRATIFIED   TOTAL    OOTPOT. 
C 

AFCTIT='     •     :     AFCNOa=«     • 
IF(ISPEC.GE.  33)        GOTO    370 
IF(TOPANS.EQ. "TOTAL')       GOTO    370 

IF(TOPANS.EQ."HIXED'.  AND.ISPEC.LE.  11)       GOTO    370 
AFCTIT=«FC   =     ' 
C    ENCODE    WRITES    THE    INTEGER    VALOE    OF    FC    INTO    AFCNOM    AS    ASCII    CHARACTERS 
ENCODE(2,110,  AFCNOH)     FC(ISPEC) 
370         WRITE(3,380)      (SPETIT(ISPEC, J) ,J=1 ,3) , AFCTIT, AFCNOH 
380         FORMAT('0', 99  {•*•)/'    *',97X,'*'/ 
*',39X,3(A5) ,35X,A5,A2,X,'*'/ 
•  •,97X,'»'/'     ',  99  C*')/'©',  99  {•-')/ 

DI AHETEH'  ,  1»X,2  ( •  :  •  ,  1 1X,  'SAWLDG  •  ,  1 1X)  ,  •  :  •  ,  1  OX  ,  '  POLPWOOD'  / 
CLASS  i',6X,'C0BIC    FOOT    VOLOME • ,5X, ':' ,6X,  ' BOARD' , 

FOOT    VOLOME', 5X, ':• ,6X,'C0BIC    FOOT    VOLOHE'/ 
',99('-')/'    •) 
DO    UOO    DBH=CLaSIZ,UO,CLASIZ 

IF(CBVLST (DBH, ISPEC) . EQ- 0. AND- PLVLST(DBH, ISPEC) .EQ.O) 

GOTO    400 
ICBVST=CBVLST (DBH, ISPEC)  +. 5 
IBDVST=BDVLST(DBH,ISPEC)+.5 
IPLVST=PLVLST  (DBH,  ISPEC)  +.5 
WRITE  (3,390)     DBH,ICBVST,IBDVST,IPLVST 
390  FORHATC    ',I8,7X,3(':',I17,  11X)) 

HOO         CONTINOE 

ICBrST=CBTST (ISPEC) ♦. 5 

IPLTST=PLTST (ISPEC) +. 5 

C    CALCOLATE    ERROR    AS    A    PERCENT 

IERCOB=0;IERBBD=0:IEBPOL=0;TVAL=0 

IF(CBTST(ISPEC)-EQ,0)       GOTO      410 

IF(EFDFCB (ISPEC) .GE-1)CALL    MDSTI (PEOBLV,EFDFCB(ISPEC) ,TVAL) 

IERCOB= (STSECB(ISPEC) *TVAL*100*TACEES/CBTST (ISPEC) ) *.S 

IF(BDTST(ISPEC) .EQ.O)       GOTO   420 

TVAL=0 

IF(EFDFBD (ISPEC) . GE- 1 ) CALL    MDSTI (PEOBLV, EFDFBD(ISPEC) , TVAL) 

IERBRD= (STSEBD(ISPEC) *TVAL* 100*TACEES/BDTST (ISPEC) )  +.5 

IF  (PLTST (ISPEC). EQ.O)       GOTO    430 

TVAL=0 

IF(EFDFPL(ISPEC).6E.  1)  CALL  MDSTI  (PEOBLV,  EFDFPL(  ISPEC)  ,TVAL) 

IERPOL= (STSEPL (ISPEC) *TVAL*100*TACRES/PLrST (ISPEC) ) +. 5 

WRITE (3,440)     ICBTST, lERCOB ,IBDTSr,IEBBRD,IPLTSr,IERPOL 

FORM AT ('0', 99 ('-' ) /'O ', 4X, 'TOTAL ',6X 

1  ,3(':',I17,'    +',I3,'X',5X)/'+*,15X,3(19X,'_',9X)/ 

2  '0',99('-')) 


1* 
2' 
3' 
U' 
5' 
6' 


1 


IBDTST=BDTST  (ISPEC) +. 5 


410 


420 


430 

440 


450  CONTINOE 


OTP12910 
DTP  12920 
OTP12930 
OTP12940 
OTP 12950 
OTP12960 
OTP 12970 
OTP12930 
OTP 12990 
OTP13000 
OTP13010 
OTP13020 
OTP13030 
OTP  13040 
OTP13050 
OTP13060 
OTP  13070 
OTP13080 
OTP13090 
OTP13100 
OTP13110 
0TP13120 
OTP13130 
0TP13140 
0TP13 150 
OTP13160 
0TP13170 
OTP13iaO 
0TP13190 
OTP13200 
0TP13210 
OTP13220 
OTP 13230 
OTP13240 
OTP13250 
OTP 13260 
OTP13270 
OTP  13280 
OTP13290 
OTP13300 
OTP13310 
OTP13320 
0TP13330 
OTP13340 
OTP13350 
OTP13360 
OTP13370 
OTP  13380 
OTP13390 
OTP13400 
OTP  13410 
OTP13420 
OTP 13430 
OTP13440 
OTP 13450 
OTP13460 
OTP 13470 
OTP13480 
OTP 13490 
OTP13500 
OTP13510 
OTP 13520 
OTP13530 
OTP  13540 
OTP13550 
OTP13560 
0TP13570 
OTP13580 
OTP  13590 
OTP13600 
OTP13610 
OTP13620 
OTP13630 
OTP 13640 
OTP13650 
OTP13660 
OTP  13670 
OTP13680 
OTP  13690 
OTP13700 
OTP13710 
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C  OTP 13720 

C***!*******  <■*************  •****«*«**  OTP13730 

C  OOTPOT  &LL  SPECIES  VOLOME  TABLE  OTP  13740 

C***«*«*4'***4<****  **********  ********  OT  PI  3750 

C  OTP  13760 

IPAGE=IPAGE  +  1  OTP13770 

WRITE (3, 100)  TITLE, IPAGE  OTP13780 

(IBITE(3,<»60)  OTP13790 

460    FORHATCO', 99  (•*')/•    •',97X,'*V  OTP13830 

1'    *',30X,'VOL0HE    TOTALS    FOE    ALL   S PECIES • ,38X, • *  V  OTP13810 

2«    *',97X,«*V'     •,99('*«)/'0',99('-')/  OTP13820 

3«     •,15X,2  Cx'.IIX.ISAHLOG'.IIX)  ,•:  •,10X,«POLPHOGDV  OTP13830 

«•  SPECIES  :«,6X,'C0BIC    FOOT    VOLO HE •, 5X ,«:•, 6X, • BOARD    ',  OTP138U0 

5'FOOT    VOLOHE'.SX, 'J'.ex, 'COBIC    FOOT    VOLOBE'/  OTP13850 

6'     ',99{'-')/'    •)  OTP13860 

DO  480  ISPEC=1,35  OTP13870 

IF(CBTST(ISPEC) . EQ. 0. AND. PLTST (ISPEC) . EQ- 0. AND.  ISPEC. Lr.33)  OTP  13880 

1         GOTO  470  OTP13890 

ICBTST=CBTST(ISPEC) +.5  ;  IBDTST=BDTST (ISPEC) +.5  OTP139aO 

IPLTST=PLTST(ISPEC) +. 5  OTP  13910 

HRITE(3,310)   (SPETIT(ISPEC,J)  ,J=1,3)  , ICBTST.IBDTST.IPLTST  OTP13920 

470    IF(ISPEC. EQ. 32. OR. ISPEC. EQ.34. OR. ISPEC. EQ.35)   HRirE(3,250)  OTP13930 

480  CONTINUE  OTP  13940 

RETURN  OTP13950 

END  OTP 13960 
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C  HDST0010 

SOBPOOTINE    HDSTI     (Q,F,X)  HDST0020 

C  f1DST0030 

C »•*»•*»*••*****•*»•»*♦*»********»***************»♦********♦**»*»***♦♦♦♦ a DS TO  OHO 

C  nDSTOOSO 

C  THIS    SOBPROGBRM    IS    COMPOSED    OF    TWO    SOBROOTINES    ADAPTED    FROH    THE  (1DST0060 

C  INTERNATIONAL    HATHEWATICAL    AND    STATISTICAL    LIBRARIES , INC.      (IBSL)  nDST0070 

C                 HDSTI      -       INVERSE   STODENT'S    T    DISTRIBOTION  HDST0080 

C                 BDNRIS    -       INVERSE    NORHAL    PROBABILITY    DISTRIBOTION    FUNCTION  HDSr0090 

C  BOTH    SUBROUTINES    HAVE    BEEN    EDITED    REHOVING    CAPABILITIES    NOT    NEEDED    nDSTOlOO 

C             IN    THIS    APPLICATION.        THE    PEHAINING    CODE    IS     REPRODUCED    HERE    IN  nDSTOIIO 

C            ACCORDANCE    WITH    THE    IMSL    POLICY    ON    RESEARCH    WORK    AS    STATED    3N  nDST0120 

C            PAGE    INTRO-22    OF    THE   JULY,     1977    LIBRARY    HANOAL.  HDST0130 

C  HDSTOmO 

Qm*m***********************************  ****************************  *****fiQST:0  150 

C  SDSTOieO 

C                                                                                            EXACT    INTEGRAL    FOP    2    D.F.  HDST0170 

IF     (ABS(F-2.0)     .GT.    .000001)        GO    TO    10  HDST0180 

X    =    SQET(2.0/(Q*(2.0-Q)  )-2.0)  (1DST0190 

GO    TO    50  nDST0200 

10    HPI    =    1.570796  MDST0210 

C                                                                                            EXACT    INTEGRAL    FOE    1    D.F.  nDST0220 

IF     (ABS(F-I.O)     .GT.     .000001)        GO   TO    20  HDST0230 

A    =    Q*HPI  HDST02U0 

X    =    COS(A)/SIN(A)  nDST0250 

GO    TO    50  HDST0260 

C                                                                                            EXPANSION    FOR    N    GREATER    THAN    2  BDST0270 

20    A    =    1.0/(F-0.5)  nDST0280 

B    =    «b.O/(A*A)  MDST0290 

C    =    (  (20700. ♦A/B-98.) 'A-ie.) *A  +  96. 36  BDST0300 

D    =    ((9«.5/(B+C)-3.0) /B+1.0) •SQRT(A*HPI) ♦F  MDST0310 

XX    =    D*Q  nDST0320 

Y  =  XX«*(2.0/F)  (1DST0330 
IF     (Y    .GT.    A*. 05)       GO    TO    «0  BDST03H0 

Y  =    { (1.0/( ((F+6.0)/(r*Y)-0.0«9*D-0.822)* (F+2.0)»3.0) ♦  BDST0350 
1            0.5/(F+i».0)  )  »Y-1.0)  »(F+1.0)/(F  +  2.0) -H.O/Y  MDST0360 

30    X    =    SQRT(F*Y)  HDST0  370 

GO  TO  50  BDST0380 

C  ASYHPTOTIC  INVERSE  EXPANSION  ABOUT    BDSr0390 

C                                      HOBHAL  BDSTOaOO 

<»0  X  =  .5  *  Q  BDSTOHIO 

CALL    BDNRIS(X,XX)  HDST0U2O 

Y  =  XX  *  XX  BDST0U30 
IF  (F  .LT.  5.)  C  =  C+0.3*{F-«.5) ♦(XX+0.6)  HDSTOUHO 
C=  ({{.05*D*XX-5.0)»XX-7.0) *XX-2.0J*XX+B+C  BDST0U50 
f    =     (((((0.«»Y  +  6.  3)  *Y*36.)  *Y  +  9i«.5) /C-Y-3.0) /B  +  1.0)  •XX  HDST0«60 

Y  =  A»Y*Y  HDST0l»70 
D    =    Y  BDST0U80 

Y  =  .05*Y*Y+Y  BDST0H90 
IF  (Y  .GT.  .002)  Y  =  EXP{D)  -  1.0  BDSTOSOO 
GO    TO    30  BDST0510 

50    RETURN  BDST0520 

END  MDST0530 
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c 
c 


SOBBOOTINE   flDNR 

DATA 

DATA 

DATA 
1 

DATA 
1 

DATA 

DATA 
1 

DATA 

DATA 
1 

DATA 

GO    TO    40 
10    IF     (X.LT.    .  15) 

IP     {X.LT.  1.85) 

SIGHA   =    -1. 

A    =    2.-X 

B   =    X-1. 

GO    TO    50 
20    Z    =    U-X 

SIGHA    =    SIGN(1. 

Z    =   ABS{Z) 

GO    TO   80 
30    SIGHA    =    1. 

A    =    X 

B    =    1<.-X 

GO    TO    50 

no    INT   =    3 
X    =    2-  *P 
GO    TO    10 


50    W    =    SQET(-ALOG (A+A*B) ) 
IF     (W.  LT.  2.5)        GO    TO    70 
IF     (W.LT-  n.  )       GO    TO    60 


IS(P,Y) 

SQBT2/1.  I»1«2136/ 

A 1,A2,A3/-. 575 17029, -1.89651 33, -.0  5496 2605/ 

BO, B1,B2,B3/-. 11 377303, -3. 2934740, -2. 3749959, 

-1.1875145/ 
CO, C1,C2,C3/-. 11 466659, -.13147744, -.236820 10, 

-050739749/ 
DO, D1,D2/-44. 279769, 21.985462,-7-5861027/ 
E0,E1,E2,E3/-. 0566842208,. 39370209, -.3 16650 10, 

.062089629/ 
F0,F1,F2/-6. 2667859,4. 6662627,-2-96  28832/ 
GO, G1,G2,G3/. 000 1851 1591, -.0020281 520, 

-.  14983844,. 010786386/ 
HO, HI, H2/. 099529751, .521 17329,-. 068883  009/ 

GO  TO    30 
GO    TO    20 


-Z) 


WI    =    1./H 
SS   =    ((G3*iri+G2 
SD    =     ((HI-t-H2)*H 
F    =    »    ■<■    H«(GO+S 
GO   TO   90 


)*WI+G1)«BI 
I+H1) ♦HI+HO 
N/SD) 


60  SB  =  ((E3«W+E2) *H*E1) *H 
SD  =  ((H  +  F2)*H  +  F1)*ll*F0 
F  =  H  +  H»(E0+SN/SD) 
GO  TO  90 


70  SH  =  {(C3»8*C2)  *B*C1)  ♦« 
SD  =  ((B+D2)*H+D1) *W+DO 
F  =  H  ♦  B»(CO+SH/SD) 
GO  TO  90 


c 

c 

80 

Z2    =    Z*Z 

F   =    Z*Z*{B0*A1« 

c 

c 

90 

Y  =    SIGHA*F 
IF     (INT.NE.  3) 

Y  =    -Y*SQBT2 

100 

BETOEH 
ESO 

INVEBSE  GAOSSIAN  ENTBY 


REDOCED  ABGOHENT  IS  IN  (.85,1.), 
OBTAIN  THE  TBANSFOEHED  VABIABLE 


W  GBEATER  THAN  4.,  APPEOX.  F  BY  A 
BATIONAL  FUNCTION  IN  1./B 


B  BETBEEN  2.5  AND  4.,  APPEOX.  F 
BY  A  RATIONAL  FONCTION  IN  B 


B    BETBEEN    1-13222    AND    2.5,    APPEOX. 
F    BY    A    BATIONAL    FaNCTION    IN    W 


Z  BETBEEN  0.  AND  -85,  APPBOX.  F 
BY  A  EATION&L  FONCTION  IN  Z 


Z2/(B1*Z2+A2/(B2+Z2-fA3/(B3  +  Z2)  )  )  ) 
FOEH   THE    SOLUTION    BY 
THE    PROPEE    SIGN 

GO    TO    100 


MOLT.    F    BY 


HDNE0010 
HDNB0020 
HDNB0030 
HDHB0040 
HDHR0050 
HDNE0060 
HDNE0070 
HDNEOOBO 
HDNE0090 
HDNR0100 
HDNE0110 
nDNE0120 
HDHBOnO 
HDNE0140 
HDHB0150 
HDNB0160 
SDNR0170 
HDNB01B0 
HDNE0190 
nDNB0200 
HDNB0210 
HDNB0220 
HDNE0230 
tlDNR0240 
HDNE0250 
HDNB0260 
HDNE0270 
HDNE02B0 
HDNE0290 
HDNB0300 
aDNB0310 
HDNB0320 
MDNB0330 
HDNB0340 
HDNR0350 
aDNB0360 
HDNS0370 
HDNB0330 
HDNE0390 
HDNB0400 
HDNB0410 
nDNB0420 
HDNE0430 
{1DNR0440 
flDNR0450 
HDNR0460 
BDNE0470 
aDNB0480 
HDNE0490 
HDNR0530 
{1DNB0510 
ODNE0520 
HDNE0530 
HDKB0543 
HDNB0550 
HDNB0560 
HDNE0570 
HDMB0580 
HDNR0590 
aDNB0600 
nDNE0610 
HDNE0620 
flDNB0630 
HDNE0640 
HDNS0650 
HDNE0660 
HDNE0670 
{IDNE0680 
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